2014.12.27训练笔记

快速幂算法

7^13

指数部分：(13)10 == 1101

1　　　　2　　　　4 　　　　8 　　　　16 ……
7^1　　7^2　　 7^4　　 　7^8 　　　7^16 ……

支持同余定理

然后，二进制每位询问（用位于运算符&），是1的，相乘：
比如上面的可以转化为

7^1 * 7^4 * 7^8 = 7^13

求7^132 mod 47

int a = 7, b = 132;
int rel = 1, cnt = 1;
while(b){
    if(b & 1){
        rel = (rel * a) % 47;
    }
    a *= a;
    a %= 47;
    b >>= 1;
}    

 

快速求素数表（去掉合数）

#include <iostream>
#include <cmath>
#include <string.h>

#define SIZE 2000000

using namespace std;

int pr[SIZE];

int main()
{

    int i = 0, j = 0;
    memset(pr, 0 ,sizeof(pr));
    int n = SIZE -1;
    for(i = 2; i <=n;i++ ){
        if(!pr[i]){
            for(j = i + i; j <= n; j+=i){
                pr[j] = 1;
            }
        }
    }
    cout << "done"<<endl;
    return 0;
}

 

文件读写

函数原型：freopen(const char* filesource, const char* mode, FILE* stream);

stream 可以是stdin , stdout;

 

多边形面积公式

S=(Σ（i=1， i <= n) (x[i] * y[i+1] - x[i+1] * y[i]) /2 ) + (x[n]*y[1] - x[1]*y[n]) / 2
S = abs(S);