这道题涉及的知识点有点多。。。

所以还是比较有意思的。

domino

描述

迈克生日那天收到一张 N*N 的表格(1 ≤ N ≤ 2000),每个格子里有一个非

负整数(整数范围 0~1000),迈克不喜欢数值太大,他手上有 K 块日字形多米

诺骨牌(1 ≤ K ≤ 8),可以完美覆盖两个相邻格子(上下左右均可)**

问迈克把 K 块骨牌全部放在表格上,使得表格可见整数和最小。

输入

第一行两个数 N,K 用空格隔开 下面 N 行 N 列为该初始表格。

输出

被 K 块骨牌挡住之后剩余数字之和。

分数分布

对于 70%数据,K≤5。

样例输入 1

3 1

2 7 6

9 5 1

4 3 8

样例输出 1

31

样例输入 2

4 2

1 2 4 0

4 0 5 4

0 3 5 1

1 0 4 1

样例输出 2

17

解释下题意,在棋盘上放上k个1*2的骨牌,使覆盖到的格子上的值得和最大。

乍一看,貌似可以贪心,枚举所有骨牌,选择前k大的骨牌。(接下来骨牌大小的定义都是它覆盖的格子的和的大小)。但第二个样例却过不到。

原因很简单,当我们选择了一个骨牌时,会导致顶多7个骨牌不能被选择,如图所示:



(这个骨牌自身也算一个)

如果我们选取的前k大个骨牌在这7个骨牌之中,则不合法,所以导致了错误。

所以我们可以将这7个骨牌全都加入在我们考虑选择的骨牌之中。于是我们考虑的骨牌最多便有了\(50\)个\((7*(8-1)+1)\)

考虑在这\(50\)个骨牌中间深搜,选择其中的\(k\)个,并且保证这些骨牌不会相互覆盖即可。

于是我们就将问题转换成了从50个物品中选出K个,使和最大,同时选择一个物品会导致多个物品不能选择的问题。这明显可以使用状压DP暴力解决。

但是直接暴力会超时,时间复杂度为\(2^{50}\),空间复杂度也实在难以接受。

可以先看看我暴力的代码。

//感觉很好理解所以没怎么打注释
#include <iostream>
#include <cstdio>
#include <vector>
#include <algorithm>
using namespace std; #define N 2010 int A[N][N],n,k,P,siz;
bool mark[51][51]; struct node {
LL x,y,sum,dir;//使用node来存储每个骨牌
node() {}; //x,y表示这个骨牌左上的坐标
node(LL X,LL Y,LL D) {
x=X;y=Y;
sum=A[x][y]+A[x+D][y+1-D];//使用dir变量来表示可能的
dir=D; //两种骨牌
}
}; vector<node> L; vector< pair<LL,LL> > getFuck(node F) {//返回一个vector生成
vector<pair<LL,LL> > D; //一个骨牌的两个坐标
D.push_back( make_pair(F.x,F.y));
D.push_back( make_pair(F.x+F.dir,F.y+1-F.dir));
return D;
} bool Fuck(node A,node B) {
vector< pair<LL,LL> > X=getFuck(A);
vector< pair<LL,LL> > Y=getFuck(B);
for(LL i=0;i<X.size();i++)//判断两个骨牌的4个覆盖区域的
for(LL j=0;j<Y.size();j++)//坐标是否重复
if(X[i]==Y[j])
return 1;
return 0;
} bool cmp(node A,node B) {
return A.sum>B.sum;
} int ans,sum;
bool vis[60]; void dfs(int pos,int haveC) {//暴力深搜,枚举选和不选当前骨牌的情况
ans=max(ans,sum);
if(pos==siz) return ;
dfs(pos+1,haveC);
if(haveC==k) return ;
for(int i=0;i<pos;i++) if(vis[i] && mark[pos][i]) return ;
vis[pos]=1;
sum+=L[pos].sum;
vis[pos]=1;
dfs(pos+1,haveC+1);
vis[pos]=0;
sum-=L[pos].sum;
} int main() {//基操
cin>>n>>k;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
cin>>A[i][j],P+=A[i][j];
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++) {
if(i<n) L.push_back(node(i,j,1));
if(j<n) L.push_back(node(i,j,0));
}
sort(L.begin(),L.end(),cmp);
//for(int i=0;i<L.size();i++) cout<<L[i].sum<<' ';
//cout<<endl;
int Sz=L.size();
siz=min(Sz,50);
for(int i=0;i<siz;i++)
for(int j=0;j<siz;j++)
if(i!=j && Fuck(L[i],L[j]))
mark[i][j]=1;
dfs(0,0);
cout<<P-ans;
return 0;

正解用的是双向DFS,但我剪了下枝,发现用位运算也可以卡过去。

#include <iostream>
#include <cstdio>
#include <vector>
#include <algorithm>
using namespace std; #define N 2010
#define LL long long
LL A[N][N],n,k,P,siz,lol;
LL mark[51]; struct node {
LL x,y,sum,dir;
node() {};
node(LL X,LL Y,LL D) {
x=X;y=Y;
sum=A[x][y]+A[x+D][y+1-D];
dir=D;
}
}; vector<node> L; vector< pair<LL,LL> > getFuck(node F) {
vector<pair<LL,LL> > D;
D.push_back( make_pair(F.x,F.y));
D.push_back( make_pair(F.x+F.dir,F.y+1-F.dir));
return D;
} bool Fuck(node A,node B) {
vector< pair<LL,LL> > X=getFuck(A);
vector< pair<LL,LL> > Y=getFuck(B);
for(LL i=0;i<X.size();i++)
for(LL j=0;j<Y.size();j++)
if(X[i]==Y[j])
return 1;
return 0;
} bool cmp(node A,node B) {
return A.sum>B.sum;
} LL ans,sum; void dfs(LL pos,LL haveC,LL cover) {
ans=max(ans,sum);
if(pos==siz) return ;
if(haveC==k) return ;
dfs(pos+1,haveC,cover);
if(!(cover & (1ll<<pos))) {
sum+=L[pos].sum;
dfs(pos+1,haveC+1,cover|mark[pos]);
sum-=L[pos].sum;
}
} int main() {
cin>>n>>k;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
cin>>A[i][j],P+=A[i][j];
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++) {
if(i<n) L.push_back(node(i,j,1));
if(j<n) L.push_back(node(i,j,0));
}
sort(L.begin(),L.end(),cmp);
lol=L[0].sum;
//for(int i=0;i<L.size();i++) cout<<L[i].sum<<' ';
//cout<<endl;
LL Sz=L.size();
siz=min(Sz,(LL)50);
for(int i=0;i<siz;i++)
for(int j=0;j<siz;j++)
if(i!=j && Fuck(L[i],L[j]))
mark[i]|=(1ll<<j);
dfs(0,0,0);
cout<<P-ans;
return 0;
}

说下我在做这题时的细节吧

1.所有变量尽量都要开long long,因为答案最终有可能爆int

2.第一部分不一定要占一半,开1<<20不会爆空间。

3.位运算要使用\(1ll<<50\).

这里附上双向DFS的代码,有兴趣的同学可以康康,注释写得挺详细的。

(估计也没人看吧hhh)

#include <iostream>
#include <cstdio>
#include <vector>
#include <algorithm>
using namespace std; #define N 2010
#define LL long long LL A[N][N],n,k,P,siz,asiz,bsiz;
LL mark[54];
LL FFuck[1<<21][9],ans; struct node {
LL x,y,sum,dir;
node() {};
node(LL X,LL Y,LL D) {
x=X;y=Y;
sum=A[x][y]+A[x+D][y+1-D];
dir=D;
}
}; vector<node> L; vector< pair<LL,LL> > getFuck(node F) {
vector<pair<LL,LL> > D;
D.push_back( make_pair(F.x,F.y));
D.push_back( make_pair(F.x+F.dir,F.y+1-F.dir));
return D;
} bool Fuck(node A,node B) {
vector< pair<LL,LL> > X=getFuck(A);
vector< pair<LL,LL> > Y=getFuck(B);
for(LL i=0;i<X.size();i++)
for(LL j=0;j<Y.size();j++)
if(X[i]==Y[j])
return 1;
return 0;
} bool cmp(node A,node B) {
return A.sum>B.sum;
} //дµÃÕæ³ó¹þ¹þ¹þ
void dfs1(LL pos,LL haveC,LL state,LL cover,LL sum) {
if(haveC>k) return ;
if(pos==asiz) {
FFuck[state][haveC]=max(FFuck[state][haveC],sum);
return ;
}
dfs1(pos+1,haveC,state,cover,sum);
if(!(cover & (1ll<<pos))) {//
LL mss=L[pos].sum;
dfs1(pos+1,haveC+1,state|(1ll<<pos),cover|mark[pos],sum+mss);
}
} void dfs2(LL pos,LL haveC,LL state,unsigned LL cover,LL sum) {
if(haveC>k) return ;
if(pos==siz) {
ans=max(ans,sum+FFuck[(~cover)&((1ll<<asiz)-1)][k-haveC]);
return ;
}
dfs2(pos+1,haveC,state,cover,sum);
if(!(cover & (1ll<<pos))) {
LL mss=L[pos].sum;
dfs2(pos+1,haveC+1,state|(1ll<<pos),cover|mark[pos],sum+mss);
}
} int main() {
cin>>n>>k;
if(n==1) {cout<<0; return 0;}
LL MM=(k-1)*7+1;
for(LL i=1;i<=n;i++)
for(LL j=1;j<=n;j++)
cin>>A[i][j],P+=A[i][j];
for(LL i=1;i<=n;i++)
for(LL j=1;j<=n;j++) {
if(i<n) L.push_back(node(i,j,1));
if(j<n) L.push_back(node(i,j,0));
}
sort(L.begin(),L.end(),cmp);
//for(LL i=0;i<L.size();i++) cout<<L[i].sum<<' ';
//cout<<endl;
LL Sz=L.size();
siz=min(Sz,MM);
asiz=min(siz/2,(LL)20);
for(LL i=0;i<siz;i++)
for(LL j=0;j<siz;j++)
if(i!=j && Fuck(L[i],L[j]))
mark[i]|=(1ll<<j);
dfs1(0,0,0,0,0);
for(LL c=0;c<=k;c++) {
for(LL i=1;i<(1ll<<asiz);i++) {
LL j=i;
while(j) {
LL t=j&-j;
FFuck[i][c]=max(FFuck[i-t][c],FFuck[i][c]);
j-=t;
}
}
}
dfs2(asiz,0,0,0,0);
cout<<P-ans;
return 0;
}

模拟赛T5 : domino ——深搜+剪枝+位运算优化的更多相关文章

  1. POJ - 3074 Sudoku (搜索)剪枝+位运算优化

    In the game of Sudoku, you are given a large 9 × 9 grid divided into smaller 3 × 3 subgrids. For exa ...

  2. Hdu3812-Sea Sky(深搜+剪枝)

    Sea and Sky are the most favorite things of iSea, even when he was a small child.  Suzi once wrote: ...

  3. php实现不用加减乘除号做加法(1、善于寻找资源:去搜为什么位运算可以实现加法,里面讲的肯定要详细一万倍)

    php实现不用加减乘除号做加法(1.善于寻找资源:去搜为什么位运算可以实现加法,里面讲的肯定要详细一万倍) 一.总结 1.善于寻找资源:去搜为什么位运算可以实现加法,里面讲的肯定要详细一万倍 二.ph ...

  4. 数独求解问题(DFS+位运算优化)

    In the game of Sudoku, you are given a large 9 × 9 grid divided into smaller 3 × 3 subgrids. For exa ...

  5. N皇后-位运算优化

    N皇后问题 时间限制: 5 Sec  内存限制: 128 MB 题目描述 魔法世界历史上曾经出现过一个伟大的罗马共和时期,出于权力平衡的目的,当时的政治理论家波利比奥斯指出:“事涉每个人的权利,绝不应 ...

  6. N皇后解法以及位运算优化

    N皇后解法以及位运算优化 观察棋盘,要求皇后之间不能处在同行同列同一条斜线,求使得每行都有一个皇后的放置方法共有多少种. 每尝试放置一个皇后,都可以把该位置所在的行.列标号用一个数组标记,含义表示该行 ...

  7. POJ-1724 深搜剪枝

    这道题目如果数据很小的话.我们通过这个dfs就可以完成深搜: void dfs(int s) { if (s==N) { minLen=min(minLen,totalLen); return ; } ...

  8. poj1190 生日蛋糕(深搜+剪枝)

    题目链接:poj1190 生日蛋糕 解题思路: 深搜,枚举:每一层可能的高度和半径 确定搜索范围:底层蛋糕的最大可能半径和最大可能高度 搜索顺序:从底层往上搭蛋糕,在同一层尝试时,半径和高度都是从大到 ...

  9. UVA 10160 Servicing Stations(深搜 + 剪枝)

    Problem D: Servicing stations A company offers personal computers for sale in N towns (3 <= N < ...

随机推荐

  1. 2.golang应用目录结构和GOPATH概念

    golang 语言有一个GOPATH的概念就是当前工作目录 [root@localhost golang_test]# tree . ├── bin │   └── hello ├── first.g ...

  2. vue.js(10)--案例--列表增加与删除

    品牌管理案例 (1)bootstrip快速布局 <div class="app"> <div class="panel panel-primary&qu ...

  3. jQuery的加法运算,val()获取的结果相加变成了字符串连接。

    加法运算 ?想必大家听到这都会不屑了,加法运算这是多么简单的运算.然而有的时候在jQuery也让人挺头疼的. 常规的是: var num1 = 123; var num2=123; var total ...

  4. python 中PIL.Image和OpenCV图像格式相互转换

    PIL.Image转换成OpenCV格式: import cv2 from PIL import Image import numpy   image = Image.open("plane ...

  5. 22pygame 安装

    实战步骤 pygame 快速体验 飞机大战 实战 确认模块 --pygame pygame 就是一个 Python 模块, 专为电子游戏设计 提示 : 学习第三方模块, 通常最好的参考资料就在官方网站 ...

  6. 1145. Hashing - Average Search Time (25)

    The task of this problem is simple: insert a sequence of distinct positive integers into a hash tabl ...

  7. alert(1) to win 7

    function escape(s) { // Pass inn "callback#userdata" var thing = s.split(/#/); if (!/^[a-z ...

  8. CH5E26 扑克牌 (计数类DP)

    $ CH~5E26~\times ~ $ 扑克牌: (计数类DP) $ solution: $ 唉,计数类DP总是这么有套路,就是想不到. 这道题我们首先可以发现牌的花色没有价值,只需要知道每种牌有 ...

  9. 碎片记录——JMeter之 http post json对象与参数化调用,以及beanshell 引用Java源码

    参考文档 http://jmeter.apache.org/usermanual/component_reference.html#samplers https://blog.csdn.net/qq_ ...

  10. 【leetcode】778. Swim in Rising Water

    题目如下: 解题思路:本题题干中提到了一个非常重要的前提:"You can swim infinite distance in zero time",同时也给了一个干扰条件,那就是 ...