ACM-ICPC 2017 Asia Urumqi,这道题目用到了逆元还有几何计算,需要推出公式
逆元模板
1.(A/B)%mod==(A%mod) * finv(B,mod) %mod
2.不互质:a/b %mod= a%(b*c) /b
费马小定理:(要求是质数,慎用)
#define ll long long
using namespace std;
ll extgcd(ll a,ll b,ll &x,ll &y)//返回a和b的最大公约数
{
if(b==0)
{
x=1;y=0;
return a;
}
ll q=extgcd(b,a%b,x,y);
ll t=x;x=y;y=t-(a/b)*y;
return q;
}
ll finv(ll a,ll m)//最小逆元
{
ll x,y;
ll g=extgcd(a,m,x,y);
x=(x%m+m)%m;
return x;
}
2.扩展欧几里得:b与mod互质
ll finv(ll a,ll b)
{
ll res=1;
b=b-2;
while(b)
{
if(b&1)
res=res*a%mod;
a=a*a%mod;
b>>=1;
}
return res;
}