幻象【NOIP2016提高A组模拟9.4】

题目

phantom是一位爱思考的哲♂学家。
最近phantom得到了森の妖精的真传。在他练功的时候, 每秒他的思绪中都有一定的概率浮现出奇♂异的幻象,持续x秒的幻象将产生x^2 的幻象值。
phantom练功发自真心,他想知道,在N秒内他期望产生的幻象值是多少.

样例输入：
第一行包含 1 个正整数 N ，表示总时间 N 秒。
第二行包含 N 个用空格隔开的在[0,100]之间的正整数，其中第i个数a[i]表示第i秒浮现幻象的概率为百分之a[i]。
3
50 50 50

样例输出：
1 个实数，四舍五入后保留一位小数，表示期望幻象值。
2.8

数据范围：
对于 40%的数据 N ≤ 10
对于 60%的数据 N ≤ 100
对于 100%的数据，N ≤ 10^6
数据规模较大，请使用效率较高的读入方式。

---

剖解题目

。。。。。。

---

思路

比赛时题目不完整，那个x^2的没有打出来，然后就一脸懵逼。实际上就算是弄出来也是差不多，毕竟实在是对期望类型的题是不熟悉的。

---

解法

40%：暴力。
60%：奇怪的dp乱推。
100%：设L[i]为第i秒幻象的持续时间的期望. 显然L[i] = (L[i-1] + 1) * a[i]% (百分号)
设f[i]表示前i秒的答案 f[i] = f[i-1] + ((L[i-1] + 1)^2 – L[i-1]^2) * a[i]%
时间复杂度为O(N)

---

代码

#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)

using namespace std;

const int maxn=1e6+2;
int n,a[maxn];
double f[maxn],L[maxn];

double sqr(double x){
    return x*x;
}
int main()
{
    scanf("%d",&n);
    fo(i,1,n) {
        scanf("%d",&a[i]);
        L[i]=(L[i-1]+1)*(double)a[i]/(double)100;
        f[i]=f[i-1]+(sqr(L[i-1]+1)-sqr(L[i-1]))*a[i]/(double)100;   
    }
    printf("%.1lf",f[n]);
}