同余与模运算

模运算:

(a+b)modn = ((amodn)+(bmodn))modn

(a-b)modn = ((amodn)-(bmodn)+n)modn

(a*b)modn = ((amodn)*(bmodn))modn

int mul_mod(int a,int b,int n)
{
    a%=n;  b%=n;
    return (int)((long long) a*b % n);
}

例1:大整数取模nmodm,你<10^100,m<10^9

1234 = ((1*10+2)*10+3)*10+4;

依次取模即可

1 int ans = 0;
2 for(int i=0;i<len;i++)
3 {
4     ans = ((long long)ans*10 +n[i])%m;
5 }


例2幂取模

快速幂

1 int pow_mod(int p,int n,int m)
2 {
3     if(n==1)
4         return p;
5     int x=pow_mod(p,n/2,m);
6     long long ans = (long long)x*x%m;
7     if(n%2)ans*=p;
8     return (int)(ans % m);
9 }

 

例3:模线性方程

输入a,b,n,解方程组a*x≡b(modn)a,b,n<10^9

对于同余符号,可以这么理解,就是a*x = y1*n+K

                                             b =y2*n+K;

所以ax-b = (y1-y2)*n。令其为ax-b = y*n,方程转化为ax-ny=b,也就是求线性方程的一个x的解

由欧几里得公式,如果b%gcd(a,n)=0时方程有解,否则无解

 

 

posted @ 2013-07-22 13:33  再见~雨泉  阅读(547)  评论(0编辑  收藏  举报