题目描述
这是在阿尔托利亚·潘德拉贡成为英灵前的事情,她正要去拔出石中剑成为亚瑟王,在这之前她要去收集一些宝石。
宝石排列在一个n*m的网格中,每个网格中有一块价值为v(i,j)的宝石,阿尔托利亚·潘德拉贡可以选择自己的起点。
开始时刻为0秒。以下操作,每秒按顺序执行
- 在第i秒开始的时候,阿尔托利亚·潘德拉贡在方格(x,y)上,她可以拿走(x,y)中的宝石。
- 在偶数秒,阿尔托利亚·潘德拉贡周围四格的宝石会消失
- 若阿尔托利亚·潘德拉贡第i秒开始时在方格(x,y)上,则在第i+1秒可以立即移动到(x+1,y),(x,y+1),(x-1,y)或(x,y-1)上,也可以停留在(x,y)上。
求阿尔托利亚·潘德拉贡最多可以获得多少价值的宝石
输入格式
第一行给出数字N,M代表行列数.N,M均小于等于100,宝石的价值不会超过10000.下面N行M列用于描述数字矩阵
输出格式
输出最多可以拿到多少价值宝石
输入输出样例
输入 #1
2 2 1 2 2 1
输出 #1
4
说明/提示
姚金宇的原创题。
代码:
include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define R register
#define inf 1e9+7
using namespace std;
const int Excalibur = 10005;//N或MAXN
int n, m, S, T, cur[Excalibur], dep[Excalibur];
struct saber {//edge,我的王
int nxt, to, v;
}rin[Excalibur<<3];//是凛哦
int Lancer[Excalibur], Fate = 1, Archar;//emmm,head和tot以及sum
int x[5] = {0, -1, 0, 1, 0}, y[5] = {0, 0, -1, 0, 1};
inline void add(int from, int to, int v) {
rin[++Fate].v = v;
rin[Fate].to = to;
rin[Fate].nxt = Lancer[from];
Lancer[from] = Fate;
}
inline bool bfs(int s, int t) {
queue<int> Rider;//q
for(R int i = 0;i <= T;++ i) cur[i] = Lancer[i], dep[i] = 0;
Rider.push(s); dep[s] = 1;
while(!Rider.empty()) {
R int vi = Rider.front();
Rider.pop();
for(R int i = Lancer[vi]; i ;i = rin[i].nxt) {
R int vc = rin[i].to;
if(!dep[vc] && rin[i].v) {
dep[vc] = dep[vi] + 1;
Rider.push(vc);
}
}
}
return dep[t];
}
int dfs(int s, int t, int flow) {
if(!flow || s == t) return flow;
int Caster = 0, Assassin;//增量
for(R int i = cur[s]; i ;i = rin[i].nxt) {
R int vc = rin[i].to;
cur[s] = i;
if(dep[vc] == dep[s] + 1 && rin[i].v) {
Assassin = dfs(vc, t, min(flow, rin[i].v));
if(!Assassin) continue;
Caster += Assassin; flow -= Assassin;
rin[i].v -= Assassin; rin[i ^ 1].v += Assassin;
if(!flow) return Caster;
}
}
return Caster;
}
int Dinic() {
int res = 0;
while(bfs(S, T)) res += dfs(S, T, inf);
return Archar - res;
}
int main() {
scanf("%d%d",&n, &m);
//以下是建图,中二到此结束
S = 0, T = n * m + 1;
for(R int i = 1;i <= n;++ i)
for(R int j = 1;j <= m;++ j) {
R int v; scanf("%d",&v);
Archar += v;
R int p = (i - 1) * m + j;
if((i + j) & 1) { add(S, p, v);add(p, S, 0); }
else { add(p, T, v);add(T, p, 0); }
}
for(R int i = 1;i <= n;++ i) for(R int j = 1;j <= m;++ j)
if((i + j) & 1)
for(R int k = 1;k <= 4;++ k) {
R int f = i + x[k], g = j + y[k];
if(f < 1 || f > n || g < 1 || g > m) continue;
R int u = (i - 1) * m + j, v = (f - 1) * m + g;
add(u, v, inf); add(v, u, 0);
}
printf("%d",Dinic());
return 0;
}