状压dp之二之三 炮兵阵地/玉米田 By cellur925

一、简单的状压dp 玉米田

题目描述

Farmer John has purchased a lush new rectangular pasture composed of M by N (1 ≤ M ≤ 12; 1 ≤ N ≤ 12) square parcels. He wants to grow some yummy corn for the cows on a number of squares. Regrettably, some of the squares are infertile and can't be planted. Canny FJ knows that the cows dislike eating close to each other, so when choosing which squares to plant, he avoids choosing squares that are adjacent; no two chosen squares share an edge. He has not yet made the final choice as to which squares to plant.

Being a very open-minded man, Farmer John wants to consider all possible options for how to choose the squares for planting. He is so open-minded that he considers choosing no squares as a valid option! Please help Farmer John determine the number of ways he can choose the squares to plant.

农场主John新买了一块长方形的新牧场,这块牧场被划分成M行N列(1 ≤ M ≤ 12; 1 ≤ N ≤ 12),每一格都是一块正方形的土地。John打算在牧场上的某几格里种上美味的草,供他的奶牛们享用。

遗憾的是,有些土地相当贫瘠,不能用来种草。并且,奶牛们喜欢独占一块草地的感觉,于是John不会选择两块相邻的土地,也就是说,没有哪两块草地有公共边。

John想知道,如果不考虑草地的总块数,那么,一共有多少种种植方案可供他选择?(当然,把新牧场完全荒废也是一种方案)

输入输出格式

输入格式:

 

第一行:两个整数M和N,用空格隔开。

第2到第M+1行:每行包含N个用空格隔开的整数,描述了每块土地的状态。第i+1行描述了第i行的土地,所有整数均为0或1,是1的话,表示这块土地足够肥沃,0则表示这块土地不适合种草。

 

输出格式:

 

一个整数,即牧场分配总方案数除以100,000,000的余数。

 

输入输出样例

输入样例#1: 
2 3
1 1 1
0 1 0
输出样例#1: 
9


首先这道题是我自己独立写出来哒没有看题解!撒花!(蒟蒻的欢呼)(虽然过程波折前后历时近两个小时)

这题是一道显然的状压dp,我们可以如法炮制地像互不侵犯把同一行不相邻1的状态先进行预处理。而这道题的特点是给出土地的
贫瘠程度,不能种在贫瘠的地方。

在赋初值(第一行)时,我们也应检查第一行的各个状态是否满足贫瘠限制(开始并没有)。之后便像最简单的普通状压Dp一样,设
f[i][j]为当前行为第i行,当前行处在第j状态的方案数。转移方程有f[i][j]+=f[i-1][k],k枚举上一行状态。

转移前我们也应检查是否满足贫瘠限制(在这里费了大工夫233),最后的方法是写了一个pd函数。有些暴力地取出二进制的每一位
进行比较。然后判断条件时不要像互不侵犯一样判断8个格子!注意审题233!

code

 1 #include<cstdio>
 2 #include<algorithm>
 3 
 4 using namespace std; 
 5 typedef long long ll;
 6 const ll moder=1e8;
 7 
 8 int n,m,cnt,FAKE;
 9 int bin[4000],pos[20];
10 int field[20][20];
11 ll ans,f[15][5000];
12 
13 void pre()
14 {
15     FAKE=(1<<m)-1;
16     for(int i=0;i<=FAKE;i++)
17      if(!(i&(i>>1))) bin[++cnt]=i;
18 }
19 
20 bool pd(int x,int y)
21 {
22     int tmp=bin[y];
23     int pos=m;
24     while(tmp)
25     {
26         int scc=tmp%2;
27         tmp>>=1;
28         if((scc)&&(!field[x][pos])) return 0;
29         pos--;
30     }
31     return 1;
32 }
33 
34 int main()
35 {
36     scanf("%d%d",&n,&m);
37     pre();
38     for(int i=1;i<=n;i++)
39      for(int j=1;j<=m;j++) scanf("%d",&field[i][j]);
40     for(int i=1;i<=cnt;i++)
41      if(pd(1,i)) (f[1][i]=1)%=moder;
42     for(int i=2;i<=n;i++)
43      for(int j=1;j<=cnt;j++)
44      {
45          if(!pd(i,j)) continue;
46         for(int k=1;k<=cnt;k++)
47         {
48              if(bin[j]&bin[k]) continue;
49              (f[i][j]+=f[i-1][k])%=moder;
50         }
51      }
52     for(int i=1;i<=cnt;i++) (ans+=f[n][i])%=moder;
53     printf("%d",ans%moder);
54     return 0;
55 }
View Code

 

二、比较难想的状压 炮兵阵地

题目描述

司令部的将军们打算在N*M的网格地图上部署他们的炮兵部队。一个N*M的地图由N行M列组成,地图的每一格可能是山地(用“H” 表示),也可能是平原(用“P”表示),如下图。在每一格平原地形上最多可以布置一支炮兵部队(山地上不能够部署炮兵部队);一支炮兵部队在地图上的攻击范围如图中黑色区域所示:

如果在地图中的灰色所标识的平原上部署一支炮兵部队,则图中的黑色的网格表示它能够攻击到的区域:沿横向左右各两格,沿纵向上下各两格。图上其它白色网格均攻击不到。从图上可见炮兵的攻击范围不受地形的影响。 现在,将军们规划如何部署炮兵部队,在防止误伤的前提下(保证任何两支炮兵部队之间不能互相攻击,即任何一支炮兵部队都不在其他支炮兵部队的攻击范围内),在整个地图区域内最多能够摆放多少我军的炮兵部队。

输入输出格式

输入格式:

 

第一行包含两个由空格分割开的正整数,分别表示N和M;

接下来的N行,每一行含有连续的M个字符(‘P’或者‘H’),中间没有空格。按顺序表示地图中每一行的数据。N≤100;M≤10。

 

输出格式:

 

仅一行,包含一个整数K,表示最多能摆放的炮兵部队的数量。

 

输入输出样例

输入样例#1: 
5 4
PHPP
PPHH
PPPP
PHPP
PHHP
输出样例#1: 
6

思路与之前并无二至,只是处理位运算有些麻烦。由于这次我们受到的限制增多,受到上下左右两行的限制,所以我们的状态就要
多加一维,记录上一行的情况,在转移的时候还要枚举上两行。

于是我们可以得到转移方程f[i][j][k]=max(f[i][j][k],f[i-1][k][l])
其中i是当前行,j是当前行状态,k是上一行状态,l是上两行状态。

这是整体思路,位运算的细节部分就写在代码里啦qwq。

code
 1 #include<cstdio>
 2 #include<algorithm>
 3 #include<iostream> 
 4 
 5 using namespace std;
 6 
 7 int n,m,cnt,ans;
 8 int bin[2000],nowk[2000],pos[300];
 9 int f[102][1030][1030];
10 char pam[150][20];
11 
12 bool mout(int x,int y)
13 {//判断当前行的当前状态是否为山地不可放 
14     if((x|pos[y])!=pos[y]) return 0;
15     return 1;
16 }
17 
18 bool pd(int x,int y)
19 {//上下两行同列是否能攻击到对方 
20     if(x&y) return 0;
21     return 1;
22 }
23 
24 void getnowk(int x)
25 {
26     int tmp=0;
27     while(x) tmp+=(x&1),x>>=1;
28     nowk[cnt]=tmp;
29 }
30 
31 void pre()
32 {
33     int fAKe=(1<<m)-1;
34     for(int i=0;i<fAKe;i++)
35      if((!(i&(i>>2)))&&(!(i&(i>>1)))) bin[++cnt]=i,getnowk(i);
36      //位运算优先级很低 需要大力加括号 
37 }
38 
39 int main()
40 {
41     scanf("%d%d",&n,&m);
42     pre();
43     for(int i=1;i<=n;i++) scanf("%s",pam[i]+1);//读入地图    
44     for(int i=1;i<=n;i++)
45      for(int j=1;j<=m;j++)
46       if(pam[i][j]=='P') pos[i]=(pos[i]<<1)|1;
47        else pos[i]=pos[i]<<1;
48     //pos[i]记录第i行地图二进制压缩的01串转为的十进制数,1代表平地
49     //0代表山地,这个方法可以记住   
50        
51     for(int i=1;i<=cnt;i++)
52      f[1][i][0]=nowk[i];//预处理第一行 
53      
54     for(int i=2;i<=n;i++)
55      for(int j=1;j<=cnt;j++)//枚举当前行状态 
56      {
57          int state_now=bin[j];
58          if(!mout(state_now,i)) continue;
59          for(int k=1;k<=cnt;k++)//枚举上一行状态 
60          {
61              int state_last=bin[k];
62              if(!mout(state_last,i-1)) continue;
63              if(!pd(state_now,state_last)) continue;
64              if(i==2)//到第3行才能用前两行递推来,这里相当于一个特判 
65              {
66                  f[i][j][k]=f[i-1][k][0]+nowk[j];
67                  continue;
68              }
69             for(int l=1;l<=cnt;l++)
70             {
71                 int state_ll=bin[l];
72                 if(!mout(state_ll,i-2)) continue;
73                 if(!pd(state_now,state_ll)) continue;
74                 if(!pd(state_last,state_ll)) continue;
75                 f[i][j][k]=max(f[i][j][k],f[i-1][k][l]+nowk[j]);
76             }
77         }
78      }
79     
80     for(int i=1;i<=cnt;i++)
81      for(int j=1;j<=cnt;j++)
82       ans=max(ans,f[n][i][j]);
83     printf("%d",ans);
84     return 0;
85 }
View Code

 

 

转载于:https://www.cnblogs.com/nopartyfoucaodong/p/9410366.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/266404.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

服务器php 启动命令_服务端的cli方式运行

既然是结合tp5&#xff0c;咱当然要借鉴tp5的优势# tp5自定义命令行这个需要自己看tp5的官方文档&#xff0c;直接搜索就能找到&#xff0c;自己添加command.php文件&#xff0c;我的如下~~~/*** Time: 13:55*/return [iss\crontab\command\IssServer,iss\crontab\command\IssC…

全国计算机等级考试题库二级C操作题100套(第47套)

第47套&#xff1a; 给定程序中&#xff0c;函数fun的功能是&#xff1a;找出形参s所指字符串中出现频率最高的字母&#xff08;不区分大小写&#xff09;&#xff0c;并统计出其出现的次数。 例如&#xff0c;形参s所指的字符串为&#xff1a;abcAbsmaxless&#xff0c;程序执…

android调用照相机拍照获取照片并做简单剪裁

1.调用系统的照相机程序 Intent intent new Intent(MediaStore.ACTION_IMAGE_CAPTURE); startActivityForResult(intent, TAKE_PHOTO_WITH_DATA); 2.在onActivityResult中处理返回的data final Bitmap photo data.getParcelableExtra("data"); 如果不需要做剪切处…

用户id可以出现在url中吗_下载Google Drive中的文件

通过Python下载Google Drive中的文件的代码如下&#xff1a;download_googledrive.pyimport

LockSupport的源码实现原理以及应用

一、为什么使用LockSupport类 如果只是LockSupport在使用起来比Object的wait/notify简单&#xff0c; 那还真没必要专门讲解下LockSupport。最主要的是灵活性。 上边的例子代码中&#xff0c;主线程调用了Thread.sleep(1000)方法来等待线程A计算完成进入wait状态。如果去掉Thr…

php网址变量怎么输出,【PHP网站】如何使用dedecms v5.7前台模版里输出变量

dedecmsv5.7是PHP网站内容管理系统&#xff0c;本篇文章将介绍如何使用dedecmsv5.7 前台模版来输出变量&#xff0c;具有一定参考意义&#xff0c;感兴趣的朋友了解一下吧。如何在PHP文件查询出来的数据赋值给前端页面展示出来&#xff1f;例如&#xff1a;PHP文件&#xff1a;…

全国计算机等级考试题库二级C操作题100套(第48套)

第48套&#xff1a; 给定程序中&#xff0c;函数fun的功能是&#xff1a;将形参s所指字符串中的数字字符转换成对 应的数值&#xff0c;计算出这些数值的累加和作为函数值返回。 例如&#xff0c;形参s所指的字符串为&#xff1a;abs5def126jkm8&#xff0c;程序执行后的输出结…

Nginx内核参数相关的优化设定

Nginx内核参数在使用的时候有不少问题需要我们解决&#xff0c;其中在优化方面就需要我们格外的注意。在下面就是对Nginx内核参数优化的详细介绍&#xff0c;希望大家有所收获。 关于Nginx内核参数的优化&#xff1a; net.ipv4.tcp_max_tw_buckets 6000 timewait的数量&#x…

python3 爬虫 requests安装_BOSS直聘招聘信息获取之爬虫工具分析

点击蓝色“不太灵光的程序员”关注我哟加个“星标”&#xff0c;每天上午 09:30&#xff0c;干货推送&#xff01;文中使用的组件库仅限于Python语言&#xff0c;由于最近收到一些同学的留言说&#xff0c;按照网上的教程一步一步的学习&#xff0c;却频繁的出现报错&#xff0…

假期周进度总计(四)

本周学习利用OEM工具创建删除表以及对表内数据进行增删改查的操作 一切均通过此工具进行&#xff0c;然后还学到了用SQL *Plus进行相应操作 本周每天4.5个小时&#xff0c;出现的错误就是点击的误操作&#xff0c;可以改正&#xff0c;下周进一步学习Oracle基本操作转载于:http…

Java集合框架(3)

Map(和Collection<E>一样都是集合框架的顶层接口) |--Hashtable:底层是哈希表数据结构&#xff0c;不可以用null对象作为键或值。它是线程同步的。 |--HashMap&#xff1a;底层是哈希表。允许使用null键null值&#xff0c;该集合是不同步的&#xff0c;效率高&#xff0c…

php函数从数组中取出指定的数目,PHP数组函数

1.array_rand()从数组中随机取出一个或多个元素(返回值是&#xff1a;随机元素的键)$arr[js,css,25,php,30];printf(%s,print_r($arr,true));//随机去2个元素$resarray_rand($arr,2);printf(%s,print_r($res,true));2.array_replace()使用后面数组元素相同 key 的值替换 array1…

存储过程 not supported yet_让我们来看看+Redis如何存储和计算一亿用户的活跃度

1前段时间&#xff0c;在网上看到一道面试题&#xff1a;如何用redis存储统计1亿用户一年的登陆情况&#xff0c;并快速检索任意时间窗口内的活跃用户数量。觉得很有意思&#xff0c;就仔细想了下 。并做了一系列实验&#xff0c;自己模拟了下 。还是有点收获的&#xff0c;现整…

HBase查询优化

1.概述 HBase是一个实时的非关系型数据库&#xff0c;用来存储海量数据。但是&#xff0c;在实际使用场景中&#xff0c;在使用HBase API查询HBase中的数据时&#xff0c;有时会发现数据查询会很慢。本篇博客将从客户端优化和服务端优化两个方面来介绍&#xff0c;如何提高查询…

全国计算机等级考试题库二级C操作题100套(第49套)

第49套&#xff1a; 给定程序中&#xff0c;函数fun的功能是&#xff1a;将形参s所指字符串中所有ASCII码值小于97 的字符存入形参t所指字符数组中&#xff0c;形成一个新串&#xff0c;并统计出符合条件的字符个数 作为函数值返回。 例如&#xff0c;形参s所指的字符串为&…

NEC SV8100电话交换机配置梓博电话计费系统

为了节约办公成本&#xff0c;规范电话使用。公司最近为NEC SV8100电话交换机系统上线了一套梓博的电话计费系统。先将配置过程分享给各位。 1、登录电话交换机系统&#xff08;默认用户名tech、密码12345678&#xff09; 2、点击系统数据配置按钮 2、通过10-01选项调整电话交换…

判断闰年 php,PHP怎么判断一年是否为闰年?

判断是否为闰年的条件是满足下列二者条件之一&#xff1a;年号能被4整除、但不能被100整除&#xff0c;二是年号能被4整除&#xff0c;又能被400整除。那么PHP怎么判断一年是否为闰年&#xff1f;下面本篇文章就来给大家介绍一下使用PHP判断一年是否为闰年的方法&#xff0c;希…

redis desktop manager_面试官:Redis分布式锁如何解决锁超时问题?

Java面试笔试面经、Java技术每天学习一点Java面试关注不迷路作者&#xff1a;wangzaiplus来源&#xff1a;https://www.jianshu.com/u/8cb4591440ca一、前言关于redis分布式锁, 查了很多资料, 发现很多只是实现了最基础的功能, 但是, 并没有解决当锁已超时而业务逻辑还未执行完…

全国计算机等级考试题库二级C操作题100套(第50套)

第50套&#xff1a; 给定程序中&#xff0c;函数fun的功能是:有NN矩阵&#xff0c;以主对角线为对称线&#xff0c;对称元素相加并将结果存放在左下三角元素中&#xff0c;右上三角元素置为0。例如&#xff0c;若N3&#xff0c;有下列矩阵&#xff1a; 1 2 3 4 5 6 7 8 9 计算…

《CLR via C#》读书笔记 之 参数

第九章 参数 2013-02-27 9.3 以传引用的方式向方法传递参数 默认情况下&#xff0c;CLR假定所有方法参数都是传值的。当传递引用类型的对象时&#xff0c;也默认是传值的&#xff0c;只不过这个值是引用&#xff08;指针&#xff09;本身。 CLR允许以传引用的方式传递参数。在C…