UVA 11600 Masud Rana

  题目大意:有一个n个点的完全图,有些路上有妖怪。现在一个人从一号点出发,每天随机走向另一个点,消灭路上的妖怪,问平均几天后所有点之间存在没有妖怪的路径。点数≤30。

  看到点这么少肯定状压,看见存不下肯定map,事实上这题数据不够强力,很多复杂度不对的东西都可以艹过去(比如我的算法)。

  首先可以缩点,把本来存在路径的点缩起来。

  设f(x)表示:当安全的点集为x时,处理完所有的点的期望次数。

  那么可以写出一个解不出来的式子:

  

  p表示走到自己集合内的概率,pi表示走到另一个点的概率。

  解不出来是因为存在自己更新自己。怎么办呢?正常人都会想着化式子吧。

 

  

  

  

  

  发现每次处理,集合都会变大。

  这个时候打一个记忆搜就可以了。边界状态是f[2^n-1]=0.00。

#include    <iostream>
#include    <cstdio>
#include    <cstdlib>
#include    <algorithm>
#include    <map>
#include    <cstring>
#include    <queue>
#include    <complex>
#include    <stack>
#define LL long long int
#define dob double
#define FILE "11600"
using namespace std;

const int N = 110;
map<int,double>Mp;
int n,m,fa[N],bin[N],U[N];

inline int gi(){
  int x=0,res=1;char ch=getchar();
  while(ch>'9'||ch<'0'){if(ch=='-')res*=-1;ch=getchar();}
  while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
  return x*res;
}

inline int find(int x){
  return x==fa[x]?x:fa[x]=find(fa[x]);
}

inline double dfs(int S){
  if(Mp.find(S)!=Mp.end())return Mp[S];
  double Ans=0.0,sz=0.0;
  for(int i=1;i<=n;++i){
    if(S&bin[i-1]){sz+=1.0;continue;}
    Ans+=dfs(S|U[i]);
  }
  Ans=(1.0*n-1.0)/(1.0*n-sz)*(Ans/(1.0*n-1.0)+1.0);
  return Mp[S]=Ans;
}

inline void solve(){
  Mp.clear();bin[0]=1;
  n=gi();m=gi();if(n==1){puts("0");return;}
  for(int i=1;i<=n;++i)fa[i]=i,bin[i]=bin[i-1]<<1,U[i]=bin[i-1];
  for(int i=1;i<=m;++i){
    int f1=find(gi()),f2=find(gi());
    if(f1!=f2)fa[f2]=f1,U[f1]|=U[f2];
  }
  for(int i=1;i<=n;++i)U[i]=U[find(i)];
  Mp[bin[n]-1]=0.0;printf("%.10lf
",dfs(U[1]));
}

int main(){
  freopen(FILE".in","r",stdin);
  freopen(FILE".out","w",stdout);
  int Case=gi();
  for(int t=1;t<=Case;++t)
    printf("Case %d: ",t),solve();
  fclose(stdin);fclose(stdout);
  return 0;
}
Masud Rana

 

 

原文地址:https://www.cnblogs.com/fenghaoran/p/7677601.html