【题解】sdoj2605闲荡（2018-08-15集训T3） 拓扑排序+概率DP+异或

题目链接

问题描述

L 饭后无聊，便在 BugTown 里闲荡。
BugTown 共有 N 栋房屋和 M 条有向道路。每栋房屋都有一个非负整数 vi 作为标
识。
BugTown 有一个特性十分神奇：从任意一个房屋离开后沿着路走再也不会回到原
地。
L 想选一个房屋作为闲荡的起点，之后，他会随机选择一条当前位置能走的道路顺
其 走过去，如此反复直到没有能走的道路。
由于极度无聊， L 发明了一个游戏以为消遣。他在闲荡的过程中记录已经过的房屋
标 识的异或和（含起点）。闲荡完后，他会得到一个数。
L 希望对每个房屋算出以它为起点能得到的数的期望值，但是他不知道怎么算，只
好 求助于你。

输入

第一行两个正整数 N; M 分别为房屋数和道路数。
第二行 N 个数为 vi。
接下来 M 行每行两个整数 ai; bi 描述一条 ai 到 bi 的道路

输出

输出共 N 行。第 i 行一个实数表示以 i 号房屋为起点时最后得到的数的期望值。 实数
四舍五入到小数点后三位。

输入样例

3 2
1 2 3
1 2
2 3

输出样例

0.000
1.000
3.000

数据规模

对于 10% 的数据， N <= 5; M <= 10。
对于 30% 的数据， N, M<= 50。
对于 70% 的数据， N, M <=1000。
对于 100% 的数据， 1 <= N,M <= 1e5; 0 <= vi <= 1e9。

---

我求异或和期望求错了……orz太蒻。先去练10道概率DP，要带异或的那种

#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
using namespace std;
int n,m;
const int N=1e5+10;
int sign[N];
int deg[N];
vector<int>vx[N];
queue<int>q;
double dp[N][30];
double ans[N];
int s[N];
int t=0;
void toposort()
{
    for(int i=1;i<=n;i++)if(!deg[i])q.push(i);
    while(!q.empty())
    {
        int u=q.front();q.pop();
        s[++t]=u;
        for(int i=0;i<vx[u].size();i++)
        {
            int v=vx[u][i];
            if(--deg[v]==0)q.push(v);
        }
    }
}
void getans()
{
    for(int i=0;i<30;i++)
        for(int j=t;j>=1;j--)
        {
            int u=s[j];
            int a=((sign[u]>>i)&1);
            if(!vx[u].size())
                dp[u][a]=1.0,dp[u][a^1]=0.0;
            else
            {
                double p=1.0/vx[u].size();
                dp[u][0]=dp[u][1]=0.0;//漏了这句 
                for(int l=0;l<vx[u].size();l++)
                {
                    int v=vx[u][l];
                    dp[u][0]+=dp[v][a]*p;
                    dp[u][1]+=dp[v][a^1]*p;//注意此处写法 
                }
            }
            ans[u]+=dp[u][1]*(1<<i);
        }
}
int main()
{
    //freopen("walk.in","r",stdin);
    //freopen("walk.out","w",stdout);
    //freopen("in.txt","r",stdin);
    scanf("%d%d",&n,&m);
    int u,v;
    for(int i=1;i<=n;i++)scanf("%d",&sign[i]);
    for(int i=1;i<=m;i++)
        scanf("%d%d",&u,&v),vx[u].push_back(v),deg[v]++;
    toposort();
    getans();
    for(int i=1;i<=n;i++)
    printf("%.3lf\n",ans[i]);
    return 0;
}