1927: [Sdoi2010]星际竞速
Time Limit: 20 Sec Memory Limit: 259 MBSubmit: 2051 Solved: 1263
[Submit][Status][Discuss]
Description
10年一度的银河系赛车大赛又要开始了。作为全银河最盛大的活动之一,夺得这个项目的冠军无疑是很多人的
梦想,来自杰森座α星的悠悠也是其中之一。赛车大赛的赛场由N颗行星和M条双向星际航路构成,其中每颗行星都
有一个不同的引力值。大赛要求车手们从一颗与这N颗行星之间没有任何航路的天体出发,访问这N颗行星每颗恰好
一次,首先完成这一目标的人获得胜利。由于赛制非常开放,很多人驾驶着千奇百怪的自制赛车来参赛。这次悠悠
驾驶的赛车名为超能电驴,这是一部凝聚了全银河最尖端科技结晶的梦幻赛车。作为最高科技的产物,超能电驴有
两种移动模式:高速航行模式和能力爆发模式。在高速航行模式下,超能电驴会展开反物质引擎,以数倍于光速的
速度沿星际航路高速航行。在能力爆发模式下,超能电驴脱离时空的束缚,使用超能力进行空间跳跃——在经过一
段时间的定位之后,它能瞬间移动到任意一个行星。天不遂人愿,在比赛的前一天,超能电驴在一场离子风暴中不
幸受损,机能出现了一些障碍:在使用高速航行模式的时候,只能由每个星球飞往引力比它大的星球,否则赛车就
会发生爆炸。尽管心爱的赛车出了问题,但是悠悠仍然坚信自己可以取得胜利。他找到了全银河最聪明的贤者——
你,请你为他安排一条比赛的方案,使得他能够用最少的时间完成比赛。
Input
第一行是两个正整数N,M。第二行N个数A1~AN,其中Ai表示使用能力爆发模式到达行星i所需的定位时间。接下
来M行,每行3个正整数ui,vi,wi,表示在编号为ui和vi的行星之间存在一条需要航行wi时间的星际航路。输入数据
已经按引力值排序,也就是编号小的行星引力值一定小,且不会有两颗行星引力值相同。
Output
仅包含一个正整数,表示完成比赛所需的最少时间。
Sample Input
1 100 100
2 1 10
1 3 1
2 3 1
Sample Output
HINT
说明:先使用能力爆发模式到行星1,花费时间1。然后切换到高速航行模式,航行到行星2,花费时间10。之
后继续航行到行星3完成比赛,花费时间1。虽然看起来从行星1到行星3再到行星2更优,但我们却不能那样做,因
为那会导致超能电驴爆炸。N≤800,M≤15000。输入数据中的任何数都不会超过106。输入数据保证任意两颗行星
之间至多存在一条航道,且不会存在某颗行星到自己的航道。
Source
分析:
就是说我们要寻找一条道路使得经过每个点一次且仅一次,但是网络流的建图无法直接反映这种道路,所以我们考虑这种道路可以转化成什么...
如果存在这种道路,那么一定是每个点入度出度都等于1,所以我们要寻找若干条道路使得每个点入度出度都为1...
很自然的我们会想到拆点,每个点拆成一个入点一个出点,S向出点连边,容量为1费用为0,入点向T连边容量为1费用为0...然后对于星际轨道假设x<y我们从x的出点向y的入点连边,容量为1费用为时间,然后从S向每个入点连边,容量为1费用为瞬移到这个星球的时间...
这样汇点的边的容量限制以及最大流保证了每个入点有且仅有一条有流量的入边...这条边要么从某个点出来要么从源点瞬移过来...
代码:
1 #include<algorithm> 2 #include<iostream> 3 #include<cstring> 4 #include<cstdio> 5 #include<queue> 6 //by NeighThorn 7 #define inf 0x3f3f3f3f 8 using namespace std; 9 10 const int maxn=2000+5,maxm=40000+5; 11 12 int n,m,S,T,cnt,w[maxm],hd[maxn],fl[maxm],to[maxm],nxt[maxm],dis[maxn],Min[maxn],vis[maxn],from[maxn]; 13 14 inline bool spfa(void){ 15 memset(dis,inf,sizeof(dis)); 16 memset(Min,inf,sizeof(Min)); 17 queue<int> q;q.push(S),vis[S]=1,dis[S]=0; 18 while(!q.empty()){ 19 int top=q.front();q.pop(),vis[top]=0; 20 for(int i=hd[top];i!=-1;i=nxt[i]) 21 if(dis[to[i]]>dis[top]+w[i]&&fl[i]){ 22 from[to[i]]=i; 23 dis[to[i]]=dis[top]+w[i]; 24 Min[to[i]]=min(Min[top],fl[i]); 25 if(!vis[to[i]]) 26 vis[to[i]]=1,q.push(to[i]); 27 } 28 } 29 return dis[T]!=inf; 30 } 31 32 inline int find(void){ 33 for(int i=T;i!=S;i=to[from[i]^1]) 34 fl[from[i]]-=Min[T],fl[from[i]^1]+=Min[T]; 35 return dis[T]*Min[T]; 36 } 37 38 inline int mcmf(void){ 39 int res=0; 40 while(spfa()) 41 res+=find(); 42 return res; 43 } 44 45 inline void add(int l,int s,int x,int y){ 46 w[cnt]=l;fl[cnt]=s;to[cnt]=y;nxt[cnt]=hd[x];hd[x]=cnt++; 47 w[cnt]=-l;fl[cnt]=0;to[cnt]=x;nxt[cnt]=hd[y];hd[y]=cnt++; 48 } 49 50 signed main(void){ 51 // freopen("in.txt","r",stdin); 52 memset(hd,-1,sizeof(hd)); 53 scanf("%d%d",&n,&m);S=0,T=n*2+1; 54 for(int i=1,x;i<=n;i++) 55 scanf("%d",&x),add(x,1,S,i),add(0,1,S,i+n),add(0,1,i,T); 56 for(int i=1,s,x,y;i<=m;i++){ 57 scanf("%d%d%d",&x,&y,&s); 58 if(x>y) 59 swap(x,y); 60 add(s,1,x+n,y); 61 } 62 printf("%d\n",mcmf()); 63 return 0; 64 }
By NeighThorn