bzoj1927 [Sdoi2010]星际竞速 费用流

Description

---

　　10年一度的银河系赛车大赛又要开始了。作为全银河最盛大的活动之一，夺得这个项目的冠军无疑是很多人的
梦想，来自杰森座α星的悠悠也是其中之一。赛车大赛的赛场由N颗行星和M条双向星际航路构成，其中每颗行星都
有一个不同的引力值。大赛要求车手们从一颗与这N颗行星之间没有任何航路的天体出发，访问这N颗行星每颗恰好
一次，首先完成这一目标的人获得胜利。由于赛制非常开放，很多人驾驶着千奇百怪的自制赛车来参赛。这次悠悠
驾驶的赛车名为超能电驴，这是一部凝聚了全银河最尖端科技结晶的梦幻赛车。作为最高科技的产物，超能电驴有
两种移动模式：高速航行模式和能力爆发模式。在高速航行模式下，超能电驴会展开反物质引擎，以数倍于光速的
速度沿星际航路高速航行。在能力爆发模式下，超能电驴脱离时空的束缚，使用超能力进行空间跳跃——在经过一
段时间的定位之后，它能瞬间移动到任意一个行星。天不遂人愿，在比赛的前一天，超能电驴在一场离子风暴中不
幸受损，机能出现了一些障碍：在使用高速航行模式的时候，只能由每个星球飞往引力比它大的星球，否则赛车就
会发生爆炸。尽管心爱的赛车出了问题，但是悠悠仍然坚信自己可以取得胜利。他找到了全银河最聪明的贤者——
你，请你为他安排一条比赛的方案，使得他能够用最少的时间完成比赛。

N≤800，M≤15000。输入数据中的任何数都不会超过106。输入数据保证任意两颗行星

Solution

---

其实我第一眼想到的是dp的，但是发现做不了
考虑费用流。拆点后建图，对于直接到达的连(s,i’)容量为1费用为ai，那些路径就连(i’,j)容量为1费用为w，剩下的就好搞了，小数据spfa能艹过去

思考一下为什么要这么连。把每个点单独出来看，那么就可以直接到达或其他点来。题目要求我们每个点只能访问一次，因此要拆点。

一不小心1a了

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <queue>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))
#define min(x,y) ((x)<(y)?(x):(y))
const int INF=0x3f3f3f3f;
const int N=20005;
const int E=640005;
struct edge{int x,y,w,c,next;}e[E];
std:: queue<int> que;
int ls[N],edCnt=1;
int dis[N],pre[N],a[N];
bool vis[N];
int n;
int read() {
    int x=0,v=1; char ch=getchar();
    for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
    for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
    return x*v;
}
void addEdge(int x,int y,int c,int w) {
    e[++edCnt]=(edge){x,y,w,c,ls[x]}; ls[x]=edCnt;
    e[++edCnt]=(edge){y,x,0,-c,ls[y]}; ls[y]=edCnt;
}
int spfa(int st,int ed) {
    while (!que.empty()) que.pop();
    que.push(st);
    rep(i,0,N-1) dis[i]=INF;
    fill(vis,0);
    dis[st]=0;
    vis[st]=1;
    while (!que.empty()) {
        int now=que.front(); que.pop();
        for (int i=ls[now];i;i=e[i].next) {
            if (e[i].w>0&&dis[now]+e[i].c<dis[e[i].y]) {
                dis[e[i].y]=dis[now]+e[i].c;
                pre[e[i].y]=i;
                if (!vis[e[i].y]) {
                    que.push(e[i].y);
                    vis[e[i].y]=1;
                }
            }
        }
        vis[now]=0;
    }
    return dis[ed]!=INF;
}
int modify(int ed) {
    int mn=INF,cost=0;
    for (int i=ed;pre[i];i=e[pre[i]].x) {
        cost+=e[pre[i]].c;
        mn=min(mn,e[pre[i]].w);
    }
    for (int i=ed;pre[i];i=e[pre[i]].x) {
        e[pre[i]].w-=mn;
        e[pre[i]^1].w+=mn;
    }
    return mn*cost;
}
int mcf(int st,int ed) {
    int ret=0;
    while (spfa(st,ed))
        ret+=modify(ed);
    return ret;
}
int main(void) {
    n=read();
    int m=read();
    rep(i,1,n) {
        a[i]=read();
        addEdge(0,i,0,1);
        addEdge(i+n,n+n+1,0,1);
        addEdge(0,i+n,a[i],1);
    }
    rep(i,1,m) {
        int x=read(),y=read(),w=read();
        if (x<y) addEdge(x,y+n,w,1);
        else addEdge(y,x+n,w,1);
    }
    int prt=mcf(0,n+n+1);
    printf("%d\n", prt);
    return 0;
}