题意:求在一个布满方格的矩阵中放置黑点的方案数,时矩阵中每一个n*n的子矩阵中的黑点数都为k。
做法:推导一下发现每一列A[i]放置的黑点数等于A[i+n],这样就可以DP了,有大神的题解,我就不多说了
#include <cstdio>
#include <iostream>
#define mod 1000000007
#define LL long long
const int LMT=102;
using namespace std;
LL c[LMT][LMT],dp[LMT][10004],val[LMT][2];
void init(void)
{
for(int i=0;i<LMT;i++)c[i][0]=c[i][i]=1;
for(int i=1;i<LMT;i++)
for(int j=1;j<i;j++)
{
c[i][j]=c[i-1][j-1]+c[i-1][j];
c[i][j]%=mod;
}
}
LL mypow(LL x,LL y)
{
LL ret=1;
while(y)
{
if(y&1)ret*=x;
ret%=mod;
x*=x;
x%=mod;
y>>=1;
}
return ret;
}
int main(void)
{
int n,k;
init();
LL p,m;
scanf("%d%I64d%d",&n,&m,&k);
dp[0][0]=1;
for(int i=0;i<=n;i++)
{
val[i][0]=mypow(c[n][i],m/n);
val[i][1]=(val[i][0]*c[n][i])%mod;
}
for(int i=1;i<=n;i++)
{
for(int j=0;j<=k;j++)
for(int t=j;t>=0&&j-t<=n;t--)
{
p=(m%n>=i);
dp[i][j]+=dp[i-1][t]*val[j-t][p];
dp[i][j]%=mod;
}
}
cout<<dp[n][k]<<endl;
return 0;
}