洛谷 1144 最短路计数

题目描述
    给出一个N个顶点M条边的无向无权图，顶点编号为1～N。问从顶点1开始，到其他每个点的最短路有几条。

输入格式：
    输入第一行包含2个正整数N，M，为图的顶点数与边数。
    接下来M行，每行两个正整数x, y，表示有一条顶点x连向顶点y的边，请注意可能有自环与重边。

输出格式：
    输出包括N行，每行一个非负整数，第i行输出从顶点1到顶点i有多少条不同的最短路，由于答案有可能会很大，你只需要输出mod 100003后的结果即可。如果无法到达顶点i则输出0。

输入样例#1：
5 7
1 2
1 3
2 4
3 4
2 3
4 5
4 5
输出样例#1：
1
1
1
2

4

解题思路：

    SPFA

    分三类： 
    没遍历过 根据dis 求出ans

    到过的，u的深度+1是v深度 更新ans；
    没到过这个点的或dis[u]+1>dis[v] 不做处理

    

代码：

var
  num,n,m,i,j,l,r,u,v,e:longint;
  p,head,f,next,dis,ans:array[1..2000000]of longint;
  line:array[1..10000000]of longint;


procedure add(x,y:longint);
begin
  inc(num);
  f[num]:=y;
  next[num]:=head[x];
  head[x]:=num;
end;


begin
  readln(n,m);
  dis[1]:=1;p[1]:=1;ans[1]:=1;
  line[1]:=1;
  for i:=1 to m do
    begin
      readln(l,r);
      add(l,r);
      add(r,l);
    end;
  l:=1;r:=1;
  while l<=r do
    begin
      u:=line[l];
      inc(l);
      e:=head[u];
      while e<>0 do
        begin
          v:=f[e];
          if p[v]=0 then
            begin
              dis[v]:=dis[u]+1;
              ans[v]:=ans[u];
              p[v]:=1;
              r:=r+1;
              line[r]:=v;
            end
          else if dis[v]=dis[u]+1 then ans[v]:=(ans[v]+ans[u])mod 100003;
          e:=next[e];
        end;
    end;
  for i:=1 to n do
    writeln(ans[i]);
end.