#include #include #include using namespace std; int F[1<<20]; pairf(long N,int p) { pairret=make_pair(0,1); while(N>=1) { long q=N/p,r=N%p; ret.second=(long)ret.second*atcoder::pow_mod(F[p-1],q,p)%p*F[r]%p; ret.first+=N=q; } return ret; } pairmul(pairret,long v,int p) { while(v%p==0)v/=p,ret.first++; ret.second=(long)ret.second*(v%p)%p; return ret; } int main() { long L,R;int M; cin>>L>>R>>M; vector >fac; { int m=M; for(int i=2;i*i<=m;i++)if(m%i==0) { fac.push_back(make_pair(i,0)); while(m%i==0)m/=i,fac.back().second++; } if(m>1)fac.push_back(make_pair(m,1)); } vectorr,m; for(pairfp:fac) { int p=fp.first,mod=1; for(int i=0;ix=f(2*L,p),y=f(L,p); for(long n=L;n<=R;n++) { if(n>L) { //x <- 2*n-1,2*n //y <- n x=mul(mul(x,2*n-1,p),2*n,p); y=mul(y,n,p); } int tmp=(long)x.second*atcoder::inv_mod((long)y.second*y.second%p,p)%p; assert(x.first-y.first*2>=0); if(x.first-y.first*2>=fp.second)tmp=0; else { for(int i=0;i