zoukankan      html  css  js  c++  java
  • TYVJ 1741 [SDOI2010]星际竞速 最小费用流

    这个题联系的时候没有做出来。最后还是lyd神犇教的。。。我还是太弱了。。

    每个点拆成一个入点(右边一排)和一个出点(左边一排),源点S向每个出点连容量1费用0的边,每个入点向汇点连容量1费用0的边,如果有高速航道(x,y),注意x<y,那么从x的出点到y的入点连容量1费用为读入权值的边。这样就能够处理所有高速航道了……

    空间跳跃处理:从S向所有入点连容量1费用为定位费用的边就行了。

    流过一个点对应的出点不代表经过了这个点,只有流过了这个点对应的的入点才算真正流过了这个点!

    View Code
     1 #include <iostream>
     2 #include <algorithm>
     3 #include <cstring>
     4 #include <cstdio>
     5 #include <cstdlib>
     6 
     7 #define N 10000
     8 #define M 1000000
     9 
    10 using namespace std;
    11 
    12 int head[N],next[M],to[M],pr[M],len[M];
    13 int q[M*5],dis[N],pre[N];
    14 bool vis[N];
    15 int n,m,cnt,S,T;
    16 
    17 inline void add(int u,int v,int r,int w)
    18 {
    19     to[cnt]=v; len[cnt]=r; pr[cnt]=w; next[cnt]=head[u]; head[u]=cnt++;
    20     to[cnt]=u; len[cnt]=0; pr[cnt]=-w; next[cnt]=head[v]; head[v]=cnt++;
    21 }
    22 
    23 inline void read()
    24 {
    25     memset(head,-1,sizeof head); cnt=0;
    26     scanf("%d%d",&n,&m);
    27     S=0; T=2*n+1;
    28     for(int i=1,a;i<=n;i++)
    29     {
    30         scanf("%d",&a);
    31         add(S,i+n,1,a);
    32     }
    33     for(int i=1,a,b,c;i<=m;i++)
    34     {
    35         scanf("%d%d%d",&a,&b,&c);
    36         if(a>b) swap(a,b);
    37         add(a,b+n,1,c);
    38     }
    39     for(int i=1;i<=n;i++) add(S,i,1,0),add(i+n,T,1,0);
    40 }
    41 
    42 inline bool spfa()
    43 {
    44     memset(dis,0x3f,sizeof dis);
    45     memset(pre,-1,sizeof pre);
    46     int h=1,t=2,sta;
    47     q[1]=S; vis[S]=true; dis[S]=0;
    48     while(h<t)
    49     {
    50         sta=q[h++]; vis[sta]=false;
    51         for(int i=head[sta];~i;i=next[i])
    52             if(len[i]&&dis[to[i]]>dis[sta]+pr[i])
    53             {
    54                 dis[to[i]]=dis[sta]+pr[i];
    55                 pre[to[i]]=i;
    56                 if(!vis[to[i]]) vis[to[i]]=true,q[t++]=to[i];
    57             }
    58     }
    59     return pre[T]!=-1;
    60 }
    61 
    62 inline void updata()
    63 {
    64     for(int i=pre[T];~i;i=pre[to[i^1]])
    65     {
    66         len[i]-=1; len[i^1]+=1;
    67     }
    68 }
    69 
    70 inline void go()
    71 {
    72     int ans=0;
    73     while(spfa()) ans+=dis[T],updata();
    74     printf("%d\n",ans);
    75 }
    76 
    77 int main()
    78 {
    79     read();
    80     go();
    81     return 0;
    82 }
    没有人能阻止我前进的步伐,除了我自己!
  • 相关阅读:
    小议sql查询返回xml数据之应用【转载】 sansan
    JScript中Date.getTime转.Net中的DateTime sansan
    iFrame 跨域高度自适应问题解决 sansan
    使用第三方应用(天气预报、Google地图之类)不影响原来页面的加载 sansan
    【转载】今天心情非常好,再发一组 Linq、 集合、数组、Lambda、QuerySyntax 的文章 sansan
    高性能ASP.NET站点构建之简单的优化措施 sansan
    Linq to SQL sansan
    高并发量网站解决方案
    JAVA开源项目[转]
    淘宝下单高并发解决方案
  • 原文地址:https://www.cnblogs.com/proverbs/p/2848043.html
Copyright © 2011-2022 走看看