大数a^b%n模板
2014-04-24 14:35
190 查看
a^b%n用到两个技巧。
首先是幂方,类似于二分。
其次是同余定理: (a * b)%n = (a%n * b%n)%n
加法也是可以用的。
吉林大学有个二进制版本的大数a^b%n的:
// a^b0 mod n
int mod_exp(int a,int b0,int n)
{
if( a>n ) a%=n;
int i,d=1,b[35];
for( i=0;i<35;++i )
{
b[i]=b0%2;
b0 /= 2;
if( b0==0 ) break;
}//b[i]b[i-1]...b[0] 为b0的二进制表示
for( ;i>=0;--i )
{
d= (d*d)%n;
if( b[i]==1 ) d=(d*a)%n;
}
return d;
}
这个是用int的,一般不用long long 因为%n把过大的数都取模了
int mod_exp(int a,int b,int n)
{
int ans=1;
while(b)
{
if(b&1)
ans = (ans * a) % n;
a = a * a % n;
b = b >> 1;
}
return ans;
}
首先是幂方,类似于二分。
其次是同余定理: (a * b)%n = (a%n * b%n)%n
加法也是可以用的。
吉林大学有个二进制版本的大数a^b%n的:
// a^b0 mod n
int mod_exp(int a,int b0,int n)
{
if( a>n ) a%=n;
int i,d=1,b[35];
for( i=0;i<35;++i )
{
b[i]=b0%2;
b0 /= 2;
if( b0==0 ) break;
}//b[i]b[i-1]...b[0] 为b0的二进制表示
for( ;i>=0;--i )
{
d= (d*d)%n;
if( b[i]==1 ) d=(d*a)%n;
}
return d;
}
这个是用int的,一般不用long long 因为%n把过大的数都取模了
int mod_exp(int a,int b,int n)
{
int ans=1;
while(b)
{
if(b&1)
ans = (ans * a) % n;
a = a * a % n;
b = b >> 1;
}
return ans;
}
相关文章推荐
- HDU 5047 Sawtooth (大数模板+找规律)
- C++大数模板
- 大数乘法模板
- Java 大数模板
- HDU 1041 递推 + 大数模板
- 大数运算模板
- 我的挫13大数模板
- 大数高精度加减乘除模板
- 高精度模板(含大数开平方)
- 大数模板
- 大数高精度运算(模板)
- HDU 6144 Arithmetic of Bomb(字符串+大数取余模板)
- 杭电OJ 1042 N!(大数阶乘模板)
- 大数除法模板(c语言)
- NEFU 748 A+B+C大数版 大数模板题
- A + B Problem II 大数相加(可以作为模板题)
- 大数乘法模板
- c++大数模板
- 用数组模拟大数加法模板
- HDU 1402 FFT模板题,求大数乘法