BZOJ1266: [AHOI2006]上学路线route

1266: [AHOI2006]上学路线route

Time Limit: 3 Sec  Memory Limit: 162 MB
Submit: 1273  Solved: 435
[Submit][Status]

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

题解：

第二问的做法比较炫酷。

把所有可以在最短路上的边拎出来，然后做一遍最小割，就OK了，这使得s到t没有一条完全有最短路组成的边，显然增大了s到t的最短路。

代码：不是WA就是T，我已经无语了。。。

#include<cstdio>
 
#include<cstdlib>
 
#include<cmath>
 
#include<cstring>
 
#include<algorithm>
 
#include<iostream>
 
#include<vector>
 
#include<map>
 
#include<set>
 
#include<queue>
 
#include<string>
 
#define inf 1000000000
 
#define maxn 100000
 
#define maxm 500000
 
#define eps 1e-10
 
#define ll long long
 
#define pa pair<int,int>
 
#define for0(i,n) for(int i=0;i<=(n);i++)
 
#define for1(i,n) for(int i=1;i<=(n);i++)
 
#define for2(i,x,y) for(int i=(x);i<=(y);i++)
 
#define for3(i,x,y) for(int i=(x);i>=(y);i--)
 
using namespace std;
 
inline int read()
 
{
 
    int x=0,f=1;char ch=getchar();
 
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
 
    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
 
    return x*f;
 
}
int  n,m,s,t,maxflow,tot=1,head[maxn],cur[maxn],h[maxn],q[maxm];
int x[maxm],y[maxm],z[maxm],w[maxm],d[2][maxn];
 
struct edge{int go,next,v;}e[2*maxm];
bool v[maxn];
 
void ins(int x,int y,int z){e[++tot].go=y;e[tot].v=z;e[tot].next=head[x];head[x]=tot;}
 
void insert(int x,int y,int z){ins(x,y,z);ins(y,x,0);}
 
bool bfs()
 
{
 
    for(int i=s;i<=t;i++)h[i]=-1;
 
    int l=0,r=1;q[1]=s;h[s]=0;
 
    while(l<r)
 
    {
 
        int x=q[++l];
 
        for(int i=head[x];i;i=e[i].next)
 
         if(e[i].v&&h[e[i].go]==-1)
 
         {
 
            h[e[i].go]=h[x]+1;q[++r]=e[i].go;
 
         }
 
    }
 
    return h[t]!=-1;
 
}
 
int dfs(int x,int f)
 
{
 
    if(x==t) return f;
 
    int tmp,used=0;
 
    for(int i=head[x];i;i=e[i].next)
 
     if(e[i].v&&h[e[i].go]==h[x]+1)
 
    {
 
        tmp=dfs(e[i].go,min(e[i].v,f-used));
 
        e[i].v-=tmp;if(e[i].v)cur[x]=i;
 
        e[i^1].v+=tmp;used+=tmp;
 
        if(used==f)return f;       
 
    }
 
    if(!used) h[x]=-1;
 
    return used;
 
}
 
void dinic()
 
{
 
    maxflow=0;
 
    while(bfs())
 
    {
 
        for (int i=s;i<=t;i++)cur[i]=head[i];maxflow+=dfs(s,inf);
 
    }
 
}
 
void spfa(int k)
 
{
 
    for(int i=1;i<=n;++i) d[k][i]=inf;
 
    memset(v,0,sizeof(v));
    s=k==0?1:n;
 
    int l=0,r=1,x;q[1]=s;d[k][s]=0;
 
    while(l<r)
 
    {
 
        x=q[++l];v[x]=0;
 
        for(int i=head[x],y;i;i=e[i].next)
 
         if(d[k][x]+e[i].v<d[k][y=e[i].go])
 
         {
 
            d[k][y]=d[k][x]+e[i].v;
            
 
            if(!v[i]){v[i]=1;q[++r]=i;}
 
         }
 
    }
 
}
inline void ins2(int x,int y,int z)
{
   e[++tot].go=y;e[tot].v=z;e[tot].next=head[x];head[x]=tot;
   e[++tot].go=x;e[tot].v=z;e[tot].next=head[y];head[y]=tot;
}
int main()
 
{
    freopen("input.txt","r",stdin);
    freopen("output.txt","w",stdout);
    n=read();m=read();
    for1(i,m)
     {
         x[i]=read();y[i]=read();z[i]=read();w[i]=read();
         ins2(x[i],y[i],z[i]);
     }
    spfa(0);spfa(1);
    tot=1;
    memset(head,0,sizeof(head));
    for1(i,m)
    {
        if(d[0][x[i]]+z[i]+d[1][y[i]]==d[0][n])insert(x[i],y[i],w[i]);
        if(d[0][y[i]]+z[i]+d[1][x[i]]==d[0][n])insert(y[i],x[i],w[i]);
    }
    s=0;t=n+1;
    insert(s,1,inf);insert(n,t,inf);
    dinic();
    printf("%d\n",d[0][n]);
    printf("%d\n",maxflow);
 
    return 0;
 
}

无奈贴了标程

#include<cstdio>
#include<cstring>
#include<algorithm>
#define INF 2100000000
#define N 1005
#define M 250005
using namespace std;
struct arr{int go,c,next,s;}a[M];
int map[N][N],dis[2][N],f[N],q[M*2],end[N];
bool flag[N];
int x,y,z,c,i,n,m,cnt,ans;
void add(int u,int v,int s,int w)
{
  a[++cnt].go=v;a[cnt].c=w;a[cnt].s=s;a[cnt].next=end[u];
  end[u]=cnt;
}
void SPFA(int o,int sta)
{
  int h=0,t=1;q[1]=sta;
  memset(flag,0,sizeof(flag));
  dis[o][sta]=0;flag[sta]=true;
  while (h<t)
  {
    int now=q[++h];
    for (int i=end[now];i;i=a[i].next)
    {
      int go=a[i].go;
      if (dis[o][now]+a[i].s<dis[o][go])
      {
        dis[o][go]=dis[o][now]+a[i].s;
        if (!flag[go]) flag[go]=true,q[++t]=go;
      }
    }
    flag[now]=false;
  }
} 
void init()
{
  for (int i=1;i<=n;i++)
    for (int j=end[i];j;j=a[j].next)
      if (dis[0][i]+a[j].s+dis[1][a[j].go]==dis[0][n])
        map[i][a[j].go]+=a[j].c;
}
bool bfs()
{
  memset(f,-1,sizeof(f));
  int h=0,t=1;q[1]=1;f[1]=1;
  while (h<t)
  {
    int now=q[++h];if (now==n) return 1;
    for (int i=1;i<=n;i++)
      if (map[now][i]&&f[i]==-1)
      {
        f[i]=f[now]+1;
        q[++t]=i;
      }
  }
  return 0;
}
int dinic(int sta,int sum)
{
  if (sta==n) return sum;int os=sum;
  for (int i=1;(i<=n)&&os;i++)
    if (map[sta][i]&&f[i]==f[sta]+1)
    {
      int Min=dinic(i,min(map[sta][i],os));
      map[sta][i]-=Min;map[i][sta]+=Min;os-=Min;
    }
  if (os==sum) f[sta]=-1;return sum-os;
}
int main()
{
  scanf("%d%d",&n,&m);
  for (i=1;i<=m;i++)
    scanf("%d%d%d%d",&x,&y,&z,&c),add(x,y,z,c),add(y,x,z,c);
  memset(dis,60,sizeof(dis));SPFA(0,1);SPFA(1,n);
  ans=0;init();
  while (bfs())
    ans+=dinic(1,INF);
  printf("%d\n%d",dis[0][n],ans);
  return 0;
}

挖坑，有空来看。。。