J Heritage of skywalkert
https://www.cnblogs.com/longl/p/9424657.html
思路:
调用题目的给的函数tang(), 生成数a[1]到a[n], 然后调用STL的库函数nth_element(原理是快排,也可以自己手搓) ,选出前100大的数,然后暴力选两个数,使得最小公倍数最大。
关于STL中的nth_element()方法的使用:通过调用nth_element(start, start+n, end) 方法可以使第n大元素处于第n位置(从0开始,其位置是下标为 n的元素),并且比这个元素小的元素都排在这个元素之前,比这个元素大的元素都排在这个元素之后,但不能保证他们是有序的。
代码:
#include<bits/stdc++.h>
using namespace std;
const int maxn =1e7+20;
unsigned x,y,z,A,B,C;
typedef unsigned long long ull;
unsigned a[maxn];
unsigned tang(){
unsigned t;
x^=x<<16;
x^=x>>5;
x^=x<<1;
t=x;
x=y;
y=z;
z=t^x^y;
return z;
}
int main()
{
int t,n;
cin>>t;
for(int cas=1;cas<=t;cas++){
cin>>n>>A>>B>>C;
x=A,y=B,z=C;
for(int i=1;i<=n;i++)
a[i]=tang();
vector<unsigned> vec;
int num=min(n,100);
nth_element(a+1,a+1+n-num,a+1+n); ///只要排一次
/// 选出前一百大的数
for(int k=1;k<=num;k++){
vec.push_back(a[n+1-k]);
}
///暴力找lcm最大
ull ans=0;
for(int i=0;i<vec.size();i++){
for(int j=i+1;j<vec.size();j++){
ull gcd = __gcd((ull)vec[i],(ull)vec[j]);
ull tp=vec[i]/gcd*vec[j];
ans = max(ans,tp);
}
}
printf("Case #%d: ",cas);
cout<<ans<<endl;
}
return 0;
}
深信服公司福利 906人发布