标签：ems code name nbsp algo end i++ ucid font
【题目链接】 http://acm.hdu.edu.cn/showproblem.php?pid=5800
【题目大意】
给出一个容量上限s，f[i][j][k][l][m]表示k和l两个物品不能选，i和j两个物品必选，最终质量为m的方案数。求这些方案数的总和。
【题解】
令dp[i][j][s1][s2]表示前i个物品填了j的体积，有s1个物品选为为必选，s2个物品选为必不选的方案数(0<=s1,s2<=2)，则有转移方程dp[i][j][s1][s2]=dp[i-1][j][s1][s2]+dp[i-1][j-a[i]][s1-1][s2]+dp[i-1][j][s1][s2-1]，边界条件为dp[0][0][0][0]=1，时间复杂度O(NS*3^2)。
【代码】
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
#define rep(i,n) for(int i=1;i<=n;i++)
typedef long long LL;
const int N=1005;
const LL mod=1e9+7;
int a[N],n,s,t,T;
LL dp[2][N][3][3],ans;
void add(LL &a,LL b){a=a+b;if(a>mod)a-=mod;}
int main(){
scanf("%d",&T);
while(T--){
scanf("%d%d",&n,&s);
rep(i,n)scanf("%d",a+i);
memset(dp,t=0,sizeof(dp));
dp[0][0][0][0]=1;ans=0;
rep(i,n){
t^=1;
for(int j=s;j>=0;j--){
for(int u=0;u<3;u++)for(int v=0;v<3;v++){
dp[t][j][u][v]=dp[t^1][j][u][v];
if(j>=a[i])add(dp[t][j][u][v],dp[t^1][j-a[i]][u][v]);
if(u&&j>=a[i])add(dp[t][j][u][v],dp[t^1][j-a[i]][u-1][v]);
if(v)add(dp[t][j][u][v],dp[t^1][j][u][v-1]);
}
}
}rep(i,s)add(ans,dp[t][i][2][2]);
ans=((ans*2)%mod)*2%mod;
printf("%lld\n",ans);
}return 0;
}
HDU 5800 To My Girlfriend（单调DP）
标签：ems code name nbsp algo end i++ ucid font
原文地址：http://www.cnblogs.com/forever97/p/hdu5800.html
