[HNOI 2007]神奇游乐园

#include <bits/stdc++.h>
#define maxn 110
using namespace std;
typedef long long ll;
int n, m;
int a[maxn][maxn];
ll ans = -1ll << 60;
#define M 2000010
#define mod 997 struct Hashmap{ ll st[M];
int h[1000], size, nxt[M];
ll f[M];
void clear(){memset(h, 0, sizeof h); size = 0;}
void push(ll hash_, ll val){
int tmp = hash_ % mod;
for(int i = h[tmp]; i; i = nxt[i]){
if(st[i] == hash_){
f[i] = max(f[i], val);
return;
}
}
int now = ++ size;
f[now] = val;
st[now] = hash_;
nxt[now] = h[tmp];
h[tmp] = now;
}
}dp[2]; int cur, code[20], ch[20]; void Decode(ll st){
for(int i = m; i >= 0; i --)
code[i] = st & 7, st >>= 3;
} ll Encode(){
ll ret = 0;
memset(ch, -1, sizeof ch);
ch[0] = 0; int cnt = 0;
for(int i = 0; i <= m; i ++){
if(ch[code[i]] == -1)ch[code[i]] = ++ cnt;
code[i] = ch[code[i]];
ret = ret << 3 | code[i];
}
return ret;
} void Shift(){
for(int i = m; i >= 1; i --)
code[i] = code[i-1];
code[0] = 0;
} inline void Change(int u, int v){
for(int i = 0; i <= m; i ++)
if(code[i] == v)
code[i] = u;
} void DP(int i, int j){
dp[cur^1].clear();
for(int k = 1; k <= dp[cur].size; k ++){
Decode(dp[cur].st[k]);
if(j == 1){if(code[m])continue;Shift();}
int Left = code[j-1], Up = code[j];
if(Left && Up){
code[j] = code[j-1] = 0;
if(Left == Up){
ll ENCODE = Encode();
if(ENCODE == 0)ans = max(ans, dp[cur].f[k] + a[i][j]);
}
else{
Change(Left, Up);
dp[cur^1].push(Encode(), dp[cur].f[k] + a[i][j]);
}
}
else if(Left || Up){
int tmp = Left ? Left : Up;
code[j-1] = 0, code[j] = tmp;
dp[cur^1].push(Encode(), dp[cur].f[k] + a[i][j]);
code[j] = 0, code[j-1] = tmp;
dp[cur^1].push(Encode(), dp[cur].f[k] + a[i][j]);
}
else{
dp[cur^1].push(Encode(), dp[cur].f[k]);
code[j] = code[j-1] = 8;
dp[cur^1].push(Encode(), dp[cur].f[k] + a[i][j]);
}
}
cur ^= 1;
} int main(){
#ifndef ONLINE_JUDGE
freopen("park.in","r",stdin);
freopen("park.out","w",stdout);
#endif
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; i ++)
for(int j = 1; j <= m; j ++)
scanf("%d", &a[i][j]);
dp[cur].clear();
dp[cur].push(0, 0);
for(int i = 1; i <= n; i ++)
for(int j = 1; j <= m; j ++)
DP(i, j); printf("%lld\n", ans);
return 0;
}

粘了一个模板上来=-=

基于连通性的动态规划,最小表示法很好用。

我们可以用一个压缩的数字表示一个连通情况,比如DP生成树,概率,棋盘上格子的情况,等等等

只要和连通性有关而且n很小时就可以用啦QAQ。

又忘了模板了QAQ

就是如果左边和上面是一个连通分量即Left == Up时,我们要合并连通分量,所以此时已经出现了一个圈了,这道题不能有多个圈,所以最后不放进去

code[j] = code[j-1] = 0.

[BZOJ 3753]Wall

扩展一下方格变成一个回路问题,然后射线法判断格子是否内部。

出现的问题是当left == up时,code[j] = code[j-1] = 0,然后再改所有的标号

当必须只能有一个回路的时候,Encode()一定要等于0.

#include <bits/stdc++.h>
#define maxn 20
using namespace std;
typedef long long ll;
int n, m; int need[maxn][maxn], a[maxn][maxn]; int cur = 0, Cnt, cnt; const int md = 997; struct Hashmap{
#define M 2000010
ll st[M];
int f[M], h[1000], nxt[M], size;
void init(){
memset(h, 0, sizeof h);
size = 0;
}
void push(ll hs, int val){
int tmp = hs % md;
for(int i = h[tmp]; i; i = nxt[i]){
if(st[i] == hs){
f[i] = max(f[i], val);
return;
}
}
int now = ++ size;
f[now] = val;
nxt[now] = h[tmp];
st[now] = hs;
h[tmp] = now;
}
}dp[2]; int code[maxn], ch[maxn]; void Decode(ll st){
for(int i = m; i >= 0; i --)
code[i] = st & 7, st >>= 3;
} ll Encode(){
ll ret = 0;
memset(ch, -1, sizeof ch);
ch[0] = 0; int cnt = 0;
for(int i = 0; i <= m; i ++){
if(ch[code[i]] == -1)
ch[code[i]] = ++ cnt;
code[i] = ch[code[i]];
ret <<= 3;
ret |= code[i];
}return ret;
} void Shift(){
for(int i = m; i; i --)
code[i] = code[i-1];
code[0] = 0;
} bool Judge(int cnt, int i, int j){
if(cnt && need[i][j] == 1)return false;
if(!cnt && need[i][j] == 2)return false;
return true;
} int ans = -0x7fffffff; void Trans(int i, int j){
dp[cur^1].init();
for(int k = 1; k <= dp[cur].size; k ++){
Decode(dp[cur].st[k]);
if(j == 1){if(code[m])continue;Shift();}
int lf = code[j-1], up = code[j], sta = 0;
for(int p = 0; p < j-1; p ++)
sta ^= (code[p] != 0);
if(lf && up){
code[j] = code[j-1] = 0;
if(lf == up){
if(Encode() == 0 && cnt == Cnt)//here!
ans = max(ans, dp[cur].f[k]);
}
else{
if(Judge(sta, i, j)){
for(int p = 0; p <= m; p ++)
if(code[p] == up)
code[p] = lf;
dp[cur^1].push(Encode(), dp[cur].f[k] + sta * a[i][j]);
}
}
}
else if(lf || up){
int tmp = lf ? lf : up;
if(Judge(sta, i, j)){
code[j-1] = 0, code[j] = tmp;
dp[cur^1].push(Encode(), dp[cur].f[k] + sta * a[i][j]);
}
sta ^= 1;
if(Judge(sta, i, j)){
code[j] = 0, code[j-1] = tmp;
dp[cur^1].push(Encode(), dp[cur].f[k] + sta * a[i][j]);
}
}
else{
if(Judge(sta, i, j))
dp[cur^1].push(Encode(), dp[cur].f[k] + sta * a[i][j]);
sta ^= 1;
if(Judge(sta, i, j)){
code[j] = code[j-1] = 15;
dp[cur^1].push(Encode(), dp[cur].f[k] + sta * a[i][j]);
}
}
}
if(need[i][j] == 2)cnt ++;
cur ^= 1;
} int main(){
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; i ++)
for(int j = 1; j <= m; j ++)
scanf("%d", &a[i][j]);
n ++, m ++;
dp[cur].init();
dp[cur].push(0, 0);
for(int i = 1; i <= n; i ++)
for(int j = 1; j <= m; j ++)
Trans(i, j);
printf("%d\n", ans); for(int i = 1; i < n; i ++)
for(int j = 1; j < m; j ++)
scanf("%d", &need[i][j]), Cnt += need[i][j] == 2;
ans = -0x7fffffff;
dp[cur].init();
dp[cur].push(0, 0);
for(int i = 1; i <= n; i ++)
for(int j = 1; j <= m; j ++)
Trans(i, j);
if(ans < -0x7ffffff)printf("Can not establish GFW.");
else printf("%d\n", ans);
return 0;
}

[插头DP自我总结]的更多相关文章

  1. 【BZOJ2310】ParkII 插头DP

    [BZOJ2310]ParkII Description Hnoi2007-Day1有一道题目 Park:给你一个 m * n 的矩阵,每个矩阵内有个权值V(i,j) (可能为负数),要求找一条回路, ...

  2. 【Learning】插头DP

    简介 插头DP(轮廓线DP)是用来解决网格图回路问题的一种算法. 插头DP解决的经典问题就是统计经过所有格子的哈密顿回路条数,某些格子有障碍. ​ 如果问题稍微进阶一点的话,不一定要求路径是回路.路径 ...

  3. 插头dp

    插头dp 感受: 我觉得重点是理解,算法并不是直接想出怎样由一种方案变成另一种方案.而是方案本来就在那里,我们只是枚举状态统计了答案. 看看cdq的讲义什么的,一开始可能觉得状态很多,但其实灰常简单 ...

  4. HDU 4113 Construct the Great Wall(插头dp)

    好久没做插头dp的样子,一开始以为这题是插头,状压,插头,状压,插头,状压,插头,状压,无限对又错. 昨天看到的这题. 百度之后发现没有人发题解,hust也没,hdu也没discuss...在acm- ...

  5. HDU 4949 Light(插头dp、位运算)

    比赛的时候没看题,赛后看题觉得比赛看到应该可以敲的,敲了之后发现还真就会卡题.. 因为写完之后,无限TLE... 直到后来用位运算代替了我插头dp常用的decode.encode.shift三个函数以 ...

  6. 插头DP专题

    建议入门的人先看cd琦的<基于连通性状态压缩的动态规划问题>.事半功倍. 插头DP其实是比较久以前听说的一个东西,当初是水了几道水题,最近打算温习一下,顺便看下能否入门之类. 插头DP建议 ...

  7. HDU 1693 Eat the Trees(插头DP、棋盘哈密顿回路数)+ URAL 1519 Formula 1(插头DP、棋盘哈密顿单回路数)

    插头DP基础题的样子...输入N,M<=11,以及N*M的01矩阵,0(1)表示有(无)障碍物.输出哈密顿回路(可以多回路)方案数... 看了个ppt,画了下图...感觉还是挺有效的... 参考 ...

  8. HDU 1693 Eat the Trees(插头DP)

    题目链接 USACO 第6章,第一题是一个插头DP,无奈啊.从头看起,看了好久的陈丹琦的论文,表示木看懂... 大体知道思路之后,还是无法实现代码.. 此题是插头DP最最简单的一个,在一个n*m的棋盘 ...

  9. HDU 4064 Carcassonne(插头DP)(The 36th ACM/ICPC Asia Regional Fuzhou Site —— Online Contest)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4064 Problem Description Carcassonne is a tile-based ...

随机推荐

  1. BZOJ2904

    找了一个晚上的资料,拼凑出来这么一个东西: 1) 如果是完全平方数返回12) 如果可以表示成形如$x^2+y^2$的形式输出2.这要求该数质因数分解后形如$4k+3$的质因数次数都是偶数.3) 如果该 ...

  2. BZOJ 3540 realtime-update 解题

    分析一下题意,大约是给定一串牛,然后找到一个跨越距离最长的牛子串使得在这个范围内白牛和花牛一样多. 白牛可以任意涂成花牛. 既然"白牛可以任意涂成花牛",那么我们需要找到一个最长的 ...

  3. 工作中常用shell之ssh登陆不用输入"yes"

    ip="192.168.5.166"ssh $ip -o StrictHostKeyChecking=no           //ssh登陆不用输入"yes" ...

  4. Valid Perfect Square

    Given a positive integer num, write a function which returns True if num is a perfect square else Fa ...

  5. SharePoint 2010 隐藏快速启动栏之使用内容编辑器webpart

    SharePoint 2010 自带的webpart里有一个叫内容编辑,在媒体和内容分类里面: 将其添加到页面后效果: 点击用于添加新内容,此时注意Ribbon菜单中的变化: 这里可以看到,你可以插入 ...

  6. Android drawable的自动缩放

    今天在写程序时发现,一张图片被自动放大了,后来发现,这张图片放在了drawable-zh文件夹下,这个文件夹没有指定屏幕密度!于是将drawable-zh改为drawable-zh-nodpi,问题解 ...

  7. Java集合框架中List接口的简单使用

    Java集合框架可以简单的理解为一种放置对象的容器,和数学中的集合概念类似,Java中的集合可以存放一系列对象的引用,也可以看做是数组的提升,Java集合类是一种工具类,只有相同类型的对象引用才可以放 ...

  8. BSD学习(BSD系统的历史和目标)

    UNIX系统的历史 unix系统的发展历程大概经历以下几个阶段: 贝尔实验室(Bell Laboratories)阶段,该实验室发明了UNIX 加州大学伯克利分校(University of Cali ...

  9. BaseServlet方法分发

    BaseServlet.java package org.guangsoft.controller; import java.io.IOException; import java.lang.refl ...

  10. icon font字体图标字库汇总

    最近在研究icon font图标字库,找了一些比较好的在线字库.大都是开源的,而且各有特色! 阿里icon font字库 http://www.iconfont.cn/ 这个是阿里妈妈M2UX的一个i ...