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