常用模板

快速幂取模 求(a^b)%c 
int q_mod(int a,int b,int c){
	int ans=1;
	a=a%c;  //先求余，缩短运算规模。(有时候不用)
	while(b){
		if(b&1)
		ans=(ans*a)%c;
		b>>=1;
		a=(a*a)%c;
	} 
	return ans;
} 

LCM(最小公倍数)
ll lcm(ll a,ll b)
{
    return a/gcd(a,b)*b;
}

素筛
a[0]=a[1]=1;  真值非素数
for(i=2;i*i<max;i++){
	if(!a[i])
	for(j=i+i;j<max;j+=i){
		a[j]=1;
	}
} 

素数判定 
int prime(int a){
	if(a==0||a==1)
	return 0;
	if(a==2)
	return 1;
	int n=sqrt(a);
	for(int i=2;i<=n;i++){
		if(a%i==0)
		return 0;
		return 1;
	}
}

完全错排 
typedef unsigned long long ull;
ull f[30];     
void list(){
	f[1]=0;f[2]=1;
	for(ull i=3;i<=30;i++){
		f[i]=(i-1)*(f[i-1]+f[i-2]);
	}
}