#include #include #include #include using namespace std; using mint=atcoder::modint; mint F[1<<20]; pairf(long N,int p) { pairret=make_pair(0L,mint(1)); while(N>=1) { long q=N/p,r=N%p; ret.second*=F[p-1].pow(q)*F[r]; ret.first+=N=q; } return ret; } pairmul(pairret,long v,int p) { while(v%p==0)v/=p,ret.first++; ret.second*=v; 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; mint::set_mod(p); 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=(x.second/y.second.pow(2)).val(); assert(x.first-y.first*2>=0); if(x.first-y.first*2>=fp.second)tmp=0; else { for(int i=0;iret=atcoder::crt(r,m); assert(ret.second==M); cout<<(ret.first-(R-L+1)*2%M+M)%M<