B - Build Roads(Gym - 103118B )完全图求最小生成树

本题给了一个生成两点间距离的函数,每次赋对一个点赋值,两边的边长位顶点值的gcd

int xorshift64()
{
	int x = seed;
	x ^= x << 13;
	x ^= x >> 7;
	x ^= x << 17;
	return seed = x;
}
int gen()
{
	return xorshift64() % (r - l + 1) + l;
}

 尝试几组 l r seed 值会发现,求得的边长大多数为1,l==r时边长为 l,所以在数据量小时跑krusal求最小生成树,数据量大时可以直接根据规律求出答案。

#include<iostream>
#include<queue>
#include<algorithm>
#include<cstring>
#include<map>
#include<set>
const int N = 1e6 + 6006;
const int M = 4 * N;
#define mod %998244353 
#define FAST ios::sync_with_stdio(false);
#define endl '\n'
#define int unsigned long long 
#define rep(i,n)for(int i=1;i<=n;++i)
#define sc(n) scanf("%d",&n);
using namespace std;
int n, l, r, seed, a[200001];
int  b[N], k[N], p[N];
struct node
{
	int u;
	int v;
	int w;
}edge[N];
typedef pair<int, int>PII;
priority_queue<PII, vector<PII>, greater<PII>>Q;
int find(int x)
{
	if (x != p[x])p[x] = find(p[x]);
	return p[x];
}
int gcd(int a, int b)
{
	return b > 0 ? gcd(b, a % b) : a;
}
int xorshift64()
{
	int x = seed;
	x ^= x << 13;
	x ^= x >> 7;
	x ^= x << 17;
	return seed = x;
}
int gen()
{
	return xorshift64() % (r - l + 1) + l;
}
int cnt = 0;
bool cmp(node a, node b)
{
	return a.w < b.w;
}
int kruskal()
{
	int res = 0;
	sort(edge + 1, edge + 1 + cnt, cmp);
	int k = 0;
	for (int i = 1; i <= cnt; i++)
	{
		if (k == n - 1)break;
		int a = edge[i].u, b = edge[i].v;
		if (find(a) != find(b))
		{
			res += edge[i].w;
			p[find(a)] = find(b);
			k++;
		}
	}
	return res;
}
signed main()
{
	FAST;
	cin >> n >> l >> r >> seed;
	if (n >= 5000)
	{
		if (l >= r)
		{
			cout << (n - 1) * l << endl;
			return 0;
		}
		else
		{
			cout << (n - 1) << endl;
			return 0;
		}
	}
	for (int i = 1; i <= n; i++)p[i] = i;
	for (int i = 1; i <= n; i++)
		a[i] = gen();
	for (int i = 1; i <= n; i++)
	{
		for (int j = i + 1; j <= n; j++)
		{
			cnt++;
			edge[cnt].u = i;
			edge[cnt].v = j;
			edge[cnt].w = gcd(a[i], a[j]);
		}
	}
	int res = kruskal();
	cout << res << endl;
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值