记录编号 285400 评测结果 AAAAAAAA
题目名称 [SDOI 2010] 星际竞速 最终得分 100
用户昵称 GravatarTenderRun 是否通过 通过
代码语言 C++ 运行时间 0.443 s
提交时间 2016-07-22 13:25:37 内存使用 7.98 MiB
显示代码纯文本
  1. #include <iostream>
  2. #include <cstring>
  3. #include <cstdio>
  4. #include <queue>
  5. using namespace std;
  6. const int maxn=2010;
  7. const int maxm=500010;
  8. const int INF=1061109567;
  9. int n,m,w[maxn],cnt,fir[maxn],to[maxm],nxt[maxm];
  10. int cap[maxm],val[maxm],path[maxn],vis[maxn],dis[maxn];
  11. queue<int>q;
  12. struct Net_Flow{
  13. Net_Flow(){cnt=1;}
  14. void add(int a,int b,int c,int v){
  15. nxt[++cnt]=fir[a];cap[cnt]=c;
  16. to[cnt]=b;val[cnt]=v;fir[a]=cnt;
  17. }
  18. void addedge(int a,int b,int c,int v){
  19. add(a,b,c,v);add(b,a,0,-v);
  20. }
  21. int Spfa(int S,int T){
  22. memset(dis,63,sizeof(dis));
  23. q.push(S);vis[S]=1;dis[S]=0;
  24. while(!q.empty()){
  25. int x=q.front();q.pop();vis[x]=0;
  26. for(int i=fir[x];i;i=nxt[i])
  27. if(cap[i]&&dis[to[i]]>dis[x]+val[i]){
  28. dis[to[i]]=dis[x]+val[i];
  29. if(!vis[to[i]])q.push(to[i]);
  30. vis[to[i]]=1;path[to[i]]=i;
  31. }
  32. }
  33. return dis[T];
  34. }
  35. int Aug(int S,int T){
  36. int p=T,f=INF;
  37. while(p!=S){
  38. f=min(f,cap[path[p]]);
  39. p=to[path[p]^1];
  40. }p=T;
  41. while(p!=S){
  42. cap[path[p]]-=f;
  43. cap[path[p]^1]+=f;
  44. p=to[path[p]^1];
  45. }
  46. return f;
  47. }
  48. int MCMF(int S,int T){
  49. int v=0,d;
  50. while((d=Spfa(S,T))!=INF)
  51. v+=d*Aug(S,T);
  52. return v;
  53. }
  54. }mcmf;
  55.  
  56. int S,T;
  57. int main(){
  58. freopen("starrace.in","r",stdin);
  59. freopen("starrace.out","w",stdout);
  60. scanf("%d%d",&n,&m);S=0;T=2*n+1;
  61. for(int i=1;i<=n;i++){
  62. scanf("%d",&w[i]);
  63. mcmf.addedge(S,i,1,0);
  64. mcmf.addedge(S,i+n,1,w[i]);
  65. mcmf.addedge(i+n,T,1,0);
  66. }
  67. for(int i=1,a,b,v;i<=m;i++){
  68. scanf("%d%d%d",&a,&b,&v);
  69. if(a>b)swap(a,b);
  70. mcmf.addedge(a,b+n,1,v);
  71. }
  72. printf("%d\n",mcmf.MCMF(S,T));
  73. return 0;
  74. }