状压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[],b[];
int main(){
    scanf("%d",&n);
    rep(i,,n-) scanf("%d",&a[i]);
    scanf("%d",&m);
    rep(i,,m-) scanf("%d",&b[i]);
    rep(i,,n-) sum[<<i]=a[i];
    rep(i,n,n+m-) sum[<<i]=-b[i-n];
    mm=<<(n+m);
    rep(i,,mm-){
        t=i&(-i);
        sum[i]=sum[t]+sum[i-t];
        rep(j,,m+n-)
            if ((<<j)&i) f[i]=max(f[i],f[i^(<<j)]);
        if (!sum[i]) ++f[i];
    }
    printf("%d\n",n+m-*f[mm-]);
    return ;
}

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[][maxn],f[][maxn],g[][maxn],d[maxn],map[][],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]==){
        sum+=(<<(x-));
        ans++;
        dfs(x+);
        sum-=(<<(x-));
        ans--;
    }
    dfs(x+);
}
int main(){
    scanf("%d%d",&n,&m);
    rep(i,,n)
        rep(j,,m) scanf("%d",&map[i][j]);
    rep(i,,n){
        tot=ans=sum=;
        y=i;
        dfs();
    }
    ans=;
    if (n==) printf("%d\n",d[]);
    else {
        rep(i,,d[]) f[][i]=;
        rep(i,,n)
            rep(j,,d[i])
                rep(k,,d[i-]){
                    if ((g[i][j]&g[i-][k])==) f[i][j]=(f[i][j]+f[i-][k])%mm;
                }
        rep(i,,d[n]) ans=(ans+f[n][i])%mm;
        printf("%d",ans);
    }
    return ;
}

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][],mm,ans;
int n,m,a[];
int main(){
    scanf("%d%d",&n,&m);
    rep(i,,n) scanf("%d",&a[i]);
    rep(i,,n-) f[<<i][i+]=;
    mm=(<<n)-;
    rep(i,,mm)
        rep(j,,n) if (<<(j-)&i)
            rep(k,,n) if ((<<(k-)|i)!=i&&abs(a[k]-a[j])>m) f[i|<<(k-)][k]+=f[i][j];
    rep(i,,n) ans+=f[mm][i];
    printf("%lld\n",ans);
    return ;
}

2073: [POI2004]PRZ

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

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

（这个比较神奇是这样的，每次&(j-1)可以发现每次去掉了原二进制数中最后一个1，也就是砍掉了这个人。这样的好处在于保证正确性然后。。因为答案是多个集合的答案加起来，如果我从头枚举转移的子集，那时间复杂度是不允许的，但是用上面那种方式转移的话，转移的时间复杂度就只有n*2^n

#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[],b[],t[],w[maxn],tim[maxn],f[maxn];
int main(){
    scanf("%d%d",&m,&n);
    rep(i,,n) scanf("%d%d",&a[i],&b[i]),t[i]=<<(i-);
    mm=(<<n)-;
    rep(i,,mm)
        rep(j,,n) if (t[j]&i) tim[i]=max(tim[i],a[j]),w[i]+=b[j];
    memset(f,/,sizeof(f)); f[]=;
    rep(i,,mm)
        for(int j=i;j;j=i&(j-))
            if (w[j]<=m) f[i]=min(f[i],tim[j]+f[i^j]);
    printf("%d\n",f[mm]);
    return ;
}