P5242 题解

_Andy_L_

已于 2022-08-26 09:41:35 修改

阅读量58

点赞数 1

分类专栏：洛谷题解文章标签：算法 c++

于 2022-08-03 12:30:35 首次发布

本文链接：https://blog.csdn.net/lzh_Andy/article/details/126138069

版权

洛谷题解专栏收录该内容

30 篇文章 0 订阅

订阅专栏

数学题一个

题目传送门

首先经过计算，我们可以得知公式：

$\because s_i=\prod_{i=1}^n(1-p_k)$

$\text{又}\because f_i=\sum_{k=1}^i \dfrac{pk}{1-pk}$

$\therefore ans(i,j)=\frac{sr}{sl-1}fr-f_{l-1}$

就这样，立马想到了暴力，时间复杂度为 $O(n^2)$ ，但是突然想到了可以用尺取法来优化。

就这样，它变成了 $O (n)$ ，皆大欢喜！

$\text{Coding Time:}$

#include<bits/stdc++.h>

using namespace std;

#define gc getchar()

inline int read()//快读
{
    int x=0; 
    bool sgn=0;
    char s=gc;
    while(!isdigit(s))sgn|=s=='-',s=gc;
    while(isdigit(s))x=(x<<1)+(x<<3)+(s-'0'),s=gc;
    return sgn?-x:x;
}

int n;
long double num[1008611], cdy = 1.0, wzy, ans;//你想见祖宗吗？

int main() 
{
    n=read();
    for (register int i = 1, x; i <= n; ++i) 
    {
        x=read();
        num[i] = ((long double)x / 1000000.0);
        ans = max(ans, num[i]);
    }
    for (register int i = 1, tail = 1; i <= n; ++i) 
    {
        while (tail <= n && cdy * wzy < cdy * (1.0 - num[tail]) * (wzy + num[tail] / (1.0 - num[tail]))) 
        {
            cdy *= (1.0 - num[tail]);
            wzy += num[tail] / (1.0 - num[tail]);
            ++tail;
        }
        ans = max(ans, cdy * wzy);
        cdy /= (1.0 - num[i]);
        wzy -= num[i] / (1.0 - num[i]);
    }
    printf("%d\n", (int)(ans * 1000000));
    return 0;//完结，撒花！
}