hdu 4336 Card Collector——最值反演

#include<cstdio>
#include<cstring>
#include<algorithm>
#define db double
using namespace std;
const int N=25,M=(1<<20)+5;
int n,ct[M],bin[N];db p[N],g[M];
void init()
{
  int n=20;
  bin[0]=1;for(int i=1;i<=n;i++)bin[i]=bin[i-1]<<1;
  for(int s=1;s<bin[n];s++)ct[s]=ct[s-(s&-s)]+1;
}
int main()
{
  init();
  while(scanf("%d",&n)==1)
    {
      for(int i=1;i<=n;i++)scanf("%lf",&p[i]);
      for(int s=1;s<bin[n];s++)
    {
      g[s]=0;
      for(int i=1;i<=n;i++)if(s&bin[i-1])g[s]+=p[i];
      g[s]=1/g[s];
    }
      db ans=0;
      for(int s=1;s<bin[n];s++)
    ans+=(ct[s]&1?1:-1)*g[s];
      printf("%.10f\n",ans);
    }
}