vfk论文上的题
记
U=2n−1
U
=
2
n
−
1
,因为
∀T⊆SU|T=U
∀
T
⊆
S
U
|
T
=
U
我们把p看成集合幂级数,那么有
ans=∑∞k=1k∗(pkU−pk−1U)
a
n
s
=
∑
k
=
1
∞
k
∗
(
p
U
k
−
p
U
k
−
1
)
若我们定义集合幂级数f的第T项为
fT=∑∞k=1k∗(pkT−pk−1T)
f
T
=
∑
k
=
1
∞
k
∗
(
p
T
k
−
p
T
k
−
1
)
将两边同时做莫比乌斯变换(子集和变换)
FT=∑∞k=1k∗(PkT−Pk−1T)
F
T
=
∑
k
=
1
∞
k
∗
(
P
T
k
−
P
T
k
−
1
)
发现右边是个等比数列求和的东西,当
PT=1
P
T
=
1
时
FT=0
F
T
=
0
,否则
FT=11−PT
F
T
=
1
1
−
P
T
求出
F
F
的每一项后做一次莫比乌斯逆变换求出
code:
#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn = 1<<20;
const double eps = 1e-9;
int n,N;
double f[maxn];
void FWT(double s[],int sig)
{
for(int i=0;i<n;i++)
{
int t=1<<i;
for(int j=0;j<N;j++) if(!(j>>i&1))
f[j+t]+=sig*f[j];
}
}
int main()
{
//freopen("tmp.in","r",stdin);
//freopen("tmp.out","w",stdout);
scanf("%d",&n); N=1<<n;
for(int i=0;i<N;i++) scanf("%lf",&f[i]);
FWT(f,1);
for(int i=0;i<N;i++) f[i]=fabs(f[i]-1)<eps?0.0:-1.0/(1.0-f[i]);
FWT(f,-1);
if(fabs(f[N-1])<eps) puts("INF");
else printf("%.8lf\n",f[N-1]);
return 0;
}