题目描述

  有一个\(n\)行\(m\)列的网格图。

  \(S\)到第一行的每一个点都有一条单向边,容量为\(\infty\)。

  最后一行的每个点到\(T\)都有一条单向边,容量为\(\infty\)。

  同一行中相邻的两个节点之间有一条无向边,\((x,y)\)和\((x,y+1)\)之间的无向边的容量为\(a_{x,y}\)。

  同一列中相邻的两个节点之间有两条有向边,\((x,y)\)到\((x+1,y)\)这条有向边的容量为\(b_{x,y}\),\((x+1,y)\)到\((x,y)\)这条有向边容量为\(\infty\)。

  求\(S\)到\(T\)的最大流。

  特别的,\(\forall i,a_{1,i}=a_{n,i}=0\)

  \(n\times m\leq 25000000\)

题解

  显然这是一个网络流。

  直接跑网络流会TLE。

  观察到这个图是一个平面图,可以把平面图网络流转化为对偶图最短路。

  怎么转化呢?

  首先你要会无向边的平面图网络流(可以百度/google)。

  有向边的连边方法和无向边的类似。

  对于一条有向边\(x\rightarrow y\),容量为\(z\)的有向边(网络流最后都是在有向边上面跑的),从\(x\rightarrow y\)这条有向边的左边对应的这个区域连一条边到右边的这个区域,权值为\(z\)。

  最后跑一次最短路就行了。

  这道题中从下往上的边的容量为\(\infty\),所以对偶图中从左往右的边的权值为\(\infty\),也就是说最短路的每一步只会向上/下/左走,这就可以DP了。

  设\(f_{i,j}\)为走到\((x,y)\)右下方那个区域的最短路

\[f_{i,j}=\min(f_{i-1,j}+a_{i,j},f_{i+1,j}+a_{i+1,j},f_{i,j+1}+b_{i+1,j})
\]

  时间复杂度:\(O(nm)\)

  我的代码中把左右反过来了

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<utility>
#include<cmath>
#include<functional>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
void sort(int &a,int &b)
{
if(a>b)
swap(a,b);
}
void open(const char *s)
{
#ifndef ONLINE_JUDGE
char str[100];
sprintf(str,"%s.in",s);
freopen(str,"r",stdin);
sprintf(str,"%s.out",s);
freopen(str,"w",stdout);
#endif
}
int rd()
{
int s=0,c;
while((c=getchar())<'0'||c>'9');
do
{
s=s*10+c-'0';
}
while((c=getchar())>='0'&&c<='9');
return s;
}
void put(int x)
{
if(!x)
{
putchar('0');
return;
}
static int c[20];
int t=0;
while(x)
{
c[++t]=x%10;
x/=10;
}
while(t)
putchar(c[t--]+'0');
}
int upmin(int &a,int b)
{
if(b<a)
{
a=b;
return 1;
}
return 0;
}
int upmax(int &a,int b)
{
if(b>a)
{
a=b;
return 1;
}
return 0;
}
int *a;
int n,m;
int A,B,Q;
int down(int x,int y)
{
return a[(x-1)*m+y];
}
int right(int x,int y)
{
return a[(n-1)*m+(x-2)*(m-1)+y];
}
ll *f;
int main()
{
open("c");
scanf("%d%d",&n,&m);
a=new int[(n-1)*m+(n-2)*(m-1)+1];
scanf("%d%d%d%d",&A,&B,&Q,&a[0]);
for(int i=1;i<=(n-1)*m+(n-2)*(m-1);i++)
a[i]=((ll)a[i-1]*A+B)%Q;
f=new ll[n];
for(int i=1;i<n;i++)
f[i]=0;
for(int i=1;i<m;i++)
{
for(int j=1;j<n;j++)
f[j]+=down(j,i);
for(int j=2;j<n;j++)
f[j]=min(f[j],f[j-1]+right(j,i));
for(int j=n-2;j>=1;j--)
f[j]=min(f[j],f[j+1]+right(j+1,i));
}
ll ans=0x7fffffffffffffffll;
for(int i=1;i<n;i++)
ans=min(ans,f[i]+down(i,m));
printf("%lld\n",ans);
return 0;
}

【XSY2849】陈姚班 平面图网络流 最短路 DP的更多相关文章

  1. 省队集训 Day3 陈姚班

    [题目大意] 给一张网格图,上往下有流量限制,下往上没有,左往右有流量限制. $n * m \leq 2.5 * 10^6$ [题解] 考场直接上最大流,50分.竟然傻逼没看出狼抓兔子. 平面图转对偶 ...

  2. 200万年薪请不到!清华姚班到底有多牛X?

    前几天,清华大学自动化系2020年大一新生的C++作业因为太难而上了热搜,该话题在知乎上的热度一度高达 1300+ 万.  在该帖子下方,有很多关于这件事的讨论,其中很多不禁赞叹"清华太牛 ...

  3. 中国 AI 天才养成计划:清华姚班和 100 个「张小龙」

    https://daily.zhihu.com/story/9653612?from=timeline&isappinstalled=0   AI财经社,专注未来,以及更好的生活 真正的 AI ...

  4. zz姚班天才少年鬲融凭非凸优化研究成果获得斯隆研究奖

    姚班天才少年鬲融凭非凸优化研究成果获得斯隆研究奖 近日,美国艾尔弗·斯隆基金会(The Alfred P. Sloan Foundation)公布了2019年斯隆研究奖(Sloan Research ...

  5. 旷视6号员工范浩强:高二开始实习,“兼职”读姚班,25岁在CVPR斩获第四个世界第一...

    初来乍到,这个人说话容易让人觉得"狂". "我们将比赛结果提交上去,果不其然,是第一名的成绩."当他说出这句话的时候,表情没有一丝波澜,仿佛一切顺理成章. 他说 ...

  6. hdu 4568 Hunter 最短路+dp

    Hunter Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Subm ...

  7. BZOJ_1003_[ZJOI2006]物流运输_最短路+dp

    BZOJ_1003_[ZJOI2006]物流运输_最短路+dp 题意:http://www.lydsy.com/JudgeOnline/problem.php?id=1003 分析: 这种一段一段的显 ...

  8. POJ 3635 Full Tank? 【分层图/最短路dp】

    任意门:http://poj.org/problem?id=3635 Full Tank? Time Limit: 1000MS   Memory Limit: 65536K Total Submis ...

  9. [USACO07NOV]牛继电器Cow Relays (最短路,DP)

    题目链接 Solution 非正解 似乎比较蛇啊,先个一个部分分做法,最短路+\(DP\). 在求最短路的堆或者队列中存储元素 \(dis_{i,j}\) 代表 \(i\) 这个节点,走了 \(j\) ...

随机推荐

  1. 钉钉JSAPI前端鉴权

    钉钉二次开发分为如下表所示三种类型的开发,只有企业内部应用才需要对JSAPI鉴权. 类型 开发方式 JSAPI鉴权 应用场景 第三方企业应用 E应用开发 不需要 用于发布到钉钉应用市场,供广大用户下载 ...

  2. JMeter 内置日期(时间)函数总结

    JMeter 内置日期(时间)函数总结   by:授客 QQ:1033553122 1. 测试环境 apache-jmeter-3.3 下载地址: http://jmeter.apache.org/c ...

  3. MongoDB副本集功能及节点属性梳理

    副本集的主要功能 副本集是MongoDB高可用的基础,其主要作用 归纳为以下几点: (1)高可用,防止设备(服务器.网络)故障.提供自动FailOver功能. (2)无需配置高可用性虚拟节点:无论是S ...

  4. Can't create/write to file '/tmp/MLjnvU95' (Errcode: 13 - Permission denied)

    今天一个同事反馈往一个MySQL数据库导入数据时,报"ERROR 1 (HY000): Can't create/write to file '/tmp/MLjnvU95' (Errcode ...

  5. LNMP时,出现502 Bad Gateway的错误提示

    因为工作需要,要在ubuntu中安装LNMP环境,在这里,php是最新版本php7.1.一切都进展得很顺利,安装完成后,在浏览器中输入http://127.0.0.1/info.php,出现了502 ...

  6. SQLServer之创建非聚集索引

    开始之前 典型实现 可以通过下列方法实现非聚集索引: UNIQUE 约束 在创建 UNIQUE 约束时,默认情况下将创建唯一非聚集索引,以便强制 UNIQUE 约束. 如果不存在该表的聚集索引,则可以 ...

  7. isinstance_issubclass

    isinstance和issubclass分别是检验是不是对象是不是类的实例化,和子类是不是这个父类的子类 class A:pass class B(A):pass a = A() print(isi ...

  8. Python爬虫【解析库之beautifulsoup】

    解析库的安装 pip3 install beautifulsoup4 初始化 BeautifulSoup(str,"解析库") from bs4 import BeautifulS ...

  9. Java注解开发与应用案例

    Java注解开发与应用案例 Annotation(注解)是JDK5.0及以后版本引入的,可以对包.类.属性.方法的描述,给被述对象打上标签,被打上标签后的类.属性.方法将被赋予特殊的“功能”:打个比喻 ...

  10. git添加/删除远程仓库

    注意:仓库只有管理员建的你才有权限上传,不然自己建的也没用,没权限上传 1.远程仓库路径查询 git remote -v 2.添加远程仓库 git remote add origin <你的项目 ...