Luogu P1291 [SHOI2002]百事世界杯之旅 // 易错的期望

最新推荐文章于 2022-06-10 16:06:44 发布

weixin_30867015

最新推荐文章于 2022-06-10 16:06:44 发布

阅读量59

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/SINXIII/p/10999041.html

本文纠正了一个关于计算拿到明星数量概率的错误公式，并提供了一个由BTC大神提出的正确算法。通过分析，展示了如何使用C++实现该算法，包括长整数运算、最大公约数计算及分数的简化。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

首先有一个式子是错的

网上有个式子是f[i] = i/n * f[i] + n-i/n * f[i-1] + 1.f[i]表示现在拿到了几个明星

这个A了但是是错的。因为只要买除i-1以外的任意一个就行，所以后面的应该是frac{n-(i-1)}{n}f[i-1]

BTC大神通过倒推写出了一个正确的式子，luogu上ButterflyDew的题解也不错。

美妙输出。

#include<bits/stdc++.h>

using namespace std ;

long long n,a,b;

long long gcd(long long a,long long b){
    return b ?  gcd(b,a%b) : a;
} 

int main(){
    scanf("%lld",&n);
    b=1;
    for(int i=1;i<=n;i++){
        b = b * i / gcd (b,i);
    }
    for(int i=1;i<=n;i++){
        a += b / i;
    }
    a *= n;
//    cout<<a<<endl<<b<<endl;
    long long p = a/b;
    a %= b;
    if(a==0){
        printf("%lld\n",p);
        return 0;
    }
    long long ya = a;
    a/=gcd(a,b);
    b/=gcd(ya,b);
//    cout<<p<<endl<<a<<endl<<b<<endl;
    long long len1=0,len2=0,sb = b,sp = p;
    while(sp){
        ++len1;
        sp/=10;
    }
    while(sb){
        ++len2;
        sb/=10;
    }
//    cout<<len1<<endl<<len2;
    for(int i=1;i<=len1;i++){
        printf(" ");
    }
    printf("%lld",a);
    printf("\n");
    if(p!=0) printf("%lld",p);
    for(int i=1;i<=len2;i++){
        printf("-");
    }
    printf("\n");
    for(int i=1;i<=len1;i++){
        printf(" ");
    }
    printf("%lld",b);
    return 0;
}