【状压dp练习】洛谷P1896 [SCOI2005]互不侵犯 P2704 [NOI2001]炮兵阵地

于 2021-04-07 19:28:47 发布
阅读量139 收藏
点赞数
分类专栏： 【算法】搜索，DP，串
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/wldcmzy/article/details/115493003
版权

可以把一种只有两种状态的一维的数据组用一个整数表示，其二进制的值可以表示每一个位置的状态。如：1010 表示第1、3个位置有，第2、4个位置没有。

洛谷 P1896 [SCOI2005]互不侵犯

init函数枚举所有的m位不存在两个相邻的1的整数，将整数的值记录在status数组中，将此数含有1的数量记录在cnt数组中
dp[i][s][n] 表示进行道第i行时，在第s个状态，使用了n个国王的情况数。

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <queue>
#include <map>
#include <set>
#include <stack>
#include <vector>
typedef std::pair<int ,int> PP;
typedef long long LL;
//const int N = 1e5+5;
LL dp[10][1025][100];
int status[1025],cnt[1025];
int n,k;
int idx;
LL ans;

void init(int v,int b,int total)
{
   if(b >=n){
      status[idx] = v;
      cnt[idx++] = total;
      return ;
   }
   init(v,b+1,total);
   init(v+(1<<b),b+2,total+1);
}

int main()
{
   scanf("%d%d",&n,&k);
   init(0,0,0);
   for(int i=0; i<idx; i++){
      dp[0][i][cnt[i]] += 1;
   }
   for(int r=1; r<n; r++){
      for(int i=0; i<idx; i++){
         for(int j=0; j<idx; j++){
            if(status[i]&status[j] || (status[i]<<1)&status[j] || (status[i]>>1)&status[j]) continue;
            for(int x=cnt[i]; x<=k; x++){
               dp[r][i][x] += dp[r-1][j][x-cnt[i]];
            }
         }
      }
   }
   ans = 0;
   for(int i=0; i<idx; i++) ans += dp[n-1][i][k];
   printf("%lld",ans);
   return 0;
}

P2704 [NOI2001]炮兵阵地

抄完了互不侵犯之后拿这道题练了练手，然后失败。

写一个错误的思路
开二维数组dp[i][s]表示进行到第i行时，当状态为s，可以摆放的最多的炮兵的数量。
初始化前两行后，当一种状态s可以摆放，则dp[i][s] = max(自己，dp[i-2][上两行的状态]+上一行状态的总人数+这一行状态的总人数) 即用上到二行为止的最大人数+上一行某状态的人数+这一行s状态的人数来更新这一行的s状态
这种思路存在一个缺陷，即只记录了单行某个状态的最大人数。
沿用这个思路，当进行到第n行时，我们用“截止到n-2行k状态的最大人数”+“n-1行取状态j的人数”+“第n行状态s的人数”来更新第n行。问题在于：我们无法保证当“截止到n-2行k状态的最大人数”取得时，用状态j当作n-2行的下一行时合法的，因为第n-1行的状态受到n-2行和n-3行限制，有可能当“截止到n-2行k状态的最大人数”取得时，j状态是不能与n-3行的状态共存的。然而在这种思路中，我们下意识地认为这种错误时合法的，所以在多数情况下无法得到正确答案。
更正思路：
开三维数组dp[i][s][last]，表示进行到第i行时，当状态为s，并且上一行状态为last时，可以摆放的最多的炮兵的数量。
错误代码：(竟然还能得20分)

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <queue>
#include <map>
#include <set>
#include <stack>
#include <vector>
//typedef std::pair<int ,int> PP;
typedef long long LL;
const int N = 105,M = 12,ER9 = 1030;
int rows[N],status[ER9],va[ER9];
LL dp[N][ER9];
char ss[12];
int n,m;
int idx;

int id[N],mp[N][ER9];

int toStatus()
{
   int re = 0;
   for(int i=0; i<m; i++){
      re += (ss[i] == 'H') ? (1<<i) : 0;
   }
   return re;
}

void init(int val,int pos,int cnt1)
{
   if(pos >= m){
      status[idx] = val;
      va[idx++] = cnt1;
      return ;
   }
   init(val,pos+1,cnt1);
   init(val+(1<<pos),pos+3,cnt1+1);
}

int main()
{
   scanf("%d%d",&n,&m);
   for(int i=0; i<n; i++){
      scanf("%s",ss);
      rows[i] = toStatus();
   }
   init(0,0,0);
   for(int i=0; i<idx; i++){
      if(status[i]&rows[0]) continue;
      mp[0][id[0]] = i;
      dp[0][id[0]] = va[i];
      id[0] += 1;
   }
   for(int i=0; i<idx; i++){
      if(status[i]&rows[1]) continue;
      mp[1][id[1]] = i;
      for(int j=0; j<id[0]; j++){
         if(status[i]&status[mp[0][j]]) continue;
         dp[1][id[1]] = std::max(dp[1][id[1]],va[i]+dp[0][j]);
      }
      id[1] += 1;
   }
   for(int r=2; r<n; r++){
      for(int i=0; i<idx; i++){
         if(status[i]&rows[r]) continue;
         mp[r][id[r]] = i;
         for(int j=0; j<id[r-1]; j++){
            for(int k=0; k<id[r-2]; k++){
               if(status[i]&status[mp[r-1][j]] || status[i]&status[mp[r-2][k]] || status[mp[r-2][k]]&status[mp[r-1][j]]) continue;
               dp[r][id[r]] = std::max(dp[r][id[r]],va[i]+va[mp[r-1][j]]+dp[r-2][k]);
            }
         }
         id[r] += 1;
      }
   }
   LL ans = 0;
   for(int i=0; i<id[n-1]; i++){
      ans = std::max(ans,dp[n-1][i]);
   }
   printf("%lld",ans);

   return 0;
}

AC代码:
看很多题解都说要用到三行的滚动数组，否则会MLE。
但是这个没有用,直接开了100行，可能是由于动态加状态节省了空间吧。破案了，题目改过，不限制空间了

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <queue>
#include <map>
#include <set>
#include <stack>
#include <vector>
//typedef std::pair<int ,int> PP;
typedef long long LL;
const int N = 105,M = 12,E10 = 1024;
int rows[N],status[E10],va[E10],dp[N][E10][E10];
//rows:把地图的每一行转化为数字，H为1，P为0，储存在rows数组
//status记录所有状态的对应的数字
//va记录所有状态对应的人数
//dp[i][j][k]表示第i行状态为j，且i-1行状态为k时，可以放置的最大炮兵数
char ss[12];
int n,m;
int idx,ans;
//idx记录status数组的大小
int id[N],mp[N][E10];
//id[i]数组记录第i行可行的总状态数
//mp[i][j]记录第i行第j个状态在status数组中对应的编号


int toStatus()
{
   int re = 0;
   for(int i=0; i<m; i++){
      re += (ss[i] == 'H') ? (1<<i) : 0;
   }
   return re;
}

void init(int val,int pos,int cnt1)
{
   if(pos >= m){
      status[idx] = val;
      va[idx++] = cnt1;
      return ;
   }
   init(val,pos+1,cnt1);
   init(val+(1<<pos),pos+3,cnt1+1);
}

int main()
{
   scanf("%d%d",&n,&m);
   for(int i=0; i<n; i++){
      scanf("%s",ss);
      rows[i] = toStatus();
   }
   init(0,0,0);
   for(int i=0; i<idx; i++){
      if(status[i]&rows[0]) continue;
      mp[0][id[0]] = i;
      dp[0][id[0]][0] = va[i];
      id[0] += 1;
   }
   for(int i=0; i<idx; i++){
      if(status[i]&rows[1]) continue;
      mp[1][id[1]] = i;
      for(int j=0; j<id[0]; j++){
         if(status[i]&status[mp[0][j]]) continue;
         dp[1][id[1]][j] = va[i] + dp[0][j][0];
      }
      id[1] += 1;
   }
   for(int r=2; r<n; r++){
      for(int i=0; i<idx; i++){
         if(status[i]&rows[r]) continue;
         mp[r][id[r]] = i;
         for(int j=0; j<id[r-1]; j++){
            for(int k=0; k<id[r-2]; k++){
               if(status[i]&status[mp[r-1][j]] || status[i]&status[mp[r-2][k]] || status[mp[r-2][k]]&status[mp[r-1][j]]) continue;
               dp[r][id[r]][j] = std::max(dp[r][id[r]][j],va[i]+dp[r-1][j][k]);
            }
         }
         id[r] += 1;
      }
   }
   for(int i=0; i<id[n-1]; i++) for(int j=0; j<id[n-2]; j++){
      ans = std::max(dp[n-1][i][j],ans);
   }
   printf("%d",ans);
   return 0;
}

换了一波滚动数组，但是还是有10.+MB

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <queue>
#include <map>
#include <set>
#include <stack>
#include <vector>
//typedef std::pair<int ,int> PP;
typedef long long LL;
const int N = 105,M = 12,E10 = 1024;
int rows[N],status[E10],va[E10],dp[3][E10][E10];
//rows:把地图的每一行转化为数字，H为1，P为0，储存在rows数组
//status记录所有状态的对应的数字
//va记录所有状态对应的人数
//dp[i][j][k]表示第i行状态为j，且i-1行状态为k时，可以放置的最大炮兵数
char ss[12];
int n,m;
int idx,ans;
//idx记录status数组的大小
int id[3],mp[3][E10];
//id[i]数组记录第i行可行的总状态数
//mp[i][j]记录第i行第j个状态在status数组中对应的编号


int toStatus()
{
   int re = 0;
   for(int i=0; i<m; i++){
      re += (ss[i] == 'H') ? (1<<i) : 0;
   }
   return re;
}

void init(int val,int pos,int cnt1)
{
   if(pos >= m){
      status[idx] = val;
      va[idx++] = cnt1;
      return ;
   }
   init(val,pos+1,cnt1);
   init(val+(1<<pos),pos+3,cnt1+1);
}

int main()
{
   scanf("%d%d",&n,&m);
   for(int i=0; i<n; i++){
      scanf("%s",ss);
      rows[i] = toStatus();
   }
   init(0,0,0);
   for(int i=0; i<idx; i++){
      if(status[i]&rows[0]) continue;
      mp[0][id[0]] = i;
      dp[0][id[0]][0] = va[i];
      id[0] += 1;
   }
   for(int i=0; i<idx; i++){
      if(status[i]&rows[1]) continue;
      mp[1][id[1]] = i;
      for(int j=0; j<id[0]; j++){
         if(status[i]&status[mp[0][j]]) continue;
         dp[1][id[1]][j] = va[i] + dp[0][j][0];
      }
      id[1] += 1;
   }
   for(int r=2; r<n; r++){
      int rm3 = r%3,rs1m3 = (r-1)%3, rs2m3 = (r-2)%3;
      memset(dp[rm3],0,sizeof(rm3));
      id[rm3] = 0;
      for(int i=0; i<idx; i++){
         if(status[i]&rows[r]) continue;
         mp[rm3][id[rm3]] = i;
         for(int j=0; j<id[rs1m3]; j++){
            for(int k=0; k<id[rs2m3]; k++){
               if(status[i]&status[mp[rs1m3][j]] || status[i]&status[mp[rs2m3][k]] || status[mp[rs2m3][k]]&status[mp[rs1m3][j]]) continue;
               dp[rm3][id[rm3]][j] = std::max(dp[rm3][id[rm3]][j],va[i]+dp[rs1m3][j][k]);
            }
         }
         id[rm3] += 1;
      }
   }
   int ns1m3 = (n-1)%3, ns2m3 = (n-2)%3;
   for(int i=0; i<id[ns1m3]; i++) for(int j=0; j<id[ns2m3]; j++){
      ans = std::max(dp[ns1m3][i][j],ans);
   }
   printf("%d",ans);
   return 0;
}

2020ICPC 江西省大学生程序设计竞赛L WZB’s Harem

题目链接https://ac.nowcoder.com/acm/contest/15806/L
题目翻译： n皇后问题的退化版——n车问题 n<=20
状压 yyds

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <string>
#include <queue>
#include <vector>
#include <stack>
#include <map>
#include <set>
typedef long long int LL;
const int INF = 0x7fffffff, SCF = 0x3fffffff;
const int N = 21, MOD =  1000000007;
int status[N][1<<N],dp[N][1<<N];
int id[N];
int board[N];
int qu[N];
int n;
int toStatus()
{
    int re = 0;
    for(int i=1; i<=n; i++){
        re = (re<<1)+qu[i];
    }
    return re;
}
int main()
{
    scanf("%d",&n);
    for(int i=1; i<=n; i++){
        for(int j=1; j<=n; j++){
            scanf("%d",&qu[j]);
        }
        board[i] = toStatus();
    }
    for(int i=0; i<(1<<n); i++){
        int tt = i, cnt = 0;
        while(tt > 0){
            cnt += (tt&1);
            tt >>= 1;
        }
        status[cnt][id[cnt]++] = i;
    }
    dp[0][0] = 1; //No.Status 0, id 0;
    for(int i=1; i<=n; i++){
        for(int k=0; k<id[i-1]; k++){
            int pre = status[i-1][k];
            for(int j=0; j<n; j++){
                if((pre>>j)&1) continue;
                if( !( (1<<j) & board[i]) ) {
                    int nxt = pre | (1<<j);
                    dp[i][nxt] = (dp[i][nxt] + dp[i-1][pre]) % MOD;
                }
            }
        }
    }
    LL ans = 1;
    for(int i=2; i<=n; i++){
        ans = ans*i%MOD;
    }
    printf("%lld",ans * dp[n][(1<<n)-1] % MOD);
}
StarsWhisper
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
题解：P2704 [NOI2001] 炮兵阵地 状压dp
轩岚的博客
09-17 157
一篇平平淡淡的题解
洛谷P1896-互不侵犯(状压dp)
Chase的博客
05-23 272
题意：洛谷P1896 在N×NN×NN×N的棋盘里面放KKK个国王，使他们互不攻击，共有多少种摆放方案。国王能攻击到它上下左右，以及左上左下右上右下八个方向上附近的各一个格子，共888个格子。(N&lt;=9,K&lt;=N∗N)(N&lt;=9,K&lt;=N*N)(N<=9,K<=N∗N) 分析： 首先暴搜肯定被否定，时间复杂度是指数级的，所以我们就...
NOI P1896 互不侵犯 状压DP
weixin_34293911的博客
08-09 96
题目描述 在N×N的棋盘里面放K个国王，使他们互不攻击，共有多少种摆放方案。国王能攻击到它上下左右，以及左上左下右上右下八个方向上附近的各一个格子，共8个格子。 注：数据有加强（2018/4/25） 输入输出格式 输入格式： 只有一行，包含两个数N，K （ 1 <=N <=9, 0 <= K <= N * N） 输出格式： 所得的方案数 ...
洛谷 1896 [SCOI 2005] 互不侵犯 (状压DP)
sdxtcqs的博客
08-29 157
https://www.luogu.org/problemnew/show/P1896 题意：求在N×N的棋盘里面放K个国王，使他们互不攻击，共有多少种摆放方案。国王能攻击到它上下左右，以及左上左下右上右下八个方向上附近的各一个格子，共8个格子。其中1&lt;=N&lt;=9,&nbsp;0&lt;=K&lt;=N∗N1&lt;=N&lt;=9,&nbsp;0&lt;=K&lt;=N∗N1 dp...
入门状压DP [SCOI2005] 互不侵犯 [USACO06NOV] 玉米田Corn Fields [NOI2001] 炮兵阵地
蒟蒻的最后的倔强
05-08 190
[SCOI2005]互不侵犯 题目描述 在N×N的棋盘里面放K个国王，使他们互不攻击，共有多少种摆放方案。国王能攻击到它上下左右，以及左上左下右上右下八个方向上附近的各一个格子，共8个格子。 输入输出格式 输入格式： 只有一行，包含两个数N，K （ 1 <=N <=9, 0 <= K <= N * N） 输出格式： 所得的方案数 输入输出样例 输入样例#1： 3 2 输出样...
[SCOI2005]互不侵犯
weixin_30342209的博客
09-10 87
嘟嘟嘟 做完[NOI2001]炮兵阵地这道题，顿时觉得此题好水~~15分钟AC 令dp[i][j][k]表示到第 i 行放了 j 个国王，该行状态为 k 时的方案数。然后dp的时候第一层循环枚举 i，第二层 j，如果j 合法的话，再枚举第三层 i - 1行的状态 k，然后如果k合法且 j 和k不冲突的话，再枚举国王数量h，于是就有： 　　　　dp[i][h][j] += dp...
hehelego#WhyNotMarkdown#SCOI2009-windy数-数位DP经典题1
07-25
[SCOI2009] windy数求十进制数$x$的计数, 满足如下约束无前导0相邻数位上的数码之差至少为2首先, 后两条约束和$a,b$没有关系, 于是, 可
[SCOI2007]压缩 个人理解分析外加部分程序
07-14
这是一个演示文档，内容是本人对SCOI2007的压缩这一题的个人理解、分析以及部分代码，近似于题解。
scoi2006动态最值minmax测试数据
04-25
【标题】"scoi2006动态最值minmax测试数据"涉及的是一个算法竞赛中的问题，其中“scoi2006”很可能指的是2006年斯洛文尼亚信息学奥林匹克竞赛（Slovenian Olympiad in Informatics）。这个题目可能要求参赛者解决一...
soci-backend
04-09
这个小程序运行SCOI Web App的后端API。要求您需要设置一些环境变量才能运行go应用。 如果没有这些设置，它将做出一些（合理的）默认假设： APP_PORT=9000DB_HOST=localhostDB_PORT=3306DB_DATABASE=socidbDB_USER=...
【浏览器插件】智译网页翻译 自动翻译 双语对照 AI对话.zip
08-21
智译网页翻译插件，由新译科技精心打造，是一款智能化的网页翻译利器。当您浏览网页时，它能够自动识别页面语言，并迅速将其翻译成您预设的目标语言。支持双语对照、自动翻译、划词翻译，让您在浏览外文网站时，信息一目了然。 AI助手LangGPT的加入，让智译更进一步。它不仅能自动总结页面内容，还能与您自由对话，高效便捷地获取所需信息。智译支持包括英语、俄语、德语、法语、西班牙语等在内的十数种语言与中文的互译，无论是欧美语言还是亚洲语言，甚至是少数民族语言，都能轻松应对。 无论是查阅资料、浏览新闻还是在线购物，智译都能助您一臂之力，节省时间，提升效率。使用起来也非常简单： 1. 自动翻译的开启与关闭：在设置页面，选择“自动翻译为”选项，即可在打开网页时自动翻译。若不需要自动翻译，关闭此选项即可，也可以点击右下角功能条上的“译”按钮手动翻译。 2. 双语对照翻译的开启与关闭：在设置页面，打开“双语对照模式”，新开或刷新页面后，翻译触发，页面将同时展示原文与译文。点击智译图标可切换查看模式，还可以自定义显示样式。 3. 划词翻译的开启与关闭：在设置页面，启用“启动划词翻译”，即可在页面上划选
毕业设计，基于SpringBoot+Vue+MySQL开发的公寓报修管理系统，源码+数据库+毕业论文+视频演示
08-21
毕业设计，基于SpringBoot+Vue+MySQL开发的公寓报修管理系统，源码+数据库+毕业论文+视频演示 现代经济快节奏发展以及不断完善升级的信息化技术，让传统数据信息的管理升级为软件存储，归纳，集中处理数据信息的管理方式。本公寓报修管理系统就是在这样的大环境下诞生，其可以帮助管理者在短时间内处理完毕庞大的数据信息，使用这种软件工具可以帮助管理人员提高事务处理效率，达到事半功倍的效果。此公寓报修管理系统利用当下成熟完善的Spring Boot框架，使用跨平台的可开发大型商业网站的Java语言，以及最受欢迎的RDBMS应用软件之一的MySQL数据库进行程序开发。公寓报修管理系统有管理员，住户，维修人员。管理员可以管理住户信息和维修人员信息，可以审核维修人员的请假信息，住户可以申请维修，可以对维修结果评价，维修人员负责住户提交的维修信息，也可以请假。公寓报修管理系统的开发根据操作人员需要设计的界面简洁美观，在功能模块布局上跟同类型网站保持一致，程序在实现基本要求功能时，也为数据信息面临的安全问题提供了一些实用的解决方案。可以说该程序在帮助管理者高效率地处理工作事务的同时，也实现了数据信息的整体化，规范化与自动化。 关键词：公寓报修管理系统；Spring Boot框架；MySQL；自动化；VUE
批量flac转mp3脚本
08-21
可以批量将.flac音频文件转换成mp3格式，高效快捷 1、将下载的flac文件放到flac目录下 2、双击flac2mp3.bat 3、等转码完成，到mp3目录下查看即可
【浏览器插件】Ctool程序开发常用工具.zip
08-21
【浏览器插件-Ctool：程序开发常用工具】 ### 功能概览 以下是Ctool插件提供的功能列表，每个功能都支持离线使用，以满足开发者在不同环境下的需求。 | 功能分类 | 支持功能 | 离线使用 | |-----------|-----------|-----------| | **哈希** | `md5`, `sha1`, `sha256`, `sha512`, `sm3`, 支持批量处理及文件哈希计算 | √ | | **加密/解密** | `AES`, `DES`, `RC4`, `Rabbit`, `TripleDes`, `sm2`, `sm4` | √ | | **BASE64编码** | 加密、解密，支持文件操作 | √ | | **URL编码** | 编码、解码 | √ | | **时间** | 时间戳转换、毫秒处理、时区转换、时间计算器 | √ | | **二维码** | 生成、解析 | √ | | **条形码** | 生成 | √ | | **汉字转拼音** | 支持声调、首字母、分隔符 | √ | | **IP地址查询** | 运营商、城市信息查询（不支
java资源数字图书制作工具EpubCheck
08-21
java资源数字图书制作工具 EpubCheck提取方式是百度网盘分享地址
2024汽车AI大模型TOP10分析报告.pdf
08-21
2024汽车AI大模型TOP10分析报告.pdf
huggingface的bert-base-uncased.zip的split的1/2
最新发布
08-21
huggingface的bert-base-uncased.zip的split的1/2
c语言学生成绩管理系统.txt
08-21
为了帮助您实现一个简单的C语言学生成绩管理系统，下面是一个基本的设计方案和示例代码。这个系统能够添加学生信息、查询成绩、修改成绩以及删除学生信息。
洛谷 P1640 BZOJ 1854 [SCOI2010]连续攻击游戏
05-14
这道题目需要进行博弈论的思考。我们可以使用动态规划来解决这道题目。 首先，我们可以定义 $f_{i,j}$ 表示当前剩下 $i$ 到 $j$ 这些数时，当前玩家最多能获得的分数。我们可以发现，在这个状态中，当前玩家只能从 ...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值