bzoj 2337 [HNOI2011]XOR和路径【高斯消元+dp】

首先,我们发现,因为是无向图,所以相连的点之间是有“依赖性”的,所以不能直接用dp求解。

因为是xor,所以按位处理,于是列线性方程组,设$ x[i] $为点i到n异或和为1的期望,因为从1到n和从n到1一样,所以选择倒着推,即,

if(deg[e[i].va]==0)

\[x[u]=\sum_{v}^{v\subset son(u)}\frac{x[v]}{deg[i]}
\]

else

\[x[u]=\sum_{v}^{v\subset son(u)}\frac{1-x[v]}{deg[i]}
\]

列n元n次方程组高斯消元求解即可

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=105,M=100005;
int n,m,h[N],cnt,in[N];
double f[N][N],ans;
struct qwe
{
int ne,to,va;
}e[M<<1];
void add(int u,int v,int w)
{
cnt++;
in[u]++;
e[cnt].ne=h[u];
e[cnt].to=v;
e[cnt].va=w;
h[u]=cnt;
}
void gaosi()
{
for(int i=1;i<=n;i++)
{
int id=i;
double mx=0.0;
for(int j=i;j<=n;j++)
if(fabs(f[j][i])>mx)
id=j,mx=fabs(f[j][i]);
if(id!=i)
for(int j=1;j<=n+1;j++)
swap(f[id][j],f[i][j]);
double t=f[i][i];
for(int j=1;j<=n+1;j++)
f[i][j]/=t;
for(int j=1;j<=n;j++)
if(j!=i)
{
double t=f[j][i];
for(int k=1;k<=n+1;k++)
f[j][k]-=t*f[i][k];
}
}
}
int main()
{
scanf("%d%d",&n,&m);
for (int i=1,x,y,z;i<=m;i++)
{
scanf("%d%d%d",&x,&y,&z);
add(x,y,z);
if (x!=y)
add(y,x,z);
}
for(int i=0;i<=30;i++)
{
memset(f,0,sizeof(f));
for(int u=1;u<=n-1;u++)
{
f[u][u]=1.0;
for(int j=h[u];j;j=e[j].ne)
{
if(e[j].va&(1<<i))
f[u][e[j].to]+=1.0/in[u],f[u][n+1]+=1.0/in[u];
else
f[u][e[j].to]-=1.0/in[u];
}
}
f[n][n]=1.0;
gaosi();
ans+=(f[1][n+1])*(1<<i);
}
printf("%.3lf\n",ans);
return 0;
}
上一篇:【BZOJ】2337: [HNOI2011]XOR和路径 期望+高斯消元


下一篇:●BZOJ 2337 [HNOI2011]XOR和路径