码迷,mamicode.com
首页 > 其他好文 > 详细

状态压缩dp poj 3254 hdu5045

时间:2014-10-07 11:59:13      阅读:157      评论:0      收藏:0      [点我收藏+]

标签:io   os   for   sp   问题   c   on   r   amp   

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


状态压缩dp poj 3254 hdu5045

标签:io   os   for   sp   问题   c   on   r   amp   

原文地址:http://blog.csdn.net/u011498819/article/details/39828305

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!