题目:
题解:
(一):
dp递推:
1.dp[i][j]的状态表示:
第i行一共有j个棋子的数目;
2.转移方程就为:
dp[i][j]=dp[i-1][j-1]*(k-j+1)+dp[i-1][j]。
前一排j-1的个数乘以这一排的方案数(k-j+1)再加上前一层j的个数。
注意:
一般要从右往左遍历比较好。
#include<bits/stdc++.h>
using namespace std;
const int mod = 100003;
int f[2005][2005];
int main()
{
int a,b,c,d,k;
scanf("%d %d %d %d %d",&a,&b,&c,&d,&k);
a+=c;b+=d;
for(int i=0;i<=a;i++)f[i][0]=1;
for(int i=1;i<=c;i++)
{
for(int j=1;j<=i&&j<=k&&j<=d;j++)
{
f[i][j]=(f[i-1][j-1]*(d-j+1)%mod+f[i-1][j])%mod;
//cout<<i<<" --- "<<j<<" :"<<f[i][j]<<endl;
}
//cout<<endl;
}
for(int i=c+1;i<=a;i++)
{
for(int j=1;j<=i&&j<=k&&j<=b;j++)
{
f[i][j]=(f[i-1][j-1]*(b-j+1)%mod+f[i-1][j])%mod;
//cout<<i<<" --- "<<j<<" :"<<f[i][j]<<endl;
}
//cout<<endl;
}
printf("%d\n",f[a][k]);
}
(二):
组合数学:
y总的做法为:
将图形切割为上下两块;
假设上面一块可以放下i个,下面一块可以放下k-i个;枚举每一种i.
对每一种i经行分析:
可以得到:
每一种的answer为C(b,i)*P(a,i) *C(d,k-i)*P(a+c-i,k-i)。
注意要对阶乘做逆元处理。
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int mod = 100003,N=2010;
int f[N],inf[N];
int ksm(int a,int n)
{
int res=1;
while(n)
{
if(n&1)res=(LL)res*a%mod;
a=(LL)a*a%mod;
n=n>>1;
}
//cout<<res<<endl;
return res;
}
int P(int a,int b)
{
if(a>b)return 0;
return (LL)f[b]*inf[b-a]%mod;
}
int C(int a,int b)
{
if(a>b)return 0;
return (LL)f[b]*inf[b-a]%mod*inf[a]%mod;
}
int main()
{
f[0] = inf[0] = 1;
for (int i = 1; i < N; i ++ )
{
f[i] = (LL)f[i - 1] * i % mod;
inf[i] = (LL)inf[i - 1] * ksm(i, mod - 2) % mod;
//cout<<f[i]<<"--"<<inf[i]<<endl;
}
//cout<<C(2,4)<<" "<<P(2,4)<<endl;
int a,b,c,d,k;
cin>>a>>b>>c>>d>>k;
int res=0;
for(int i=0;i<=k;i++)
{
res=(res+(LL)C(i,b)*P(i,a)%mod*C(k-i,d)%mod*P(k-i,a+c-i)%mod)%mod;
}
cout<<res<<endl;
}