bzoj 1266 1266: [AHOI2006]上學路線route

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=1;
int scr,decc,ans,maxn;
int lev[M],cur[M];
int dis[N],dis1[N],vis[N];
int to1[M*2],net1[M*2],cap1[M*2],head1[M*2];
int to[M*2],net[M*2],cap[M*2],val[M*2],head[M*2];
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]=0;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,0,sizeof(vis));
    memset(dis,0x3f,sizeof(dis));
    while(!que.empty())    que.pop();
    que.push(s);
    vis[s]=1;dis[s]=0;
    while(!que.empty()){
        int now=que.front();
        que.pop();
        vis[now]=0;
        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]]=1;
                    que.push(to[i]);
                }
            }
    }
}
void spfa1(int s){
    memset(vis,0,sizeof(vis));
    memset(dis1,0x3f,sizeof(dis));
    while(!que.empty())    que.pop();
    que.push(s);
    vis[s]=1;dis1[s]=0;
    while(!que.empty()){
        int now=que.front();
        que.pop();
        vis[now]=0;
        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]]=1;
                    que.push(to[i]);
                }
            }
    }
}
bool bfs(){
    for(int i=scr;i<=decc;i++){
        lev[i]=-1;
        cur[i]=head1[i];
    }
    while(!que.empty())    que.pop();
    que.push(scr);
    lev[scr]=0;
    while(!que.empty()){
        int now=que.front();
        que.pop();
        for(int i=head1[now];i;i=net1[i])
            if(lev[to1[i]]==-1&&cap1[i]>0){
                lev[to1[i]]=lev[now]+1;
                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=0,detal;
    for(int & i=cur[now];i;i=net1[i])
        if(lev[to1[i]]==lev[now]+1&&cap1[i]){
            detal=dinic(to1[i],min(flow-rest,cap1[i]));
            if(detal){
                cap1[i]-=detal;
                cap1[i^1]+=detal;
                rest+=detal;
                if(rest==flow)    break;
            }
        }
    if(rest!=flow)    lev[now]=-1;
    return rest;
}
int main(){
    scanf("%d%d",&n,&m);
    scr=1;decc=n;
    for(int i=1;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(1);
    cout<<dis[n]<<endl;
    spfa1(n);
    for(int i=1;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;
}