Scout YYF I POJ - 3744

题目链接：http://poj.org/problem?id=3744

感觉很有必要先说下编译器的问题，以下链接主要说下g++，c++编译器在输入输出的不同：

https://blog.csdn.net/febr2/article/details/52068357

这个博客里说的很清楚，感谢大佬。

再者相同的编译器在不同OS里会有不同的情况。

说下题目吧：

概率dp+矩阵快速幂优化。

如果要安全走过每个mine，那就必须走到mine（第i个地雷）后面的位子，因为你一次只能走一个单位或两个单位。

我们用array[i]表示i-th mine 的位子，f(array[i]-1)表示安全到达array[i]-1的概率，则走过i的概率：dp[i]=（1-p）*(f[array[i]-1]);

接下来就是求 f(array[i]-1):

 我们可以看成从array[i-1]+1到array[i]-1的概率与dp[i-1]的积，模拟斐波那契大数的求法：

 f(n)=q*（f(n-1)）+(1-q）*f(n-2);

 f(n-1)=1*f(n-1)+0*f(n-2);

•            f(n)              q   ,   1-q          f(n-1)
•   =>                =                        *                 

                 f(n-1)           1  ，  0            f(n-2)

接下快速矩阵幂 就行了；

还有条更简单的思路就是你可以把整个路分成n段分别求每段的安全走过的概率 ，再乘起来就可以的。

ACode：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<iomanip>
#define LL long long 
using namespace std;
const int M=10e7+4;
bool divs[M];
double p;
double fast_pow(double p1,int n)
{
    double a[2][2]={p,1-p,1.0,0.0},f[2]={p1,0.0},temp[2][2];
	while(n)
	{
		if(n&1)
		{
			double t0=f[0],t1=f[1];
			f[0]=a[0][0]*t0+a[0][1]*t1;
			f[1]=a[1][0]*t0+a[1][1]*t1;
		}
		n>>=1;
		for(int i=0;i<2;++i)
		{
			for(int j=0;j<2;++j)
			{
			    temp[i][j]=a[i][j];	
			}
		}
		for(int i=0;i<2;++i)
		{
			for(int j=0;j<2;++j)
			{
				a[i][j]=0.0;
				for(int k=0;k<2;++k)
				{
					
					a[i][j]+=temp[i][k]*temp[k][j];
					 
				}
			}
		}
		//cout<<endl;
	}
	return f[0];
} 
int main()
{
	int n;
	double p1,dp[15];
    LL  array[15]; 
	while(scanf("%d%lf",&n,&p)!=-1)
	{
	    for(int i=1;i<=n;++i)
	    {
	    	scanf("%lld",&array[i]);
		}
		sort(array+1,array+n+1);
		if(array[1]==1.0)
		{
			printf("%0.7lf\n",0.0);
			continue;
		}
		dp[0]=1.0;
		array[0]=0;
		for(int i=1;i<=n;++i)
		{
		    if(array[i]==array[i-1])
		    {
		    	dp[i]=dp[i-1];
		    	continue;
			}
			if(array[i]==array[i-1]+1)
			{
				dp[n]=0.0;
				break;
			}
			dp[i]=(1-p)*fast_pow(dp[i-1],array[i]-2-array[i-1]);
		}
		printf("%0.7lf\n",dp[n]);//c++可以但g++不能，%f 都可以
	}
}