牛客网暑期ACM多校训练营（第六场）- J Heritage of skywalkert (随机概率)

原创于 2018-08-06 21:08:00 发布 · 324 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

其他同时被 2 个专栏收录

106 篇文章

订阅专栏

暴力

11 篇文章

订阅专栏

本文介绍了一种解决求最大最小公倍数(LCM)问题的方法。通过记录前30个最大的随机数，并使用快速选择算法(nth_element)来优化搜索过程，最终找到两数使得它们的LCM最大。

https://www.nowcoder.com/acm/contest/144/J

题意：

给你一个操作，每次操作会产生一个a[i]。

让你在n次操作之后，找到两个a[i]，使得他们的lcm最大。

POINT：

这个操作相当于随机数，我们找出前30个最大数。

然后暴力找。

可以用nth_element

#include <stdio.h>
#include <iostream>
#include <string.h>
#include <algorithm>
#include <queue>
#include <cmath>
using namespace std;
#define LL unsigned long long
const int N = 15;
unsigned int x,y,z;
 
 
unsigned int tang()
{
    x^=x<<16;
    x^=x>>5;
    x^=x<<1;
    LL t=x;
    x=y;
    y=z;
    z=t^x^y;
    return z;
}
 
LL gcd(LL x,LL y)
{
    return y==0?x:gcd(y,x%y);
}
 
LL f(LL x,LL y)
{
    return 1LL*x/gcd(x,y)*y;
}
LL a[31];
 
int main()
{
 
    int T;
    int n;
    cin>>T;
    int cas=0;
    while(T--){
        cin>>n>>x>>y>>z;
        LL ans=0;
        for(int i=1;i<=30;i++)
            a[i]=0;
        for(int i=1;i<=n;i++){
            tang();
            if(z>a[1]){
                a[1]=z;
                sort(a+1,a+31);
            }
        }
        for(int i=1;i<=30;i++){
            for(int j=i;j<=30;j++){
                if(a[j]==0||a[i]==0) continue;
                    ans=max(ans,f(a[i],a[j]));
            }
        }
 
        printf("Case #%d: %llu\n",++cas,ans);
 
    }
 
    return 0;
}