(P2774 方格取数问题)<状压DP>

来源:互联网 发布:淘宝网店如何装修店铺 编辑:程序博客网 时间:2024/06/08 14:55

这是一道用状压DP似乎AC无望的题。。。
毕竟题解都是网络流。。。#>#
不过用这个题学会了DP滚动数组的巧妙写法

题目描述

在一个有 m*n 个方格的棋盘中,每个方格中有一个正整数。现要从方格中取数,使任意 2 个数所在方格没有公共边,且取出的数的总和最大。试设计一个满足要求的取数算法。对于给定的方格棋盘,按照取数要求编程找出总和最大的数。

输入输出格式

输入格式:
第 1 行有 2 个正整数 m 和 n,分别表示棋盘的行数和列数。接下来的 m 行,每行有 n 个正整数,表示棋盘方格中的数。

输出格式:
程序运行结束时,将取数的最大总和输出

输入输出样例

输入样例#1:
3 3
1 2 3
3 2 3
2 3 1
输出样例#1:
11

与“互不侵犯”不同的地方:

  1. 很明显的是,求的方向不同,所以DP方程不同,但其他方面的思想都基本一致。本题要把每种方案的结果算出来取最大值。
  2. 初始化的一点小小的不同:每种方案对应的二进制数不仅要用于判断条件(无公共边),还有用于计算每种方案的代数和
while(t){       num[a][st]+=p[a][(t&1)*k];//k代表二进制的第k位,如移到第5位时,k为5       t>>=1;       k++;}
  1. 滚动数组优化:注意到每一行的f都由上一行转移而来,所以可以用滚动数组。但是不能像背包问题那样无脑,因为无论你怎样调整循环的顺序,总是无法保证f数组依然保存上一行的结果。所以需要用另一个数组记录每一行DP完成后的值,下一行DP时调用这个数组里的值,而不是直接调用f。但是这样做需要每一行DP完成后将值赋到另一个数组里,复杂度高。一个比较巧妙的办法是,将f多开一维,新开的这一维只有0和1,代表每一行的准确值和计算值,即原来的f数组和新开的数组。每一次算完后,将0或1与1异或,这样0变成1,1变成0,实现两个数组的相互调用。

只有72分的代码

#include<iostream>#include<cstring>#include<cstdio>#include<algorithm>#include<cmath>using namespace std;int n,m;//正整数 m 和 n,分别表示棋盘的行数和列数int p[21][21];int s[21][2000],num[2000][2000];int st;long long f[2][2000];void _state(){    for(int a=1;a<=m;++a){        st=0;        for(int i=1;i<(1<<n);++i){            if(i&(i<<1)) {//                cout<<a<<" "<<i<<endl;                continue;            }            int t=i;            int k=1;            while(t){                num[a][st]+=p[a][(t&1)*k];                t>>=1;                k++;            }            s[a][st++]=i;        }    }/*    for(int i=0;i<=n;++i){        for(int j=0;j<=st;++j) cout<<num[i][j]<<"("<<s[i][j]<<")"<<" ";        cout<<endl;    }*/}void DP(){    long long ans=0;    int k=1,dm;//    for(int i=1;i<st;++i) f[1][i]+=num[1][i];    for(int i=1;i<=m+5;++i){    ////        dm=k^1;            for(int a=0;a<=st;++a){            for(int b=0;b<=st;++b){                if(s[i][a]&s[i-1][b]) continue;                f[dm][a]=max(f[k][b]+num[i][a],f[dm][a]);                ans=max(ans,f[dm][a]);                //f[k][a]=max(f[dm][a],f[k][a]);//这里不能赋值            }        }        k=dm;    }    printf("%lld",ans);}int main(){    scanf("%d%d",&m,&n);    for(int i=1;i<=m;++i)        for(int j=1;j<=n;++j)            scanf("%d",&p[i][j]);    _state();    DP();    return 0;}
0 0
原创粉丝点击