请柬(invite)

【问题描述】

在电视时代,没有多少人观看戏剧表演。Malidinesia古董喜剧演员意识到这一事实，他们想宣传剧院,尤其是古色古香的喜剧片。他们已经打印请帖和所有必要的信息和计划。许多学生被雇来分发这些请柬。每个学生志愿者被指定一个确切的公共汽车站，他或她将留在那里一整天,邀请人们参与。

　　这里的公交系统是非常特殊的:所有的线路都是单向的，连接两个站点。公共汽车离开起始点，到达目的地之后又空车返回起始点。　

学生每天早上从总部出发，乘公交车到一个预定的站点邀请乘客。每个站点都被安排了一名学生。在一天结束的时候,所有的学生都回到总部。现在需要知道的是，学生所需的公交费用的总和最小是多少。

【输入文件】

第1行有两个整数n、m(1<=n,m<=1000000)，n是站点的个数，m是线路的个数。

然后有m行，每行描述一个线路，包括3个整数，起始点，目的地和价格。

总部在第1个站点，价钱都是整数，且小于1000000000。

【输出文件】

输出一行，表示最小费用。

【样例输入】

4 6

1 2 10

2 1 60

1 3 20

3 4 10

2 4 5

4 1 50

【样例输出】

210

【样例解释】

学生各自从总部被派遣到2,3,4站点，然后又回到总部

1-2-4-1:10+5+50=65

1-3-4-1:20+10+50=80

1-2-4-1:10+5+50=65

65+80+65=210

【注意】

此题数据规模较大，需要使用较为高效的算法，此题不设小规模数据分数。

建立双向图然后两次单源最短路，用dijkstra+heap练练手

program zhen;type   sonlink=^node;   node=record      e:longint;      v:int64;      next:sonlink;   end;var   son:array[0..1000000,1..2]of sonlink;   dist:array[0..1000000]of int64;   heap,hpos:array[0..1000000]of longint;   i,j,n,m,len:longint;   ans:int64;procedure init_data;var   i,x,y,v:longint;   p:sonlink;begin   readln(n,m);   for i:=1 to m do     begin       readln(x,y,v);       new(p);         p^.v:=v;         p^.e:=y;         p^.next:=son[x,1];son[x,1]:=p;       new(p);         p^.v:=v;         p^.e:=x;         p^.next:=son[y,2];son[y,2]:=p;     end;end;procedure swap(var x,y:longint);var   t:longint;begin   t:=x;x:=y;y:=t;end;procedure sifup(s:longint);var   p:longint;begin   while s>1 do     begin       p:=s shr 1;       if dist[heap[p]]>dist[heap[s]] then         begin           hpos[heap[p]]:=s;hpos[heap[s]]:=p;           swap(heap[p],heap[s]);         end       else         break;       s:=p;     end;end;procedure sifdown(p:longint);var   s:longint;begin   while p shl 1<=len do     begin       s:=p shl 1;       if (s<len) and (dist[heap[s+1]]<dist[heap[s]]) then inc(s);       if dist[heap[p]]>dist[heap[s]] then         begin           hpos[heap[p]]:=s;hpos[heap[s]]:=p;           swap(heap[p],heap[s]);         end       else         break;       p:=s;     end;end;procedure extract_min;begin   heap[1]:=heap[len];   hpos[heap[1]]:=1;   dec(len);   sifdown(1);end;procedure dijkstra(x:longint);var   i,j,k:longint;   t:int64;   p:sonlink;begin   dist[1]:=0;   for i:=2 to n do dist[i]:=1000000000000;   len:=n;   for i:=1 to n do     begin       hpos[i]:=i;       heap[i]:=i;     end;   for i:=1 to n do     begin       k:=heap[1];       extract_min;       p:=son[k,x];       while p<>nil do         begin           t:=dist[k]+p^.v;           if t<dist[p^.e] then             begin               dist[p^.e]:=t;               sifup(hpos[p^.e]);             end;           p:=p^.next;         end;     end;   for i:=2 to n do     ans:=ans+dist[i];end;begin   assign(input,'invite.in');reset(input);   assign(output,'invite.out');rewrite(output);   init_data;   ans:=0;   dijkstra(1);   dijkstra(2);   writeln(ans);   close(input);close(output);end.