BZOJ1877（洛谷P2153）[SDOI2009]晨跑

费用流

BZOJ题目传送门

洛谷题目传送门

考虑这么建图：
①因为每个点只能访问一次，因此考虑进行拆点，把每个点拆成入点x与出点y（1和n除外），在x与y之间建一条容量为1，距离为0的边。
②对于从1开始的点，建一条1–>i.x容量为1的边。
③对于指向n的点，建一条i.y–>n容量为1的边。
④对于其他点u,v,建一条u.y–>v.x容量为1的边。

最后跑一遍费用流即可，最大流的答案即为天数，费用即为路程和。

代码：

#include<cstdio>#include<cstring>#include<algorithm>#define MAXN 400#define MAXM 40000using namespace std;struct edge{    int next,to,dis,v,flow;};struct father{    int x,e;};int n,m,k;int h[MAXN+5],dis[MAXN+5],rem[MAXN+5],que[2*MAXM+5];bool f[MAXN+5];edge ed[2*MAXM+5];father fa[MAXN+5];int spfa(){    memset(f,false,sizeof(f));    memset(dis,0x3f3f3f3f,sizeof(dis));    dis[1]=0; que[1]=1; rem[1]=0x7fffffff;    int r=0,w=1;    while (r<w){        int x=que[++r];        f[x]=false;        for (int i=h[x];~i;i=ed[i].next)            if (ed[i].v>ed[i].flow&&dis[ed[i].to]>dis[x]+ed[i].dis){                int v=ed[i].to;                dis[v]=dis[x]+ed[i].dis;                fa[v].x=x; fa[v].e=i;                rem[v]=min(rem[x],ed[i].v-ed[i].flow);                if (!f[v]){ f[v]=true; que[++w]=v;  }            }    }    if (dis[n]==0x3f3f3f3f) return 0;    return rem[n];}void change(int sum){    int now=n;    while (now!=1){        int e=fa[now].e;        ed[e].flow+=sum; ed[e^1].flow-=sum; now=fa[now].x;    }}void mcmf(){    int ans=0,cst=0;    while (1){        int sum=spfa();        if (!sum) { printf("%d %d\n",ans,cst); return;  }        ans++; cst+=sum*dis[n]; change(sum);    }}void addedge(int x,int y,int z,int w){    ed[k].next=h[x]; ed[k].to=y; ed[k].v=z; ed[k].dis=w; h[x]=k++;    ed[k].next=h[y]; ed[k].to=x; ed[k].v=0; ed[k].dis=-w; h[y]=k++;}int main(){    scanf("%d%d",&n,&m);    memset(h,-1,sizeof(h));    for (int i=1;i<=m;i++){        int u,v,d;        scanf("%d%d%d",&u,&v,&d);        if (u!=1&&v!=n) addedge(u+n,v,1,d);        if (u==1&&v!=n) addedge(u,v,1,d);        if (u!=1&&v==n) addedge(u+n,v,1,d);        if (u==1&&v==n) addedge(u,v,1,d);    }    for (int i=2;i<n;i++)        addedge(i,i+n,1,0);    mcmf();    return 0;}