额,首先要推个公式。。就是一个T1循环后面跟着m个TYPE2循环的复杂度(显示语句次数)
这公式刚开始也不会推了,后来随便整,发现了个规律
ANS = C(m,0)+C(m,1)+...+C(m+n-1,n-1)
组合数学当中有个公式,C(M,0)+C(M,1)+....+C(m+n-1,n-1) = C(m+n-1,n)
接下来就是求C(M+N-1,N)%MOD
出题人太坑爹了,模数竟然能够分解成97和3761599。。。。实在缺德
这样的话就是求组合数模两个素数的乘积,我们令这个结果为x,组合数为C
设 C% p1 =m1 ,C%p2 =m2
则 x = m1 (mod p1),x= m2(mod p2)
用中国剩余定理求x即可,求模97的余数用LUCAS定理,求大素数的余数直接搞逆元
LUCAS定理一直搞错了,是这样 L(N,M) = C(N%P,M%P) * L(N/P,M/P) ......附代码:
#include <iostream>
#include <cstdio>
#define ll long long
#define ss(a) scanf("%d",&a)
#define N 1100005
#define p1 97
#define p2 3761599
#define mod (p1*p2)
using namespace std;
int a[20],c[100][100];
ll y1,y2,mul[N];
void extend(ll a0,ll b0,ll &x,ll &y)
{
if (!b0)
{
x=1;y=0;
}
else
{
extend(b0,a0%b0,x,y);
ll temp=x;x=y;
y=temp-a0/b0*y;
}
}
ll niyuan(ll t,ll p)
{
ll x,y;
extend(t,p,x,y);
return (x+p)%p;
}
ll lucas(ll x,ll y)
{
if(y>x) return 0;
ll re=1,dx,dy;
dx = x%p1;
dy = y%p1;
re*=c[dx][dy];
dx = x/p1;
dy = y/p1;
if(dx || dy)
{
re*=lucas(dx,dy);
}
re%=p1;
return re;
}
ll C(int n,int m)
{
ll r1,r2,re;
r1=lucas(n,m);
r2=mul[n]*niyuan(mul[m]*mul[n-m]%p2,p2)%p2;
re=(r1*y1%mod+r2*y2%mod)%mod;
return re;
}
void init()
{
y1=p2*niyuan(p2,p1)%mod;
y2=p1*niyuan(p1,p2)%mod;
int i,j;
for (i=1;i<100;i++)
for (j=0;j<=i;j++)
{
if (j==0) c[i][j]=1;
else if (j>i-j) c[i][j]=c[i][i-j];
else c[i][j]=(c[i-1][j-1]+c[i-1][j])%p1;
}
mul[0]=1;
for (i=1;i<N;i++) mul[i]=(mul[i-1]*i)%p2;
}
int main()
{
int i,j,T,cas,n,k,m;
ss(T);
init();
for (cas=1;cas<=T;cas++)
{
ss(n);ss(m);ss(k);
for (i=1;i<=k;i++) ss(a[i]);
a[k+1]=m;
ll res=1;
for (i=1;i<=k;i++)
{
int t=a[i+1]-a[i];
if (t==1) res=(res*n)%mod;
else res=res*C(n+t-1,t)%mod;
}
printf("Case #%d: %d\n",cas,res);
}
return 0;
}