题目
原题链接:点这里
题解
算法思路
- 还是没有发散性思维,比赛时候没有想到这个
- 很明显的随机数函数,还在纠结有没有规律,这个东西一旦敲定,就要下定论,别把时间花在无意义的地方
- 质数的分布密度……
代码实现
#include<bits/stdc++.h>
#include<unordered_map>
#include<unordered_set>
using namespace std;
#define el '\n'
#define cl putchar('\n')
#define pb push_back
#define eb emplace_back
#define fir first
#define sec second
#define int long long
typedef long long ll;
typedef pair<int,int> pii;
typedef vector<int> vci;
typedef map<int,int> mii;
typedef mii::iterator mii_it;
const int N=1e4+10,M=6e7+10;
int T,n,m,x,y,k;
int a[N];//a也需要开longlong
int L=1,R=10000;
unsigned long long seed;
//这个生成a[i]的函数,没有规律,就是一个真的随机数
unsigned long long xorshift(){
unsigned long long x=seed;
x^=x<<13;
x^=x>>7;
x^=x<<17;
return seed=x;
}
int gen(){
return xorshift()%(R-L+1)+L;
}
struct eage{
int u,v,c;
}e[M];
int fa[N];
int find(int x){
if(x==fa[x])return x;
else return fa[x]=find(fa[x]);
}
int cmp(eage x,eage y){
return x.c<y.c;
}
int gcd(int a,int b){
return !b?a:gcd(b,a%b);
}
signed main() {
cin.tie(0);
cout.tie(0);
cin>>n>>L>>R>>seed;
for(int i=1;i<=n;i++){
a[i]=gen();
}
if(L==R)cout<<R*(n-1);//所以边都相同,这个需要放在最前面(以后记得考虑if else的顺序)
else if(n>10000)cout<<n-1;//数量那么大,很难保证随机数里面不生成质数(质数密度),所以每条边跟质数都是1
else{//n小的时候就直接kruscal算答案
long long ans=0,cnt=1;
for(int i=1;i<=n;i++){
for(int j=i+1;j<=n;j++){
k++;
e[k].u=i;
e[k].v=j;
e[k].c=gcd(a[i],a[j]);
}
}
for(int i=1;i<=n;i++)fa[i]=i;
sort(e+1,e+k+1,cmp);
for(int i=1;i<=k;i++){
int x=find(e[i].u);
int y=find(e[i].v);
if(x!=y){
fa[x]=y;
cnt++;
ans+=e[i].c;
if(cnt==n)break;
}
}
cout<<ans;
}
}