计蒜客--蒜头君王国 概率dp

参考自https://blog.csdn.net/nbl97/article/details/79968761 

#include <bits/stdc++.h>
#define MAXN 100005 
using namespace std;
int n; // 输入
double p;
double result[MAXN]; // i个节点联通的概率 
 
// 排列组合公式 : C(n,m)=n!/(m!(n-m)!)
int C(int n,int m) {
	int ret = 1;
	for(int i = n;i > m;i--) {
		ret *= i;
	}
	for(int i = n-m;i >= 1;i--) {
		ret /= i;
	}
	
	return ret;
}

// 快速幂 
double Pow(double x,int n) {
	double ret = 1.0;
	while(n) {
		if(n & 1) {
			// n为奇数
			ret *= x; 
		}		
		x *= x;
		n >>= 1;
	}
	
	return ret;
}

int main() {
	
	scanf("%d%lf",&n,&p);
	result[1] = 1;
	result[2] = p;
	for(int i = 3;i <= n;i++) {
		result[i] = 1.0;
		for(int k = 1;k < i;k++) {
			result[i] -= C(i-1,k-1)*result[k]*Pow(1.0-p,k*(i-k)) ;				
		}		
	}
		
	printf("%.3lf\n",result[n]);
		
	return 0;
}