HDU 4832 Chess (动态规划+组合数学)
2015-10-31 16:15
351 查看
关键点是行和列是独立的。
移动K步可以视为x方向移动i步(y向不动),y方向移动K-i步(x向不动)。分别预处理出x和y方向上移动i和K-i步时的种数,然后使用排列组合枚举i。相乘即可。
注意,计算组合数时可能使用了除法这将使模运算出错。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
#include<map>
#include<string>
#include<set>
#include<queue>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn=1010;
const LL mod=9999991;
LL C[maxn][maxn];
LL dp1[maxn][maxn],dp2[maxn][maxn];
LL dp[maxn][maxn];
const int mv[4]= {1,2,-1,-2};
void init()
{
C[0][0]=1;
for(int i=1; i<=1000; ++i)
{
C[i][0]=C[i-1][0];
for(int j=1; j<=i; ++j)
{
C[i][j]=C[i-1][j-1]+C[i-1][j];
C[i][j]%=mod;
}
}
}
int main()
{
init();
int T,kase=0;
scanf("%d",&T);
while(T--)
{
int n,m,K,x,y;
scanf("%d%d%d%d%d",&n,&m,&K,&x,&y);
memset(dp1,0,sizeof(dp1));
memset(dp2,0,sizeof(dp2));
memset(dp,0,sizeof(dp));
dp1[0][0]=dp1[0][x]=1;
for(int i=1; i<=K; ++i)
{
for(int j=1; j<=n; ++j)
{
for(int k=0; k<4; ++k)
{
int s=j+mv[k];
if(1<=s&&s<=n)
dp1[i][j]+=dp1[i-1][s];
dp1[i][j]%=mod;
}
dp1[i][0]+=dp1[i][j];
dp1[i][0]%=mod;
}
}
dp2[0][0]=dp2[0][y]=1;
for(int i=1; i<=K; ++i)
{
for(int j=1; j<=m; ++j)
{
for(int k=0; k<4; ++k)
{
int s=j+mv[k];
if(1<=s&&s<=m)
dp2[i][j]+=dp2[i-1][s];
dp2[i][j]%=mod;
}
dp2[i][0]+=dp2[i][j];
dp2[i][0]%=mod;
}
}
LL ans=0;
for(int i=0; i<=K; ++i)
{
int j=K-i;
ans+=C[K][i]*dp1[i][0]%mod*dp2[j][0]%mod;
ans%=mod;
}
printf("Case #%d:\n%I64d\n",++kase,ans);
}
return 0;
}
移动K步可以视为x方向移动i步(y向不动),y方向移动K-i步(x向不动)。分别预处理出x和y方向上移动i和K-i步时的种数,然后使用排列组合枚举i。相乘即可。
注意,计算组合数时可能使用了除法这将使模运算出错。
#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
#include<map>
#include<string>
#include<set>
#include<queue>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn=1010;
const LL mod=9999991;
LL C[maxn][maxn];
LL dp1[maxn][maxn],dp2[maxn][maxn];
LL dp[maxn][maxn];
const int mv[4]= {1,2,-1,-2};
void init()
{
C[0][0]=1;
for(int i=1; i<=1000; ++i)
{
C[i][0]=C[i-1][0];
for(int j=1; j<=i; ++j)
{
C[i][j]=C[i-1][j-1]+C[i-1][j];
C[i][j]%=mod;
}
}
}
int main()
{
init();
int T,kase=0;
scanf("%d",&T);
while(T--)
{
int n,m,K,x,y;
scanf("%d%d%d%d%d",&n,&m,&K,&x,&y);
memset(dp1,0,sizeof(dp1));
memset(dp2,0,sizeof(dp2));
memset(dp,0,sizeof(dp));
dp1[0][0]=dp1[0][x]=1;
for(int i=1; i<=K; ++i)
{
for(int j=1; j<=n; ++j)
{
for(int k=0; k<4; ++k)
{
int s=j+mv[k];
if(1<=s&&s<=n)
dp1[i][j]+=dp1[i-1][s];
dp1[i][j]%=mod;
}
dp1[i][0]+=dp1[i][j];
dp1[i][0]%=mod;
}
}
dp2[0][0]=dp2[0][y]=1;
for(int i=1; i<=K; ++i)
{
for(int j=1; j<=m; ++j)
{
for(int k=0; k<4; ++k)
{
int s=j+mv[k];
if(1<=s&&s<=m)
dp2[i][j]+=dp2[i-1][s];
dp2[i][j]%=mod;
}
dp2[i][0]+=dp2[i][j];
dp2[i][0]%=mod;
}
}
LL ans=0;
for(int i=0; i<=K; ++i)
{
int j=K-i;
ans+=C[K][i]*dp1[i][0]%mod*dp2[j][0]%mod;
ans%=mod;
}
printf("Case #%d:\n%I64d\n",++kase,ans);
}
return 0;
}
相关文章推荐
- C++动态规划之最长公子序列实例
- C++动态规划之背包问题解决方法
- C#使用动态规划解决0-1背包问题实例分析
- 动态规划
- C++ 动态规划
- DP(动态规划) 解游轮费用问题
- 动态规划的用法——01背包问题
- 动态规划的用法——01背包问题
- 《收集苹果》 动态规划入门
- 《DNA比对》蓝桥杯复赛试题
- 《背包问题》 动态规划
- 初学ACM - 半数集(Half Set)问题 NOJ 1010 / FOJ 1207
- 初学ACM - 组合数学基础题目PKU 1833
- 关于爬楼梯的动态规划算法
- 动态规划 --- hdu 1003 **
- DP问题各种模型的状态转移方程
- 0-1背包解题过程
- 背包问题
- USACO 3.2.2:Stringsobits
- 字符串编辑距离