题目链接: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 都可以
}
}