题解:非常简单。
第i+1道题作对的概率,是
1)第i道题选择了1~a(i+1)中的一个选项,并且作对了i+1题。
2)若ai>a(i+1) 第i道题选择了a(i+1)+1~ai中的一个选项并做对了i+1题。显然这部分概率为0。
如果ai>a(i+1),那么“第i道题选择了1~a(i+1)中的一个选项”的概率是ai/a(i+1)。
否则ai<=a(i+1),“第i道题选择了1~a(i+1)中的一个选项”的概率就是1.
综合起来就是min(ai,a(i+1))/ai。
在“第i道题选择了1~a(i+1)中的一个选项”的前提下,作对i+1题的概率为1/(a(i+1))。
那么根据贝叶斯公式,P(作对i+1题)=P(第i道题选择了1~a(i+1)中的一个选项)*P(在第i道题选择了1~a(i+1)中的一个选项的前提下作对第i+1题)
因此P(i+1)=min(ai,a(i+1))/ai*1/a(i+1)=min(ai,a(i+1))/(ai*a(i+1))
经过化简P(i+1)=1/max(ai,a(i+1))。
所以O(n)扫一遍即可。
//BZOJ 2134
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define lint long long
#define mod 100000001
#define MAXN 10000010
using namespace std;
int a[MAXN];
double ans;
int main()
{
int n,A,B,C;scanf("%d%d%d%d%d",&n,&A,&B,&C,&a[1]);
for(int i=2;i<=n;i++) a[i]=((lint)a[i-1]*A+B)%mod;
for(int i=1;i<=n;i++) a[i]=a[i]%C+1;a[n+1]=a[1];
for(int i=1;i<=n;i++) ans+=1.0/max(a[i],a[i+1]);
printf("%.3f\n",ans);return 0;
}