1266: [AHOI2006]上学路线route

Time Limit: 3 Sec  Memory Limit: 162 MB
Submit: 2356  Solved: 841
[Submit][Status][Discuss]

Description

可可和卡卡家住合肥市的东郊,每天上学他们都要转车多次才能到达市区西端的学校。直到有一天他们两人参加了学校的信息学奥林匹克竞赛小组才发现每天上学的乘车路线不一定是最优的。 可可:“很可能我们在上学的路途上浪费了大量的时间,让我们写一个程序来计算上学需要的最少时间吧!” 合肥市一共设有N个公交车站,不妨将它们编号为1…N的自然数,并认为可可和卡卡家住在1号汽车站附近,而他们学校在N号汽车站。市内有M条直达汽车路线,执行第i条路线的公交车往返于站点pi和qi之间,从起点到终点需要花费的时间为ti。(1<=i<=M, 1<=pi, qi<=N) 两个人坐在电脑前,根据上面的信息很快就编程算出了最优的乘车方案。然而可可忽然有了一个鬼点子,他想趁卡卡不备,在卡卡的输入数据中删去一些路线,从而让卡卡的程序得出的答案大于实际的最短时间。而对于每一条路线i事实上都有一个代价ci:删去路线的ci越大卡卡就越容易发现这个玩笑,可可想知道什么样的删除方案可以达到他的目的而让被删除的公交车路线ci之和最小。 [任务] 编写一个程序:  从输入文件中读取合肥市公交路线的信息;  计算出实际上可可和卡卡上学需要花费的最少时间;  帮助可可设计一个方案,删除输入信息中的一些公交路线,使得删除后从家到学校需要的最少时间变大,而被删除路线的ci和最小;向输出文件输出答案。

Input

输入文件中第一行有两个正整数N和M,分别表示合肥市公交车站和公交汽车路线的个数。以下M行,每行(第i行,总第(i+1)行)用四个正整数描述第i条路线:pi, qi, ti, ci;具体含义见上文描述。

Output

输出文件最多有两行。 第一行中仅有一个整数,表示从可可和卡卡家到学校需要的最短时间。 第二行输出一个整数C,表示Ci之和

Sample Input

6 7
1 2 1 3
2 6 1 5
1 3 1 1
3 4 1 1
4 6 1 1
5 6 1 2
1 5 1 4

Sample Output

2
5

HINT

2<=N<=500, 1<=M<=124 750, 1<=ti, ci<=10 000
合肥市的公交网络十分发达,你可以认为任意两个车站间都可以通过直达或转车互相到达,当然如果在你提供的删除方案中,家和学校无法互相到达,那么则认为上学需要的最短为正无穷大:这显然是一个合法的方案。

Source

思路:跑spfa+网络流最小割。跑正反两边spfa,然后拎出一条边,判断边权+左端点到1号节点的最短路+右端点到n号节点的最短路(离1号节点近的点为左端点,反之为右端点)是否等于1到n的最短路,如果相等左右端点连边,权值为这条边代价。然后跑dinic。

错因:网络流建图错误,不应该是跑spfa时建图,因为spfa有松弛操作,可能边数建多了。

#include<queue>
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define N 501
#define M 124751
using namespace std;
queue<int>que;
struct nond{
int x,y,z,zz;
}edge[M];
int n,m;;
int tot,tot1=;
int scr,decc,ans,maxn;
int lev[M],cur[M];
int dis[N],dis1[N],vis[N];
int to1[M*],net1[M*],cap1[M*],head1[M*];
int to[M*],net[M*],cap[M*],val[M*],head[M*];
void add1(int u,int v,int w){
to1[++tot1]=v;net1[tot1]=head1[u];cap1[tot1]=w;head1[u]=tot1;
to1[++tot1]=u;net1[tot1]=head1[v];cap1[tot1]=;head1[v]=tot1;
}
void add(int u,int v,int w,int z){
to[++tot]=v;net[tot]=head[u];cap[tot]=w;val[tot]=z;head[u]=tot;
to[++tot]=u;net[tot]=head[v];cap[tot]=w;val[tot]=z;head[v]=tot;
}
void spfa(int s){
memset(vis,,sizeof(vis));
memset(dis,0x3f,sizeof(dis));
while(!que.empty()) que.pop();
que.push(s);
vis[s]=;dis[s]=;
while(!que.empty()){
int now=que.front();
que.pop();
vis[now]=;
for(int i=head[now];i;i=net[i])
if(dis[to[i]]>dis[now]+cap[i]){
dis[to[i]]=dis[now]+cap[i];
if(!vis[to[i]]){
vis[to[i]]=;
que.push(to[i]);
}
}
}
}
void spfa1(int s){
memset(vis,,sizeof(vis));
memset(dis1,0x3f,sizeof(dis));
while(!que.empty()) que.pop();
que.push(s);
vis[s]=;dis1[s]=;
while(!que.empty()){
int now=que.front();
que.pop();
vis[now]=;
for(int i=head[now];i;i=net[i])
if(dis1[to[i]]>dis1[now]+cap[i]){
dis1[to[i]]=dis1[now]+cap[i];
if(!vis[to[i]]){
vis[to[i]]=;
que.push(to[i]);
}
}
}
}
bool bfs(){
for(int i=scr;i<=decc;i++){
lev[i]=-;
cur[i]=head1[i];
}
while(!que.empty()) que.pop();
que.push(scr);
lev[scr]=;
while(!que.empty()){
int now=que.front();
que.pop();
for(int i=head1[now];i;i=net1[i])
if(lev[to1[i]]==-&&cap1[i]>){
lev[to1[i]]=lev[now]+;
que.push(to1[i]);
if(to1[i]==decc) return true;
}
}
return false;
}
int dinic(int now,int flow){
if(now==decc) return flow;
int rest=,detal;
for(int & i=cur[now];i;i=net1[i])
if(lev[to1[i]]==lev[now]+&&cap1[i]){
detal=dinic(to1[i],min(flow-rest,cap1[i]));
if(detal){
cap1[i]-=detal;
cap1[i^]+=detal;
rest+=detal;
if(rest==flow) break;
}
}
if(rest!=flow) lev[now]=-;
return rest;
}
int main(){
scanf("%d%d",&n,&m);
scr=;decc=n;
for(int i=;i<=m;i++){
int u,v,w,z;
scanf("%d%d%d%d",&u,&v,&w,&z);
edge[i].x=u;
edge[i].y=v;
edge[i].z=w;
edge[i].zz=z;
add(u,v,w,z);
}
spfa();
cout<<dis[n]<<endl;
spfa1(n);
for(int i=;i<=m;i++){
int u=edge[i].x,v=edge[i].y;
if(dis[u]>dis[v]) swap(u,v);
if(dis[u]+dis1[v]+edge[i].z==dis[n])
add1(u,v,edge[i].zz);
}
while(bfs())
ans+=dinic(scr,0x7f7f7f7f);
cout<<ans;
}

最新文章

  1. Liferay7 BPM门户开发之46: 集成Activiti用户、用户组、成员关系同步
  2. Servlet作业2-将表单提交的商品信息输出到页面中
  3. 首个攻击该Mac OS系统的恶意软件——KeRanger
  4. LightOJ1122 Digit Count(DP)
  5. 使用StarUML创建类图
  6. struts2漏洞原理及解决办法
  7. VS调试Libevent流程
  8. 代码片段 - JavaScript 求时间差
  9. SPRING 配置文件和类名
  10. 根据input 标签取value属性的值
  11. Sql Server中三种字符串合并方法的性能比较
  12. HDU5918(KMP)
  13. MySQL InnoDB下关于MVCC的一个问题的分析
  14. rest service下载文件与上传
  15. So you want to be a computational biologist?
  16. dubbo-001--前言
  17. 287find-the-duplicate-number
  18. js里面声明变量时候的注意事项
  19. css用法大全
  20. 20155230 实验二《Java面向对象程序设计》实验报告

热门文章

  1. android 经典博客
  2. python时间戳
  3. 从发请求到AJAX到同源政策
  4. Codeforces Round #445
  5. 9) 十分钟学会android--使用Fragment建立动态UI
  6. 6) 十分钟学会android--Activity的生命周期之启动与销毁
  7. Windows Phone - 按钮/button 控件
  8. 【从零开始】【Java】【0】装软件些
  9. 关于《Python核心编程》第2版和第3版
  10. jq遍历 下拉框