數(shù)論
最大公約數(shù)(GCD)/最小公倍數(shù)(LCM)
/*非遞歸版本求最大公約數(shù)*/
int gcd(int a,int b)
{
if(0==b) return a;
while(b>0)
{
int temp=a%b;
a=b;
b=temp;
}
return a;
//while(b^=a^=b^=a%=b)
}
/*遞歸版本求最大公約數(shù)*/
int gcd(int a,int b)
{
return b?gcd(b,a%b):a;
}
/*求最小公倍數(shù)*/
int lcm(int a,int b)
{
return a*b/gcd(a,b);
}
素?cái)?shù)判斷及打表
/*判斷n為素?cái)?shù)返回1,合數(shù)返回0*/
int IsPrime(int n)
{
if(n==2) return 1;
if(n%2==0||n<2) return 0;
int l=sqrt(n+1);
for(int i=3;i<=l;i+=2)
if(n%i==0) return 0;
return 1;
}
/*素?cái)?shù)表打表 歐拉篩法*/
/*primes[] [2,N] 之間的素?cái)?shù) primes[0]第0位素?cái)?shù)2*/
/*isprime[] */
const int maxn=1000007;
int primes[maxn];
int isprime[maxn];
void euler_sieve()
{
int tot=0;
memset(isprime,1,sizeof(isprime));
isprime[0]=isprime[1]=0;
for(int i=2;i<=maxn;i++)
{
if(isprime[i]) primes[tot++]=i;
for(int j=0;j<tot;j++)
{
if(i*primes[j]>maxn) break;
isprime[i*primes[j]]=0;
if(i%primes[j]==0) break;
}
}
}
快速冪/乘取模
/*快速冪取模 (a^b)%p*/
int quickpow(int a,int b,int p)
{
int ret=1;
a%=p;
while(b)
{
if(b&1) ret=(a*ret)%p;
a=(a*a)%p;
b>>=1;
}
return ret;
}
/*快速乘法取模 (a*b)%p*/
int quickmul(int a,int b,int p)
{
int ans=0;
while(b)
{
if(b&1) {b--;ans=(ans+a)%p;}
b>>=1;
a=(a+a)%p;
}
return ans;
}
拓展歐幾里得
/*拓展歐幾里德 ax+by=1*/
int exgcd(int a,int b,int &x,int &y)
{
if(b==0)
{
x=1;
y=0;
return a;
}
int ans=exgcd(b,a%b,y,x);
y-=a/b*x;
return ans;
}
求乘法逆元
/*求乘法逆元(數(shù)論倒數(shù)) ax=1(modb) -> ax-1=by -> ax+by=1*/
int modinverse(int a,int b)
{
int x,y;
int d=exgcd(a,b,x,y);
return d==1?(x%b+b)%b:-1;
}
中國(guó)剩余定理(解一元線性同余方程組)
/*中國(guó)剩余定理(互質(zhì)) 一元線性同余方程組 x=a(modm)*/
/*設(shè)ai,Mi為除該數(shù)外的摸數(shù)乘積,ti為Mi模mi數(shù)論倒數(shù)*/
/*ans(ai*ti*Mi)*/
/*a[],p[]模,n,長(zhǎng)度*/
int crt(int a[],int p[],int n)
{
int muls=1;
int ret=0;
for(int i=0;i<n;i++) muls*=p[i];
for(int i=0;i<n;i++)
{
int x,y;
int mi=muls/p[i];
exgcd(mi,p[i],x,y);
ret=(ret+x*a[i]*mi)%muls;
}
return (ret+muls)%muls;
}
/*非互質(zhì),無(wú)解返回-1*/
int crt(int a[],int p[],int n)
{
if(n==1)
{
return p[0]>a[0]?a[0]:-1;
}
int x,y,d;
for(int i=1;i<n;i++)
{
if(p[i]<=a[i]) return -1;
d=exgcd(p[0],p[i],x,y);
if((a[i]-a[0])%d!=0) return -1;
int t=p[i]/d;
x=((a[i]-a[0])/d*x%t+t)%t;
a[0]=x*p[0]+a[0];
p[0]=p[0]*p[i]/d;
a[0]=(a[0]%p[0]+p[0])%p[0];
}
return a[0];
}
求解方程ax=b(modn)
/*求解ax=b(modn), 解的個(gè)數(shù)為gcd(a,n)，返回x為vector<>*/
verctor<int> lmodeq(int a,int b,int n)
{
int x,y;
int d=exgcd(a,n,x,y);
vector<int> ans;
ans.clear();
if(b%d==0)
{
x=(x%n+n)%n;
x%=(n/d);
ans.push_back(x*(b/d)%(n/d));
for(int i=1;i<d;i++)
ans.push_back((ans[0]+i*n/d)%n);
}
return ans;
}
莫比烏斯函數(shù)求解
/*求莫比烏斯函數(shù)mu
mu={
1 ,u=1
(-1)^k ,能分解成k個(gè)不同的質(zhì)因數(shù)因子
0 ,u包含平方因子
}
*/
int Mobius(int n)
{
int cnt,k=0;
for(int i=2;i*i<=n;i++)
{
if(n%i) continue;
cnt=0;
k++;
while(n%i==0)
{
n/=i;
cnt++;
}
if(cnt>=2) return 0;
}
if(n) k++;
return k%2?-1:1;
}
組合數(shù)求解(逆元求組合數(shù)/Lucas定理)
/*求組合數(shù):數(shù)較小mod較大時(shí),使用逆元.數(shù)較大mod較小時(shí)使用lucas*/
/*逆元求組合數(shù)*/
int exgcd(int a,int b,int &x,int &y)
{
if(b==0)
{
x=1;
y=0;
return a;
}
int ans=exgcd(b,a%b,y,x);
y-=a/b*x;
return ans;
}
int fac(int n,int p)
{
int sum=1;
for(int i=1;i<=n;i++)
{
sum=(num*i)%p;
}
return sum;
}
int comb(int n,int m,int p)
{
int a=fac(m)*fac(n-m)%p;
int x,y;
exgcd(a,p,x,y);
return ((fac(n)*x)%p+p)%p;
}
/*Lucas定理求組合數(shù)*/
int Lucas(int n,int m,int p)
{
if(m==0) return 1;
return comb(n%p,m%p,p)*Lucas(n/p,m/p,p)%p;
}
/*組合數(shù)打表 c(n,m)*/
int initcomb()
{
for(int i=0;i<MAX;i++)
{
Co[i][0]=Co[i][i]=1;
for(int j=1;j<i;j++)
{
Co[i][j]=(Co[i-1][j]+Co[i-1][j-1])%p;
}
}
}
