G - Research Productivity Index（期望dp）

作为一个刚入门的蒟蒻，在神仙队友的帮助下，此处艾特含哥@会飞的小蛇 写出了第一个期望dp，写此博客，引以纪念。

题目传送门

题目大概意思是有n篇论文，每一篇通过的概率给出，提交b篇通过a篇的得分是a^a/b，特殊的，如果a=0，那么得分是0。
输入两行，第一行是论文总数n，第二行有n个整数，第i个整数表示第i篇论文通过的概率，以概率的100倍给出。

显然，我们肯定是“择优选取”，即通过概率大的先拿。

初始化
初始化的时候我们需要找到一个边界条件，即已知或者可推理得到的值。
这个题，我们用dp数组去存概率，dp[i][j]里存的是提交i篇通过j篇的值，那么我们可以推理得到，dp[0][0]=1；这是显而易见的，提交0篇通过0篇，他是肯定成立的，然后我们枚举一遍i，把提交i篇但是都不通过的值算出。

    dp[0][0]=1;
    for(int i=1;i<=n;i++)
    {
        dp[i][0]=dp[i-1][0]*(100-a[i])*1.0/100;
    }

状态转移
这题的难点在于计算概率，我们对于提交i篇通过j篇，可以由两种状态转移过来，那就是之前交了i-1篇，通过了j-1篇，现在再交一篇且通过、之前交了i-1篇，通过了j篇，现在再交一篇，没通过。这里需要注意的是后者的一个边界，就是当我需要计算交了i篇通过i篇的时候，第二个状态是交i-1篇通过了i篇，这是不符合常理的，然而，我们用一个二维数组去存概率的时候并不用考虑这个，因为在上文中说的初始化和上文提到的状态转移中，都不会改变到类似dp[i-1][i]这种值，那么它就是0，所以对于我们的状态转移并无影响。

    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=i;j++)
        {
            dp[i][j]=dp[i-1][j-1]*a[i]*1.0/100+dp[i-1][j]*(100-a[i])*1.0/100;
        }
    }

最后就是我们枚举提交的篇数，然后计算期望啦。

    double res=0;
    for(int i=1;i<=n;i++)
    {
        double temp=0;
        for(int j=1;j<=i;j++)
        {
            temp+=(dp[i][j]*pow(j,j*1.0/i));
        }
        res=max(res,temp);
    }

蒟蒻AC代码

#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<string.h>
#include<map>
#include<cmath>
#include<queue>
#include<vector>
using namespace std;
#define int long long
const int N=1e3+5;
#define ll long long
int a[N];
double dp[N][N];
bool cmp(int a,int b)
{
    return a>b;
}
signed main()
{
    int n;
    cin>>n;
    for(int i=1;i<=n;i++)
        cin>>a[i];
    sort(a+1,a+n+1,cmp);
    dp[0][0]=1;
    for(int i=1;i<=n;i++)
    {
        dp[i][0]=dp[i-1][0]*(100-a[i])*1.0/100;
//        dp[i][i]=dp[i-1][i-1]*a[i]*1.0/100;
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=i;j++)
        {
            dp[i][j]=dp[i-1][j-1]*a[i]*1.0/100+dp[i-1][j]*(100-a[i])*1.0/100;
        }
    }
    double res=0;
    for(int i=1;i<=n;i++)
    {
        double temp=0;
        for(int j=1;j<=i;j++)
        {
            temp+=(dp[i][j]*pow(j,j*1.0/i));
        }
        res=max(res,temp);
    }
    printf("%.9lf\n",res);
}