ACWing 状态压缩DP （完结撒花 day40

伊吹琳琳子

已于 2022-07-27 19:25:42 修改

阅读量227

点赞数

文章标签： c++ 算法蓝桥杯动态规划

于 2022-07-26 02:19:21 首次发布

本文链接：https://blog.csdn.net/qq_23852555/article/details/125976872

版权

确实如y总所说

状态压缩dp真不难

其可以一眼看出的范围

和必须暴力枚举的思路也不是很难想

你问我为啥今天才三题

不知道感觉假期中摆烂还是很正常吧！

明天想出去玩了咕咕咕

我这样真的能拿au吗估计挺难的吧哈哈

AcWing 292. 炮兵阵地

说hard 也不是很 hard

感觉和昨天一样的套路

本来以为四层循环会爆掉

好吧并没有而且我这里用的 g[i-2]都没有越界很神秘

还学到了滚动数组

很厉害呢先照常写写过了再加个& 就可以了呢！

#include <bits/stdc++.h>
using namespace std;
const int N = 110;
const int M = 1<<10;
const int mod = 1e8;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int g[N],f[2][M][M],n,m,cnt[M];
vector<int>state;
bool check(int x){
    for(int i=0;i<m;i++){
        if((x>>i&1)&&((x>>(i+1)&1)||(x>>(i+2)&1)))return false;
    }
    return true;
}
int count(int x){
    int res=0;
    for(int i=0;i<m;i++)if(x>>i&1)res++;
    return res;
}
signed main(){
    fast
    cin>>n>>m;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            char x;cin>>x;
            if(x=='H')g[i]+=1<<(j-1);
        }
    }
    for(int i=0;i<1<<m;i++){
        if(check(i)){
            state.push_back(i);
            cnt[i]=count(i);
        }
    }
    for(int i=1;i<=n;i++){
        for(auto a:state){
            for(auto b:state){
                for(auto c:state){//j i k i-1
                    if(i>=2&&(a&b)==0&&(a&c)==0&&(b&c)==0&&(g[i]&a)==0&&(g[i-1]&b)==0&&(g[i-2]&c)==0)
                        f[i&1][a][b]=max(f[i&1][a][b],f[i-1&1][b][c]+cnt[a]);
                }
            }
        }
    }
    int ans=0;
    for(int i=0;i<1<<m;i++)
        for(int j=0;j<1<<m;j++)
            ans=max(ans,f[n&1][i][j]);
    cout<<ans<<endl;
    return ~~(0^_^0);
}

524. 愤怒的小鸟

卧槽怎么还有卡精度的题啊

说实话这道题 so hard

思路：

我们知道有n个点那么我们可以搞出n2条抛物线（注意有些抛物线可以只过一个点

然后就类似与重复覆盖问题

一般我们都是搜索

不过今天用了状压优化相当于

因为我们这个与顺序无关所以我们可以根据每个状态不同考虑刷表枚举一个没有被打死的猪即可

思路挺简单的可是就是细节有点多

path存的是过了i j 两点的状态卡精度 a要<0 枚举两个点时不能在一条直线上

#include <bits/stdc++.h>
using namespace std;
const int N = 20;
const int M = 1<<18;
const int mod = 1e8;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define x first
#define y second
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int f[M],path[N][N],n,m;//path[]存的是过i j点线的状态
pair<double,double>q[N];
int cmp(double x, double y)
{
    if (fabs(x - y) < 1e-6) return 0;
    if (x < y) return -1;
    return 1;
}
signed main(){
    fast
    int t;cin>>t;
    while(t--){
        cin>>n>>m;
        for(int i=0;i<n;i++)cin>>q[i].x>>q[i].y;
        memset(path,0,sizeof path);
        for(int i=0;i<n;i++){
            path[i][i]=1<<i;
            for(int j=0;j<n;j++){
                double x1=q[i].x,y1=q[i].y,x2=q[j].x,y2=q[j].y;
                double a=((y1/x1)-(y2/x2))/(x1-x2);
                if(!cmp(x1,x2))continue;
                if (cmp(a,0) >= 0) continue;
                double b=(y1/x1)-a*x1;
                int state=0;
                for(int k=0;k<n;k++){
                    double x=q[k].x,y=q[k].y;
                    if(!cmp(a*x*x+b*x,y))state+=1<<k;
                }
                path[i][j]=state;
            }
        }
        memset(f,0x3f3f,sizeof f);
        f[0]=0;
        for(int i=0;i+1<1<<n;i++){
            int x=0;
            for(int j=0;j<n;j++){
                if((i>>j&1)==0){
                    x=j;
                    break;
                }
            }
            for(int j=0;j<n;j++){
                f[i|path[x][j]]=min(f[i|path[x][j]],f[i]+1);
            }
        }
        cout<<f[(1<<n)-1]<<Endl;
    }
    return ~~(0^_^0);
}

P3959 [NOIP2017 提高组] 宝藏

hard

我们都知道暴力怎么打

每个点作为起点爆搜一次

所以状压dp也就自然而然出来了

我们因为有这个条件

新开发一条道路的代价是 \mathrm{L} \times \mathrm{K}L×K。其中 LL 代表这条道路的长度，KK 代表从赞助商帮你打通的宝藏屋到这条道路起点的宝藏屋所经过的宝藏屋的数量（包括赞助商帮你打通的宝藏屋和这条道路起点的宝藏屋）

所以要维护一维dis 自然还要维护一维点最后再维护一维状态

状态转移自然：

f[j][1<<(j-1)|s][depth+1]=sum+dis[i]*g[i][j];

#include <bits/stdc++.h>
using namespace std;
const int N = 15;
const int M = 1<<12;
const int mod = 1e8;
#define int long long
#define endl '\n'
#define Endl '\n'
#define _ 0
#define inf 0x3f3f3f3f3f3f3f3f
#define fast ios::sync_with_stdio(false);cin.tie(nullptr);
int max(int x,int y){return x>y?x:y;}
int min(int x,int y){return x<y?x:y;}
int n,m,ans=inf,g[N][N],dis[N],f[N][M][N];
void dfs(int s,int sum,int depth){
    if(sum>=ans)return;
    if(s==(1<<n)-1){ans=sum;return;}
    for(int i=1;i<=n;i++){
        if(!(1<<(i-1)&s))continue;
        for(int j=1;j<=n;j++){
            if(!(1<<(j-1)&s)&&g[i][j]<2e9){
                if(f[j][1<<(j-1)|s][depth+1]<=sum+dis[i]*g[i][j])continue;
                f[j][1<<(j-1)|s][depth+1]=sum+dis[i]*g[i][j];
                dis[j]=dis[i]+1;
                dfs(1<<(j-1)|s,f[j][1<<(j-1)|s][depth+1],depth+1);
            }
        }
    }
}
signed main(){
    fast
    cin>>n>>m;
    memset(g,0x3f3f,sizeof g);
    for(int i=1;i<=m;i++){
        int x,y,w;cin>>x>>y>>w;
        g[x][y]=g[y][x]=min(w,g[x][y]);
    }
    for(int i=1;i<=n;i++){
        memset(dis,0,sizeof dis);
        memset(f,0x3f3f,sizeof f);
        dis[i]=1;dfs(1<<(i-1),0,0);
    }
    cout<<ans<<endl;
    return ~~(0^_^0);
}

话说这还是一道ac率不高的紫题欸