【BZOJ1579】Revamping Trails（分层图，最短路，堆）

题意：每天,农夫John需要经过一些道路去检查牛棚N里面的牛. 农场上有M(1<=M<=50,000)条双向泥土道路,编号为1..M.

道路i连接牛棚P1_i和P2_i (1 <= P1_i <= N; 1 <= P2_i<= N). John需要T_i (1 <= T_i <= 1,000,000)时间单位用道路i从P1_i走到P2_i或者从P2_i 走到P１_i

他想更新一些路经来减少每天花在路上的时间.具体地说,他想更新K (1 <= K <= 20)条路经，将它们所须时间减为０．

帮助FJ选择哪些路经需要更新使得从１到N的时间尽量少.

N<=10000

思路：裂点建分层图后跑最短路

C++ dijksta+堆板子

 #include<cstdio>
 #include<cstring>
 #include<string>
 #include<cmath>
 #include<iostream>
 #include<algorithm>
 #include<map>
 #include<set>
 #include<queue>
 #include<vector>
 using namespace std;
 typedef long long ll;
 typedef unsigned int uint;
 typedef unsigned long long ull;
 typedef pair<int,int> PII;
 typedef vector<int> VI;
 #define fi first
 #define se second
 #define MP make_pair
 #define N   11000
 #define M   51000
 #define eps 1e-8
 #define pi  acos(-1)
 #define oo  1e9
 priority_queue<pair<int,int> > q;
 int num[N][],head[N*],vet[M*],len[M*],nxt[M*],
     dis[N*],vis[N*],tot;

 int add(int a,int b,int c)
 {
     nxt[++tot]=head[a];
     vet[tot]=b;
     len[tot]=c;
     head[a]=tot;
 }

 int main()
 {
     int n,m,K;
     scanf("%d%d%d",&n,&m,&K);
     int s=;
     for(int i=;i<=n;i++)
      for(int j=;j<=K;j++) num[i][j]=++s;
     for(int i=;i<=m;i++)
     {
         int x,y,z;
         scanf("%d%d%d",&x,&y,&z);
         for(int j=;j<K;j++)
         {
             add(num[x][j],num[y][j+],);
             add(num[y][j],num[x][j+],);
         }
         for(int j=;j<=K;j++)
         {
             add(num[x][j],num[y][j],z);
             add(num[y][j],num[x][j],z);
         }
     }
     memset(vis,,sizeof(vis));
     memset(dis,0x3f,sizeof(dis));
     q.push(MP(,num[][])); dis[num[][]]=;
     while(!q.empty())
     {
         int u=q.top().se;
         q.pop();
         if(vis[u]) continue;
         vis[u]=;
         int e=head[u];
         while(e)
         {
             int v=vet[e];
             if(dis[u]+len[e]<dis[v])
             {
                 dis[v]=dis[u]+len[e];
                 q.push(MP(-dis[v],v));
             }
             e=nxt[e];
         }
     }
     int ans=oo;
     for(int i=;i<=K;i++) ans=min(ans,dis[num[n][i]]);
     printf("%d\n",ans);
     return ;
 }