hdu 6143

来源:互联网 发布:网络部咨询员提成方案 编辑:程序博客网 时间:2024/06/14 19:27

在m个字母中任取
构成2n长度的名字,前n长度和后n长度的出现的字母不能重复
难点应该在去重
定义数组a[i]为n长度时必定有i个字母的数量

先考虑n长度<=i字母的时候 易得pow(i,n)
对于数组a来说,上式 会出现1个字母,2个字母…n-1个字母的情况
所以a[i]=pow(i,n)-ΣC(i,j)*a[j] (1<=j< n )
而原题的话,就是2n长度,利用组合数求出m个字母取出的不同情况,然后组合起来乘一下就好了

#include<bits/stdc++.h>using namespace std;typedef long long LL;const LL M=1e9+7;#define N 2005LL c[N][N],sum,ans,a[N],midd;LL pow_mod(LL a,LL b){    LL s=1;    while(b){        if(b&1) s=s*a%M;        b>>=1;        a=a*a%M;    }    return s;}int main(){    int t,i,j,n,m;    scanf("%d",&t); for(i=1;i<N;i++){        c[i][0]=c[i][i]=1;        for(j=1;j<i;j++) c[i][j]=(c[i-1][j-1]+c[i-1][j])%M;    }    while(t--){        scanf("%d%d",&n,&m);        a[1]=1; sum=c[m][1];        for(i=2;i<m;i++){            a[i]=pow_mod(i,n)%M;            for (j=1;j<i;j++)            {                a[i]=a[i]-c[i][j]*a[j]%M;                a[i]=a[i]%M;                if (a[i]<0) a[i]=(a[i]+M)%M;            }        }        ans=0;        for (i=1;i<m;i++)            for (j=i;i+j<=m;j++)        {            midd=c[m][i]%M*c[m-i][j]%M*a[i]%M*a[j]%M;            ans=(ans+midd)%M;            if (i!=j) ans=(ans+midd)%M;        }        printf("%lld\n",ans);    }    return 0;}
原创粉丝点击