卢卡斯定理模板
#define ll long long ll Exgcd(ll a,ll b,ll& x,ll& y){ if(a%b==0){ x=0,y=1; return b; } ll r,tx,ty; r=Exgcd(b,a%b,tx,ty); x=ty; y=tx-a/b*ty; } ll Comp(ll a,ll b,ll m){ if(a<b) return 0; if(a==b) return 1; if(b>a-b) b=a-b; ll ans=1,ca=1,cb=1; for(int i=0;i<b;i++){ ca=ca*(a-i)%m; cb=cb*(b-i)%m; } ll x,y; Exgcd(cb,m,x,y); x=(x%m+m)%m; ans=ca*x%m; return ans; } //a^b mod m ll Lucas(ll a,ll b,ll m){ ll ans=1; while(a&&b){ ans=(ans*Comp(a%m,b%m,m))%m; a/=m; b/=m; } return ans; }