题意分析

首先 要求起点终点不连通

再结合数据范围 就是最小割了

首先我们可以建一个图出来

如果\(x\)可以到\(y\)的话

那么我们就从\(x\)向\(y\)连一条代价为\(h[x]-h[y]+1\)的边 代表不联通的代价

可是如果存在以下情况呢

如果我们选择切断\(c\)到\(d\)的边的话

实际上我们也切断了\(a\)到\(c\)以及\(b\)到\(c\)

所以我们可以这么建

然后跑最大流就可以了

由于起点以及终点不可以被修改

所以我忽视了别的点到起点的连边

同时别的点到终点的连边边权都是\(inf\)

CODE:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#include<string>
#include<queue>
#include<map>
#include<stack>
#include<list>
#include<set>
#include<deque>
#include<vector>
#include<ctime>
#define ll long long
#define inf 0x7fffffff
#define N 88
#define IL inline
#define M 5008611
#define D double
#define ull unsigned long long
#define R register
using namespace std;
template<typename T>IL void read(T &_)
{
T __=0,___=1;char ____=getchar();
while(!isdigit(____)) {if(____=='-') ___=0;____=getchar();}
while(isdigit(____)) {__=(__<<1)+(__<<3)+____-'0';____=getchar();}
_=___ ? __:-__;
}
/*-------------OI使我快乐-------------*/
int T;ll ans;
int n,m,tot=1,sx,sy,tx,ty;
int hei[N][N],to[M],nex[M],head[N*N*5],w[M];
int dep[5*N*N],cur[5*N*N];
vector<pair<int,int> > G[N*N];
queue<int> Q;
bool vis[N][N];
int tox[6]={0,0,0,1,-1},toy[6]={0,1,-1,0,0};
IL int id(int x,int y){return (x-1)*m+y;}
IL bool safe(int x,int y){return x>=1&&x<=n&&y>=1&&y<=m;}
IL void add(int x,int y,int z)
{to[++tot]=y;nex[tot]=head[x];head[x]=tot;w[tot]=z;
swap(x,y);to[++tot]=y;nex[tot]=head[x];head[x]=tot;w[tot]=0;}
IL bool bfs()
{
for(R int i=1;i<=5*n*m;++i) dep[i]=0;
Q.push(id(sx,sy));dep[id(sx,sy)]=1;
for(;!Q.empty();)
{
int u=Q.front();Q.pop();
for(R int i=head[u];i;i=nex[i])
{
int v=to[i];
if(w[i]>0&&dep[v]==0)
{
dep[v]=dep[u]+1;Q.push(v);
}
}
}
return dep[id(tx,ty)]!=0;
}
IL int dfs(int now,int res)
{
if(now==id(tx,ty)||res==0) return res;
for(R int &i=cur[now];i;i=nex[i])
{
int v=to[i];
if(w[i]>0&&dep[v]==dep[now]+1)
{
int have=dfs(v,min(w[i],res));
if(have>0)
{
w[i]-=have;w[i^1]+=have;
return have;
}
}
}
return 0;
}
IL void Dinic()
{
while(bfs())
{
for(R int i=1;i<=5*n*m;++i) cur[i]=head[i];
int d=dfs(id(sx,sy),inf);
while(d) ans+=d,d=dfs(id(sx,sy),inf);
}
}
int main()
{
freopen("ski.in","r",stdin);
freopen("ski.out","w",stdout);
read(T);
while(T--)
{
read(n);read(m);read(sx);read(sy);read(tx);read(ty);
tot=1;ans=0;memset(head,0,sizeof head);
for(R int i=1;i<=n;++i)
for(R int j=1;j<=m;++j)
read(hei[i][j]);
if(abs(sx-tx)+abs(sy-ty)==1) {puts("-1");continue;}
for(R int i=1;i<=n;++i)
{
for(R int j=1;j<=m;++j)
for(R int k=1;k<=4;++k)
{
int nowx=i+tox[k],nowy=j+toy[k];
if(nowx==sx&&nowy==sy) continue;
if(safe(nowx,nowy)&&hei[nowx][nowy]<=hei[i][j])
G[id(nowx,nowy)].push_back(make_pair(hei[i][j]-hei[nowx][nowy]+1,id(i,j)));
}
}
for(R int i=1;i<=n*m;++i)
{
if(G[i].size())
sort(G[i].begin(),G[i].end());
}
// for(R int i=1;i<=n*m;++i)
// {
// if(G[i].size())
// {
// for(R int j=0;j<(int)G[i].size();++j)
// printf("%d -- %d --> %d\n",G[i][j].second,G[i][j].first,i);
// }
// }
for(R int i=1;i<=n*m;++i)
{
int now=0;
if(id(tx,ty)==i)
{
for(R int j=(int)G[i].size()-1;j>=0;--j)
{
// printf("%d -- %d --> %d\n",G[i][j].second,G[i][j].first,i);
++now;
add(G[i][j].second,i+now*n*m,inf);
add(i+now*n*m,i+(now-1)*n*m,inf);
}
}
else
{
for(R int j=(int)G[i].size()-1;j>=0;--j)
{
// printf("%d -- %d --> %d\n",G[i][j].second,G[i][j].first,i);
++now;
add(G[i][j].second,i+now*n*m,inf);
add(i+now*n*m,i+(now-1)*n*m,G[i][j].first);
}
}
}
Dinic();
printf("%lld\n",ans);
for(R int i=1;i<=n*m;++i) G[i].clear();
}
fclose(stdin);
fclose(stdout);
return 0;
}

HEOI 2019 RP++

考试题 T2的更多相关文章

  1. hdu3089 Josephus again|快速约瑟夫环

    题目链接:戳我 貌似是高一昨天的考试题T2?????感觉挺好玩的就搞了搞qwqwq 其实是HDU上面的题啦.... 对于普通的约瑟夫问题,大概是n个人围成一个环,从1开始报数,数到k,那个人出队,最后 ...

  2. [Noip2016]蚯蚓 D2 T2 队列

    [Noip2016]蚯蚓 D2 T2 Description 本题中,我们将用符号[c]表示对c向下取整,例如:[3.0」= [3.1」=[3.9」=3.蛐蛐国最近蚯蚓成灾了!隔壁跳 蚤国的跳蚤也拿蚯 ...

  3. T2 Func<in T1,out T2>(T1 arg)

    委托调用方法的4种方式. using System; using System.Collections.Generic; namespace ConsoleApplication1 { delegat ...

  4. Hotelling T2检验和多元方差分析

    1.1 Hotelling T2检验 Hotelling T2检验是一种常用多变量检验方法,是单变量检验的自然推广,常用于两组均向量的比较. 设两个含量分析为n,m的样本来自具有公共协方差阵的q维正态 ...

  5. bzoj4034: [HAOI2015]T2

    4034: [HAOI2015]T2 Time Limit: 10 Sec  Memory Limit: 256 MB Submit: 2684  Solved: 843 Description 有一 ...

  6. 【BZOJ 4517】【SDOI 2016 Round1 Day2 T2】排列计数

    本蒟蒻第一次没看题解A的题竟然是省选$Round1$ $Day2$ $T2$ 这道组合数学题. 考试时一开始以为是莫队,后来想到自己不会组合数的一些公式,便弃疗了去做第三题,,, 做完第三题后再回来看 ...

  7. NOIP欢乐模拟赛 T2 解题报告

    小澳的坐标系 (coordinate.cpp/c/pas) [题目描述] 小澳者表也,数学者景也,表动则景随矣. 小澳不喜欢数学,可数学却待小澳如初恋,小澳睡觉的时候也不放过. 小澳的梦境中出现了一个 ...

  8. Action<T1, T2>委托

    封装包含两个参数的方法委托,没有返回值. 语法 public delegate void Action<in T1, in T2>( T1 arg1, T2 arg2 ) 类型参数 in ...

  9. NOIP2013普及组 T2 表达式求值

    OJ地址:洛谷P1981 CODEVS 3292 正常写法是用栈 #include<iostream> #include<algorithm> #include<cmat ...

随机推荐

  1. 转)bash快捷键

    粗体表示推荐,也许对每个人不同. Ctrl-A 相当于HOME键,用于将光标定位到本行最前面 Ctrl-E 相当于End键,即将光标移动到本行末尾 Ctrl-B 相当于左箭头键,用于将光标向左移动一格 ...

  2. golang之指针

    先上代码 package main import "fmt" type myInt int //匿名函数 //值的接收 //指针的接收 //*myint表示myInt的指针类型 / ...

  3. C语言中的undefined behavior

    参考: http://www.cnblogs.com/aoaoblogs/archive/2010/08/31/1813982.html

  4. code1105 过河

    dp方程很简单: f[i] = min{ f[i-j] } + stone[i] 但是数据10^9太大了,超时超空间,这样只能过30% 来自:http://blog.csdn.net/w1996070 ...

  5. 在Linux下使用logrotate管理日志(转)

    原文地址:http://www.tuicool.com/articles/ieAnMjN logrotate是日志循环管理工具,可以分割日志文件,删除旧的日志文件,创建新的日志文件,循环管理日志从而节 ...

  6. 安装gcc及其依赖

    在gcc-4.8.2和gcc-4.1.2基础上编译gcc-5.2.0,有可能会遇到一些问题. 要想成功编译gcc,则在编译之前需要安装好它的至少以下三个依赖: gmp mpfr mpc 而mpc又依赖 ...

  7. Linux 基础教程 37-进程命令

    pidof     我们知道每个小孩一出生就会一个全国唯一的编号来对其进行标识,用于以后上学,办社保等,就是我们的身份证号.那么在Linux系统中,用来管理运行程序的标识叫做PID,就是大家熟知的进程 ...

  8. Tweak和app交互方案【进程通信】

    Core Foundation DEMO:Tweak端: CFNotificationCenterAddObserver(CFNotificationCenterGetDarwinNotifyCent ...

  9. JSP和servlet之间的传值(总结的很全面)

    转自:http://blog.csdn.net/ssy_shandong/article/details/9328985 1.从一个jsp页面跳转到另一个jsp页面时的参数传递     (1)使用re ...

  10. 在DOS命令行窗口中显示系统环境环境变量

    (这是一个小技巧) 示例命令: echo %path% path是系统环境变量,使用百分号包围起来 http://www.cnblogs.com/danzhang 张洪君 微软ALM MVP