状压dp初探

写了几道状压。。。然后就一直在颓废。。。

2064: 分裂

http://www.lydsy.com/JudgeOnline/problem.php?id=2064

初始的为正,最后的为负,假设我们能找到k组凑成0的话,答案就是n+m-2*k。于是状压。。其实我一点都不会。。

#include<cstdio>

#include<cstring>

#include<iostream>

#include<algorithm>

#define rep(i,l,r) for (int i=l;i<=r;i++)

#define maxn 1<<22

using namespace std;

int t,mm,f[maxn],sum[maxn],ans,n,m,a[50],b[50];

int main(){

    scanf("%d",&n);

    rep(i,0,n-1) scanf("%d",&a[i]);

    scanf("%d",&m);

    rep(i,0,m-1) scanf("%d",&b[i]);

    rep(i,0,n-1) sum[1<<i]=a[i];

    rep(i,n,n+m-1) sum[1<<i]=-b[i-n];

    mm=1<<(n+m);

    rep(i,1,mm-1){

        t=i&(-i);

        sum[i]=sum[t]+sum[i-t];

        rep(j,0,m+n-1)

            if ((1<<j)&i) f[i]=max(f[i],f[i^(1<<j)]);

        if (!sum[i]) ++f[i];   

    }

    printf("%d\n",n+m-2*f[mm-1]);

    return 0;

}

1725: [Usaco2006 Nov]Corn Fields牧场的安排

http://www.lydsy.com/JudgeOnline/problem.php?id=1725

f[i][j]表示第i行,情况是j的答案。于是裸状压。

#include<cstdio>

#include<cstring>

#include<iostream>

#include<algorithm>

#define rep(i,l,r) for (int i=l;i<=r;i++)

#define maxn 5000

#define mm 100000000

int y,sum,tot,ans,ff[20][maxn],f[20][maxn],g[20][maxn],d[maxn],map[20][20],n,m;

using namespace std;

void dfs(int x){

    if (x>m) {

        d[y]=++tot;

        g[y][tot]=sum;

        ff[y][tot]=ans;

        return;

    }

    if (map[y][x]==1){

        sum+=(1<<(x-1));

        ans++;

        dfs(x+2);

        sum-=(1<<(x-1));

        ans--;

    }

    dfs(x+1);

}

int main(){

    scanf("%d%d",&n,&m);

    rep(i,1,n)

        rep(j,1,m) scanf("%d",&map[i][j]);

    rep(i,1,n){

        tot=ans=sum=0;

        y=i;

        dfs(1);

    }

    ans=0;

    if (n==1) printf("%d\n",d[1]);

    else {

        rep(i,1,d[1]) f[1][i]=1;

        rep(i,2,n)

            rep(j,1,d[i])

                rep(k,1,d[i-1]){

                    if ((g[i][j]&g[i-1][k])==0) f[i][j]=(f[i][j]+f[i-1][k])%mm;

                }  

        rep(i,1,d[n]) ans=(ans+f[n][i])%mm;

        printf("%d",ans);

    }

    return 0;

}

1231: [Usaco2008 Nov]mixup2 混乱的奶牛

http://www.lydsy.com/JudgeOnline/problem.php?id=1231

#include<cstdio>

#include<cstring>

#include<iostream>

#include<algorithm>

#include<cmath>

#define rep(i,l,r) for (int i=l;i<=r;i++)

#define ll long long

#define maxn 70000

ll f[maxn][20],mm,ans;

int n,m,a[20];

int main(){

    scanf("%d%d",&n,&m);

    rep(i,1,n) scanf("%d",&a[i]);

    rep(i,0,n-1) f[1<<i][i+1]=1;

    mm=(1<<n)-1;

    rep(i,0,mm)

        rep(j,1,n) if (1<<(j-1)&i)

            rep(k,1,n) if ((1<<(k-1)|i)!=i&&abs(a[k]-a[j])>m) f[i|1<<(k-1)][k]+=f[i][j];

    rep(i,1,n) ans+=f[mm][i];  

    printf("%lld\n",ans);

    return 0;

}

2073: [POI2004]PRZ

http://www.lydsy.com/JudgeOnline/problem.php?id=2073

主要就是枚举子集转移的时候比较神奇。其实也没什么。。。

#include<cstdio>

#include<cstring>

#include<iostream>

#include<algorithm>

#include<cmath>

#define rep(i,l,r) for (int i=l;i<=r;i++)

#define ll long long

#define maxn 70000

using namespace std;

int mm,n,m,a[20],b[20],t[20],w[maxn],tim[maxn],f[maxn];

int main(){

    scanf("%d%d",&m,&n);

    rep(i,1,n) scanf("%d%d",&a[i],&b[i]),t[i]=1<<(i-1);

    mm=(1<<n)-1;

    rep(i,0,mm)

        rep(j,1,n) if (t[j]&i) tim[i]=max(tim[i],a[j]),w[i]+=b[j];

    memset(f,127/3,sizeof(f)); f[0]=0;

    rep(i,1,mm)

        for(int j=i;j;j=i&(j-1))

            if (w[j]<=m) f[i]=min(f[i],tim[j]+f[i^j]);

    printf("%d\n",f[mm]);

    return 0;

}

时间: 2024-10-10 21:33:51

状压dp初探的相关文章

状压DP初探&#183;总结

2018过农历新年这几天,学了一下状态压缩动态规划,现在先总结一下.   状态压缩其实是一种并没有改变dp本质的优化方法,阶段还是要照分,状态还是老样子,决策依旧要做,转移方程还是得列,最优还是最优,无后还是无后,所以它比较好理解.   状压,顾名思义就是要将一些状压想办法压缩起来(可以压,也可以删).其中这些状态都满足相似性和数量很多.这样才好压而且压得有意义.常见于一般的方格题,网络题等等.   所以一般基础的状压就是将一行的状态压成一个数,这个数的二进制形式反映了这一行的情况.所以位运算可

ZOJ3305Get Sauce 状压DP,

状压DP的题目留个纪念,首先题意一开始读错了,搞了好久,然后弄好了,觉得DFS可以,最后超时,修改了很久还是超时,没办法看了一下n的范围,然后觉得状压可以,但是没有直接推出来,就记忆化搜索了一下,可是一直错,莫名奇妙,然后没办法看了一下题解,发现了下面这个比较好的方法,然后按照这个方程去推,然后敲,也是WA了好多把,写的太搓了,没人家的清楚明了,唉~也算是给自己留个纪念,状压一直做的都不太好~唉~还好理解了, 参考了  http://blog.csdn.net/nash142857/articl

poj 2411 Mondriaan&#39;s Dream(状压DP)

Mondriaan's Dream Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 12232   Accepted: 7142 Description Squares and rectangles fascinated the famous Dutch painter Piet Mondriaan. One night, after producing the drawings in his 'toilet series

(状压dp)uva 10817 Headmaster&#39;s Headache

题目地址 1 #include <bits/stdc++.h> 2 typedef long long ll; 3 using namespace std; 4 const int MAX=1e5+5; 5 const int INF=1e9; 6 int s,m,n; 7 int cost[125]; 8 //char sta[MAX]; 9 string sta; 10 int able[125]; 11 int dp[125][1<<8][1<<8]; 12 in

HDU5816 Hearthstone(状压DP)

题目 Source http://acm.hdu.edu.cn/showproblem.php?pid=5816 Description Hearthstone is an online collectible card game from Blizzard Entertainment. Strategies and luck are the most important factors in this game. When you suffer a desperate situation an

HDU 4336 容斥原理 || 状压DP

状压DP :F(S)=Sum*F(S)+p(x1)*F(S^(1<<x1))+p(x2)*F(S^(1<<x2))...+1; F(S)表示取状态为S的牌的期望次数,Sum表示什么都不取得概率,p(x1)表示的是取x1的概率,最后要加一因为有又多拿了一次.整理一下就可以了. 1 #include <cstdio> 2 const int Maxn=23; 3 double F[1<<Maxn],p[Maxn]; 4 int n; 5 int main() 6

Travel(HDU 4284状压dp)

题意:给n个城市m条路的网图,pp在城市1有一定的钱,想游览这n个城市(包括1),到达一个城市要一定的花费,可以在城市工作赚钱,但前提有工作证(得到有一定的花费),没工作证不能在该城市工作,但可以走,一个城市只能工作一次,问pp是否能游览n个城市回到城市1. 分析:这个题想到杀怪(Survival(ZOJ 2297状压dp) 那个题,也是钱如果小于0就挂了,最后求剩余的最大钱数,先求出最短路和 Hie with the Pie(POJ 3311状压dp) 送披萨那个题相似. #include <

BZOJ 1087: [SCOI2005]互不侵犯King( 状压dp )

简单的状压dp... dp( x , h , s ) 表示当前第 x 行 , 用了 h 个 king , 当前行的状态为 s . 考虑转移 : dp( x , h , s ) = ∑ dp( x - 1 , h - cnt_1( s ) , s' ) ( s and s' 两行不冲突 , cnt_1( s ) 表示 s 状态用了多少个 king ) 我有各种预处理所以 code 的方程和这有点不一样 ------------------------------------------------

BZOJ 1072 排列 状压DP

题意:链接 方法:状压DP? 题解:这题其实没啥好写的,不算很难,推一推就能搞出来. 首先看到这个问题,对于被d整除这个条件,很容易就想到是取余数为0,所以想到可能状态中刚开始含有取余数. 先说我的第一个想法,f[i][j]表示选取i个数且此时的mod为j,这样的思想是第一下蹦出来的,当时想的就是在线来搞最终的答案.不过转瞬即发现,这TM不就是暴力吗魂淡!并没有什么卵用,于是开始想这个状态可不可以做什么优化. 显然第二维的j并不需要太大的优化,暂且先将其搁置一边,来考虑第一维的i怎么优化. 把滚