近来感觉状态压缩dp的强大性(灵活利用了二进制运算很关键)。。。于是做了俩提来看看。。毕竟队友是专业的dp,我只是管中窥豹下而已。。日后有机会再与之玩耍玩耍。。。ps:如果上天再给我一次机会,当年我愿意选择状态dp而不是网络流(只针对目前比赛出题潮流)
经典问题,不相邻/禁点方案数问题。poj3254
#include<iostream> #include<cstdio> using namespace std; int n,m; int dp[5000][15]; int yu[5000]; int numstate=0; int fib[15]; void init() //n行m列,状态一行推一行 { scanf("%d%d",&n,&m); int maxs=(1<<m); for(int i=0;i<maxs;i++) { if(i&(i<<1)); //预处理掉每行相邻的状态(注意这里的右移用的好) else yu[numstate++]=i; } int tx; for(int i=0;i<n;i++) { for(int j=0;j<m;j++) { scanf("%d",&tx); if(tx==0) //fib【i】记录第i行禁止放的位子 fib[i]=fib[i]|(1<<(m-j-1)); } } } void jhh() { for(int i=0;i<numstate;i++) //首行 { if((yu[i]&fib[0])); else dp[yu[i]][0]=1; } for(int i=1;i<n;i++) { for(int j=0;j<numstate;j++) //到当前行该(合法)状态的种树 =上行合法状态之和 { if(fib[i]&yu[j])continue; for(int k=0;k<numstate;k++) { if(yu[j]&yu[k]||yu[k]&fib[i-1])continue; //合法而且上下不冲突 dp[yu[j]][i]=(dp[yu[j]][i]+dp[yu[k]][i-1])%100000000; } } } int ans=0; for(int i=0;i<numstate;i++) { ans=(ans+dp[yu[i]][n-1])%100000000; } printf("%d\n",ans); } int main() { init(); jhh(); }
hdu 5045, 求哪种排列最优。把复杂度从n!--》2^n*n。
#include<iostream> #include<cstdio> #include<algorithm> #include<cstring> using namespace std; int n,m; double dp[1040][12]; //dp【state】【I】做到第i题的状态的目前最优情况。 int numstate=0; double a[12][1005]; int main() { int T; scanf("%d",&T);int cnt=1; while(T--) { scanf("%d%d",&n,&m); for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) scanf("%lf",&a[i][j]); double sums=0; for(int ii=0;ii<m/n;ii++) //分成n段处理 { for(int yy=0;yy<1029;yy++) //注意点double 型数组 除了0外不能用memset初始化 for(int xx=0;xx<12;xx++) dp[yy][xx]=-1; dp[0][0]=0; int maxst=1<<n; for(int i=0;i<n;i++) for(int j=0;j<maxst;j++) { if(dp[j][i]==-1)continue; for(int k=1,kk=n;k<maxst;k=k<<1,kk--) { if((j&k)==0) //&优先级比比较运算符的低啊!! dp[j|k][i+1]=max(dp[j|k][i+1],dp[j][i]+a[kk][i+1+n*ii]); } } sums=sums+dp[maxst-1][n]; } for(int yy=0;yy<1029;yy++) for(int xx=0;xx<12;xx++) dp[yy][xx]=-1; dp[0][0]=0; int maxst=1<<n; for(int i=0;i<m%n;i++) for(int j=0;j<maxst;j++) { if(dp[j][i]==-1)continue; for(int k=1,kk=n;k<maxst;k=k<<1,kk--) { if((j&k)==0) dp[j|k][i+1]=max(dp[j|k][i+1],dp[j][i]+a[kk][i+1+m/n*n]); } } double maxss=-1; for(int i=0;i<maxst;i++) maxss=max(maxss,dp[i][m%n]); sums=sums+maxss; printf("Case #%d: %.5lf\n",cnt++,sums); } }
原文地址:http://blog.csdn.net/u011498819/article/details/39828305