1877: [SDOI2009]晨跑
1877: [SDOI2009]晨跑
Time Limit: 4 Sec Memory Limit: 64 MB
Submit: 2007 Solved: 1085
[Submit][Status][Discuss]
Description
Input
Output
Sample Input
1 2 1
1 3 1
2 4 1
3 4 1
4 5 1
4 6 1
2 5 5
3 6 6
5 7 1
6 7 1
Sample Output
HINT
对于30%的数据,N ≤ 20,M ≤ 120。
对于100%的数据,N ≤ 200,M ≤ 20000。
Source
【思路】
拆点。我们可以将每个路口拆成两个点(i)和(i+N)。由于Ai与Bi之间有长度为C的街道,则在(Ai+N)和(Bi)之间添加一条容量为1,费用为c的边。然后对于每个Ai,添加一条(Ai,Ai+n)的边,容量为1(保证每个路口仅仅经过一次),费用为0的边。最大流即可。
要注意的是1和N不算是路口,所以(1与n+1)(n与2n)之间的边容量为INF。
#include<cstdio>
#include<iostream>
#define R register
#define inf 0x3f3f3f3f
using namespace std;
int read(){
R int x=;bool f=;
R char ch=getchar();
while(ch<''||ch>''){if(ch=='-')f=;ch=getchar();}
while(ch>=''&&ch<=''){x=(x<<)+(x<<)+ch-'';ch=getchar();}
return f?x:-x;
}
const int N=,M=1e5+;
struct node{
int u,v,next,cap,cost;
}e[M];int tot=;
int n,m,S,T,ans1,ans2,head[N],pree[N],prev[N],flow[N],dis[N],q[M*];
bool vis[N];
void add(int x,int y,int cap,int cost){
e[++tot].u=x;e[tot].v=y;e[tot].cap=cap;e[tot].cost=cost;e[tot].next=head[x];head[x]=tot;
e[++tot].u=y;e[tot].v=x;e[tot].cap=;e[tot].cost=-cost;e[tot].next=head[y];head[y]=tot;
}
bool spfa(){
for(int i=S;i<=T;i++) vis[i]=,dis[i]=inf;
int h=,t=;
q[t]=S;dis[S]=;vis[S]=;
while(h!=t){
int x=q[++h];vis[x]=;
for(int i=head[x];i;i=e[i].next){
int v=e[i].v;
if(e[i].cap>&&dis[v]>dis[x]+e[i].cost){
dis[v]=dis[x]+e[i].cost;
pree[v]=i;
if(!vis[v]){
vis[v]=;
q[++t]=e[i].v;
}
}
}
}
return dis[T]<inf;
}
void augment(){
int x=inf;
for(int i=pree[T];i;i=pree[e[i].u]) x=min(x,e[i].cap);
ans1++;
for(int i=pree[T];i;i=pree[e[i].u]){
ans2+=x*e[i].cost;
e[i].cap-=x;
e[i^].cap+=x;
}
}
int main(){
n=read();m=read();S=;T=n<<;
for(int i=,x,y,z;i<=m;i++){
x=read();y=read();z=read();
add(x+n,y,,z);
}
for(int i=;i<n;i++) add(i,i+n,,);
add(,S+n,inf,);
add(n,T,inf,);
while(spfa()) augment();
printf("%d %d",ans1,ans2);
return ;
}
/*zjk'code 跑的比我快,存一下
#include<cstdio>
#include<iostream>
#define N 410
#define M 41000
#define inf 1000000000
using namespace std;
int head[N],dis[N],q[N*10],vis[N],fa[N],n,m,S,T,cnt=1,timer,ans;
struct node{
int v,pre,f,w;
};node e[M];
void add(int u,int v,int f,int w){
e[++cnt].v=v;e[cnt].f=f;e[cnt].w=w;e[cnt].pre=head[u];head[u]=cnt;
e[++cnt].v=u;e[cnt].f=0;e[cnt].w=-w;e[cnt].pre=head[v];head[v]=cnt;
}
bool spfa(){
for(int i=1;i<=T;i++)dis[i]=inf;
int h=0,t=1;
q[1]=S;dis[S]=0;vis[S]=1;
while(h<t){
int x=q[++h];vis[x]=0;
for(int i=head[x];i;i=e[i].pre){
int v=e[i].v;
if(dis[v]>dis[x]+e[i].w&&e[i].f){
dis[v]=dis[x]+e[i].w;
fa[v]=i;
if(!vis[v]){
q[++t]=v;
vis[v]=1;
}
}
}
}
return dis[T]!=inf;
}
void updata(){
int tmp=T;
while(tmp!=S){
int l=fa[tmp],v=e[l].v;
e[l].f--;e[l^1].f++;
tmp=e[l^1].v;
}
ans+=dis[T];
}
int main(){
scanf("%d%d",&n,&m);
S=0,T=n*2+1;
add(S,1,inf,0);add(n*2,T,inf,0);
for(int i=2;i<n;i++)add(i,i+n,1,0);
add(1,1+n,inf,0);add(n,n+n,inf,0);
for(int i=1;i<=m;i++){
int x,y,z;scanf("%d%d%d",&x,&y,&z);
add(x+n,y,1,z);
}
while(spfa()){
timer++;updata();
}
printf("%d %d",timer,ans);
return 0;
}*/
1877: [SDOI2009]晨跑的更多相关文章
- BZOJ 1877: [SDOI2009]晨跑 费用流
1877: [SDOI2009]晨跑 Description Elaxia最近迷恋上了空手道,他为自己设定了一套健身计划,比如俯卧撑.仰卧起坐等 等,不过到目前为止,他坚持下来的只有晨跑. 现在给出一 ...
- BZOJ 1877: [SDOI2009]晨跑( 最小费用最大流 )
裸的费用流...拆点, 流量限制为1, 最后的流量和费用即答案. ------------------------------------------------------------------- ...
- BZOJ 3438 小M的作物 & BZOJ 1877 [SDOI2009]晨跑
我由衷地为我的朋友高兴.哈哈,yian,当你nick name破百上千时,再打“蒟蒻”就会被打的. 好的,说正事吧.请注意,这还是题解.但我发现,网络流实在是太套路了(怪不得这两年几乎销声匿迹).我们 ...
- 【BZOJ】1877: [SDOI2009]晨跑(最小费用最大流)
http://www.lydsy.com/JudgeOnline/problem.php?id=1877 费用流做多了,此题就是一眼题. 拆点表示只能经过一次,容量为1,费用为0. 然后再连边即可,跑 ...
- bzoj 1877 [SDOI2009]晨跑(最小费用最大流)
Description Elaxia最近迷恋上了空手道,他为自己设定了一套健身计划,比如俯卧撑.仰卧起坐等 等,不过到目前为止,他坚持下来的只有晨跑. 现在给出一张学校附近的地图,这张地图中包含N个十 ...
- 1877. [SDOI2009]晨跑【费用流】
Description Elaxia最近迷恋上了空手道,他为自己设定了一套健身计划,比如俯卧撑.仰卧起坐等 等,不过到目前为止,他 坚持下来的只有晨跑. 现在给出一张学校附近的地图,这张地图中包含N个 ...
- BZOJ 1877 [SDOI2009]晨跑(多条不交叉最短路)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=1877 [题目大意] 找出最多有几条点不重复的从1到N的路,并且要求在满足这个条件的情况 ...
- bzoj 1877: [SDOI2009]晨跑
#include<cstdio> #include<iostream> #include<cstring> #define M 6009 #define inf 2 ...
- BZOJ 1877: [SDOI2009]晨跑(费用流)
看到要求两个量就下意识的想到了费用流= =,先把一个点拆成两个点就能够解决一个的只经过一次的限制 CODE: #include<cstdio>#include<iostream> ...
随机推荐
- 11gRAC报错CRS-4535, CRS-4000解决
环境:AIX6.1 + Oracle11.2.0.4 RAC(2 nodes) 1.故障现象 2.定位问题 3.处理问题 1.故障现象 使用crsctl查看集群各资源状态,在任一节点都会直接报错CRS ...
- code
using System;using System.Threading; namespace ThreadLocalTest{ public class MyObject { ...
- Markdown通用的常用语法说明
前言 Markdown 是一种轻量级的 标记语言,语法简洁明了.学习容易,还具有其他很多优点,目前被越来越多的人用来写作使用. Markdown具有一系列衍生版本,用于扩展Markdown的功能(如表 ...
- asp.net记录错误日志的方法
1.说明 在调试发布后的asp.net项目时有可能会遇到意想不到的错误,而未能及时的显示.这就需要记录日志来跟踪错误信息,所以写了个简单的记录信息的方法,记录简单的文本信息也可以使用.此方法是以生成文 ...
- Windows安装RabbitMQ集群的几个注意点
记录一下RabbitMQ在windows平台下安装的几个注意点- -,好记性不如烂笔头 安装过程与Linux安装一致,教程参照官网集群配置:此处只列举出几个注意点: 1. erlang的版本需要一致, ...
- JVM调优总结
堆大小设置JVM 中最大堆大小有三方面限制:相关操作系统的数据模型(32-bt还是64-bit)限制:系统的可用虚拟内存限制:系统的可用物理内存限制.32位系统下,一般限制在1.5G~2G:64为操作 ...
- SAP CRM 客户控制器与数据绑定
当用户从视图离开时,视图将失去它的数据.解决这个问题,需要引入客户控制器(Custom Controller)(译者注:SAP CRM客户端中,不同地方的Custom Controller会翻译为“客 ...
- [Android]Android MVP&依赖注入&单元测试
以下内容为原创,欢迎转载,转载请注明 来自天天博客:http://www.cnblogs.com/tiantianbyconan/p/5422443.html Android MVP&依赖注入 ...
- dSYM 文件分析工具
来到新公司后,前段时间就一直在忙,前不久 项目 终于成功发布上线了,最近就在给项目做优化,并排除一些线上软件的 bug,因为项目中使用了友盟统计,所以在友盟给出的错误信息统计中能比较方便的找出客户端异 ...
- Unable to simultaneously satisfy constraints.
在进行版本的迭代更新时,新功能需求需要对主页面的UI进行重新的布局,但是,报了错误,出了好多约束方面的问题: Unable to simultaneously satisfy constraints. ...