记录编号	118682	评测结果	AAAAAAAAAAAAAAAAAAAA
题目名称	[東方S3] 铃仙•优昙华院•稻叶	最终得分	100
用户昵称	HouJikan	是否通过	通过
代码语言	C++	运行时间	2.865 s
提交时间	2014-09-08 21:02:24	内存使用	9.49 MiB

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <list>
#include <vector>
#include <ctime>
#include <iterator>
#include <functional>
#define pritnf printf
#define scafn scanf
#define For(i,j,k) for(int i=(j);i<=(k);(i)++)
using namespace std;
typedef long long LL;
typedef unsigned int Uint; 
const int INF=0x7fffffff;
//==============struct declaration==============
struct node
{
  int now,from;
  int time;
  node (int now=0,int from=0,int time=0):now(now),from(from),time(time){}
};
//==============var declaration=================
const int MAXN=55;
int n,m,t;
bool inq[510][MAXN][MAXN];//t秒时在i节点，上一秒在j节点
bool path[MAXN][MAXN];
int  out[MAXN]; 
double poss[510][MAXN][MAXN]; 
//==============function declaration============

//==============main code=======================
int main()
{  
  string FileName="reisen";//程序名 
  string FloderName="COGS";//文件夹名 
  freopen((FileName+".in").c_str(),"r",stdin);
  freopen((FileName+".out").c_str(),"w",stdout);
#ifdef DEBUG  
  system(("cp C:\\Users\\Administrator\\Desktop\\"+FloderName+"\\standard.cpp C:\\Users\\Administrator\\Desktop\\"+FloderName+"\\submit.txt").c_str());
  clock_t Start_Time=clock();
#endif    
  scanf("%d%d%d",&n,&m,&t);
  memset(out,false,sizeof(out));
  memset(path,false,sizeof(path));
  For(i,1,m)
  {
    int s,e;
    scanf("%d%d",&s,&e);
    path[s][e]=true;
    out[s]++;
  }
  For(i,1,n)
  {
    out[i]++;
    path[i][i]=true;
  } 
  memset(poss,0,sizeof(poss));
  memset(inq,false,sizeof(inq));
  queue <node> Q;
  poss[0][1][0]=1;
  Q.push(node(1,0,0));
  while (!Q.empty())
  {
    node x=Q.front();Q.pop();
    if (x.time==t) break;
    inq[x.time][x.now][x.from]=false;
    double next;
    if (!path[x.now][x.from])
      next=1.0/out[x.now];
    else 
      next=1.0/(out[x.now]-1.0);
    For(i,1,n)
      if (i!=x.from&&path[x.now][i])
      {
        if (i!=x.now)
        {
          poss[x.time+1][i][x.now]+=poss[x.time][x.now][x.from]*next;
          if (!inq[x.time+1][i][x.now])
          {
            inq[x.time+1][i][x.now]=true;
            Q.push(node(i,x.now,x.time+1));
          }
        }
        else
        {
          poss[x.time+1][x.now][x.from]+=poss[x.time][x.now][x.from]*next;
          if (!inq[x.time+1][x.now][x.from])
          {
            inq[x.time+1][x.now][x.from]=true;
            Q.push(node(i,x.from,x.time+1));
          }
        } 
      }
  }
  For(i,1,n)
  {
    double sum=0;
    For(j,0,n)
      sum+=poss[t][i][j];
    printf("%.3lf\n",sum*100);
  }
#ifdef DEBUG  
  clock_t End_Time=clock();
  cout<<endl<<endl<<"Time Used: "<<End_Time-Start_Time<<endl;
#endif    
  return 0;
}
//================fuction code====================