标签:head cur int P4206 NOI2005 聪聪 1005 include dis
预处理 nxt[i,j] 表示i到j最近距离i的下一个位置
最后记忆化搜索就好
#include<cstdio>
#include<queue>
#include<cstring>
#include<iostream>
using namespace std;
int cur,n,m,s,t;
int head[1005],p[1005];
int dis[1005][1005],nxt[1005][1005];
bool vis[1005],visit[1005][1005];
double f[1005][1005];
struct EDGE{
int t,next;
}e[2005];
#define INF 0x3f3f3f3f
void add(int a,int b)
{
cur++;
e[cur].t=b;
e[cur].next=head[a];
head[a]=cur;
}
queue < int > q;
void SPFA(int *dis,int *nxt,int s)
{
dis[s]=0;
q.push(s);
while (!q.empty())
{
int u=q.front();q.pop();
vis[u]=false;
for (int h=head[u];h!=-1;h=e[h].next)
{
int v=e[h].t;
if (dis[u]+1<dis[v])
{
dis[v]=dis[u]+1;
if (!vis[v])
{
vis[v]=true;
q.push(v);
}
}
}
}
}
double DFS(int u,int v)
{
if (visit[u][v]) return f[u][v];
if (u==v) return 0;
int fir=nxt[u][v];
int sec=nxt[fir][v];
if (fir==v||sec==v) return 1;
f[u][v]=1;
for (int h=head[v];h!=-1;h=e[h].next)
{
int w=e[h].t;
f[u][v]+=DFS(sec,w)/(p[v]+1);
}
f[u][v]+=DFS(sec,v)/(p[v]+1);
visit[u][v]=true;
return f[u][v];
}
int main()
{
scanf("%d%d%d%d",&n,&m,&s,&t);
memset(head,-1,sizeof head);
for (int i=1;i<=m;i++)
{
int a,b;
scanf("%d%d",&a,&b);
add(a,b);
add(b,a);
p[a]++;p[b]++;
}
for (int i=1;i<=n;i++)
{
for (int j=1;j<=n;j++)
dis[i][j]=nxt[i][j]=INF;
}
for (int i=1;i<=n;i++)
{
SPFA(dis[i],nxt[i],i);
}
for (int i=1;i<=n;i++)
for (int h=head[i];h!=-1;h=e[h].next)
{
int t=e[h].t;
for (int j=1;j<=n;j++)
if (dis[i][j]-1==dis[t][j])
{
nxt[i][j]=min(nxt[i][j],t);
}
}
printf("%.3lf",DFS(s,t));
return 0;
}
标签:head,cur,int,P4206,NOI2005,聪聪,1005,include,dis 来源: https://www.cnblogs.com/wzxbeliever/p/16178416.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。