NOI2005聪明与可可

来自蒟蒻XXJ的学习记录

可以看《浅析竞赛中一类数学期望问题的解决方法》汤可因

#include<bits/stdc++.h>
#define mem(i,j) memset(i,j,sizeof(i))
#define GO(i,here) for(int i=head[here];i!=-1;i=nex[i])
using namespace std;
const int MAXN=1024;
const int MAXE=2048;
inline int in(){
    int a(0);char c=getchar();
    while(c<'0'||c>'9') c=getchar();
    while(c>='0'&&c<='9') a=(a<<1)+(a<<3)+c-'0',c=getchar();
    return a;
}
int n,e,c,m;
int head[MAXN],to[MAXE],nex[MAXE],top1;
int p[MAXN][MAXN],vis[MAXN];
double f[MAXN][MAXN],t[MAXN];
queue<int> q;
void bfs(int s){
    q.push(s);mem(vis,-1);
    vis[s]=0;
    while(!q.empty()){
        int here=q.front();q.pop();
        GO(i,here){
            if(vis[to[i]]==vis[here]+1) p[s][to[i]]=min(p[s][to[i]],p[s][here]);
            if(vis[to[i]]>=0) continue;
            p[s][to[i]]=p[s][here];
            if(s==here) p[s][to[i]]=to[i];
            vis[to[i]]=vis[here]+1;
            q.push(to[i]);
        }
    }
}
void add(int a,int b){
    nex[top1]=head[a];head[a]=top1;to[top1++]=b;
    nex[top1]=head[b];head[b]=top1;to[top1++]=a;
}
double dp(int a,int b){
    if(f[a][b]) return f[a][b];
    if(a==b) return 0;
    if(p[a][b]==b||p[p[a][b]][b]==b) return f[a][b]=1;
    double tot=dp(p[p[a][b]][b],b);
    GO(i,b){
        int tmp=to[i];
        tot+=dp(p[p[a][b]][b],to[i]);
    }
    f[a][b]=tot/(t[b]+1.0)+1.0;
    return f[a][b];
}
void input(){
    mem(head,-1);
    n=in();e=in();
    c=in();m=in();
    for(int i=1;i<=e;i++){
        int a,b;
        a=in();b=in();
        add(a,b);
        ++t[a];++t[b];
    }
    for(int i=1;i<=n;i++)
        bfs(i);
}
void xxj(){
    cout<<fixed<<setprecision(3)<<dp(c,m)<<endl;
}
int main(){
    input();
    xxj();
    return 0;
}

转载于:https://www.cnblogs.com/Xiaojianxiang/p/6644209.html