记录编号 232494 评测结果 AAAAAAAAAA
题目名称 [NOI 2005]聪聪与可可 最终得分 100
用户昵称 Gravatarmikumikumi 是否通过 通过
代码语言 C++ 运行时间 0.343 s
提交时间 2016-03-01 18:11:52 内存使用 11.29 MiB
显示代码纯文本
#include<cstdio>
#include<deque>
#include<iostream>
#include<cstring>
using namespace std;
const int SIZEN=1000;
const double zero=1e-6;
int N,E;
deque<int> s[SIZEN];
int d[SIZEN]={0};
int C,M;
void read()
{
	scanf("%d%d",&N,&E);
	scanf("%d%d",&C,&M);
	int fr,to;
	for(int i=1;i<=E;i++) 
	{
		scanf("%d%d",&fr,&to);
		s[fr].push_back(to);
		s[to].push_back(fr);
		d[fr]++;
		d[to]++;
	}
}
int p[SIZEN][SIZEN]={0};
int dis[SIZEN][SIZEN];
void bfs(int x)
{
	deque<int> Q;
	dis[x][x]=0;
	Q.push_back(x);
	while(!Q.empty())
	{
		int v=Q.front();Q.pop_front();
		int tem=p[x][v];
		for(int i=0;i<s[v].size();i++)
		{
			int to=s[v][i];
			if(dis[x][to]==-1||(dis[x][v]+1==dis[x][to]&&tem<p[x][to]))
			{
				dis[x][to]=dis[x][v]+1;
				if(!tem) p[x][to]=to;
				else p[x][to]=tem;
				Q.push_back(to);
			}
		}
	}
}
double f[SIZEN][SIZEN]={0};
double DP(int x,int y)
{
	if(f[x][y]>zero) return f[x][y];
	if(x==y) return 0;
	if(p[x][y]==y||p[p[x][y]][y]==y) return f[x][y]=1;
	for(int i=0;i<s[y].size();i++)
	{
		int to=s[y][i];
		f[x][y]+=DP(p[p[x][y]][y],to);
	}
	f[x][y]+=DP(p[p[x][y]][y],y);
	f[x][y]/=(d[y]+1.0);
	f[x][y]+=1.0;
	return f[x][y];
}
int main()
{
	freopen("cchkk.in","r",stdin);
	freopen("cchkk.out","w",stdout);
	read();
	memset(dis,-1,sizeof(dis));
	for(int i=1;i<=N;i++) bfs(i);
	double ans=DP(C,M);
	printf("%.3lf\n",ans);
	return 0;
}