BZOJ2134 单选错位

传送门

题解：非常简单。

第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;
}