近来感觉状态压缩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); } }