雀巢咖啡杯~.~(二) 难得感觉到了最短路的神奇

   

       同样大萎的.......第二试，打裸的程序比一试更令人#.#

     

     第一题题意是给出n个集合，第i个集合的元素自然数y满足y= s[i]+k*d[i] ，k属于自然数集， 且y<= e[i]，满足每个集合的个数的总数不超过10^8, n 不超过200000， e不超过10^9 个，保证有唯一一个数n个集合中出现奇数次，要求找出这个数；

     算法是好想的， 奇数只有一个，二分要求的数，o（m）求出在区间【1，m】的数的和，并判奇偶就行了。

     当然，亦可以ws的把集合中的数xor 起来.............

    第三题不解释了............

 ============================以上的都是not important 的===================================

       给定一个V<=1000 E<=v^2 的无向图，要求你在图中选取若V-1条边，构成一棵树，使得在树中每个节点到1节点的最短距离，等于该节点在原图中到1节点的最短距离相同，并问构造的方案数。

       记得曾经一道最小生成树计数，让me 感受到了kruscal 最小生成树的美妙性质。

       其实，最短路和最小生成树..............意会吧，大家都懂的。

       最短路的性质也很美妙，首先，把dijstra 决策的边弄出来，是一个树(考场上一直在纠结树怎么构造....@.@)，

而且，已经标记永久化的节点，所有的东西都确定了，其他的节点不管怎么乱搞都不会对它有影响了。

       所以，处理一下dijstra，在标记永久化i点的时候/顺便求一下/有几条边连向i的边/可以在到达i的最短路上（绕口），在用乘法原理统计之...............

        

      继续陶醉ing..............

  贴个代码吧：

program lmd;
const mo=1 shl 31-1;
var
   c,next,l,sum:array[0..1000000]of longint;
   dist,way:array[0..3000]of longint;
   yes:array[0..3000]of boolean;
   k,i,j,bj,x,y,z,n,m,top:longint;
   ans:int64;
procedure inf;
begin
    assign(input,'castle.in');
    assign(output,'castle.out');
    reset(input); rewrite(output);
end;
procedure ouf;
begin
     close(Input); close(output);
end;
procedure link(x,y,z:longint);
begin
   inc(top);
   next[top]:=l[x];
   l[x]:=top;
   sum[top]:=y;
   c[top]:=z;
end;
procedure init;
begin
     read(n,m);
     for i:=1 to m do
        begin
           read(x,y,z);
           link (x,y,z);
           link (y,x,z);
        end;
     fillchar(dist,sizeof(dist),127);
     fillchar(yes,sizeof(yes),true);
     dist[1]:=0;
     way[1]:=1;
end;
procedure dijstra;
begin
    for i:=1 to n do
      begin
         bj:=0;
         for j:=1 to n do
          if yes[j] and (dist[j]<dist[bj]) then
            bj:=j;
         yes[bj]:=false;
         if bj<>0then
           begin
              k:=l[bj];
              while k<>0 do
                begin
                  if dist[sum[k]]>dist[bj]+c[k] then
                    begin
                      dist[sum[k]]:=dist[bj]+c[k];
                      way[sum[k]]:=1;
                    end
                  else if  dist[sum[k]]=dist[bj]+c[k] then
                     inc(way[sum[k]],1);
                  k:=next[k];
                end;
           end;
      end;
end;
procedure print;
begin
   ans:=1;
   for i:=1 to n do
     ans:=(ans*way[i]) mod mo;
   write(ans);
end;
begin
   inf; init;
   dijstra;
   print;
   ouf;
end.