Don't fake it till you made it

Fake it till you become it

Hackerrank Random Number Generator

题目链接:
https://www.hackerrank.com/challenges/random-number-generator-1
题目大意:
Pi>=0
1=Σni=1Pi
最大化 Σni=1Pi(1Pi)i

一开始的感觉很对 但是发现拉格朗日出现负数就不会处理
后来问了uoj的大爷们才发现自己对拉格朗日理解不深。。。

套用拉格朗日
f(x)=Σni=1Pi(1Pi)i
φ(x)=Σni=1Pi1
F(x,φ)=Σni=1Pi(1Pi)iλ(Σni=1Pi1)

F(x,φ)xi=i2Piiλ=0
iλ=2Pii
Pi=0.5λ/(2i)

Pi=0.5λ/(2i)
Σni=1Pin0.5λ0.5(Σni=11/n)=1
λ=(n2)/(Σni=11/n)
λn/ln(n)
然后就会出现Pi<0
考虑为什么会出现这种情况 感性理解后发现是由于有权重 i 的存在,使越靠后的Pi越接近于0.5(上面推出来的对称轴) 那么我们需要使这部分小于0的Pi 变成0 同时满足此时目标最小
继续套用拉格朗日 二分找到一个分界点mid 使任意i<midPi=0 并满足上述条件
然后就是O(Tlogn)

#include<cstdio>
#include<cstring>
using namespace std;

double f[1000001];


int main()
{
    int n;
    n=1000000;
    int i;
    for(int i=1;i<=n;i++)f[i]=f[i-1]+1.0/i;
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        if(n==1){puts("0.00000000000");continue;}
        if(n==2){puts("0.50000000000");continue;}
        int l=1,mid,r=n;
        while(l<r)
        {
            mid=l+r>>1;
            double F=(n-mid-1)/(f[n]-f[mid-1]);
            if(0.5-F/2/mid>=0)r=mid;
            else l=mid+1;
        }
        r=n;
        double F=(r-l-1)/(f[r]-f[l-1]);
        double As=(r-l+1)*1ll*(r+l)/4.0-F*F/2*(f[r]-f[l-1]);
        printf("%.10f\n",As/(2*n-1));
    }
    return 0;
}
阅读更多
版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/liutian429073576/article/details/53541886
个人分类: 拉格朗日乘数法
想对作者说点什么? 我来说一句

随机函数发生器_C程序

2013年05月10日 1011B 下载

CUDA 分子动力学模拟相关

2008年09月10日 11KB 下载

Random Number Generator v1.0

2008年06月09日 7KB 下载

没有更多推荐了,返回首页

不良信息举报

Hackerrank Random Number Generator

最多只允许输入30个字

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭