HDU6237 A Simple Stone Game (思维 数学)

题目链接: A Simple Stone Game

大致题意

给定 n n n堆石子, 第 i i i堆有 a i a_i ai个.

你可以执行任意次操作, 每次选择编号为 i i i j j j的石子堆, 表示把一个石子从第 i i i堆挪到第 j j j堆.

问: 至少执行多少次操作, 才能使得 g c d ( { a 1 , a 2 , . . , a n } ) > 1 gcd(\{ a_1, a_2, .., a_n \}) > 1 gcd({a1,a2,..,an})>1

解题思路

思维

我们简化一下题意, 相当于最后 a 1 , a 2 , . . , a n a_1, a_2, .., a_n a1,a2,..,an都有公因子 x x x. ( x x x可以是 g c d ( { a 1 , a 2 , . . , a n } ) gcd(\{ a_1, a_2, .., a_n \}) gcd({a1,a2,..,an})的因子)

可以证明, 我们只需要枚举所有的 x ∈ { 质 数 } x \in \{ 质数 \} x{} 即可.

因为如果公因子 x x x不是质数, 可以写成 x = p 1 k 1 ⋅ p 2 k 2 . . . x = p_1^{k_1} · p_2^{k_2}... x=p1k1p2k2...的形式, 我们枚举 p i p_i pi同样是合法的, 且是更优的.


s u m = ∑ i − 1 n a i sum = \displaystyle\sum_{i - 1}^na_i sum=i1nai, 我们只需要检查 s u m sum sum的所有质因子即可.

对于每个质因子 p p p而言, 我们可以设 b [ ] b[] b[], 其中 b i = a i b_i = a_i % p bi=ai, 表示 a i a_i ai需要移除多少个石子才能有因子 p p p.

我们较容易想到, 如果对于某个 b i b_i bi而言, 如果比较小, 我们从这堆移出石子到其他堆更优, 反之应当把石子移入当前堆.

因此我们不妨对于 b [ ] b[] b[]从小到大排序, 用双指针 l = 1 , r = n l = 1, r = n l=1,r=n来表示, 我需要从第 l l l堆移石子到第 r r r堆.

可以保证最后一定可以使得 ∀ b i % p = 0 \forall b_i \% p = 0 bi%p=0.

因为sum % p = 0, 则 ∑ i − 1 n b i   %   p = 0 \displaystyle\sum_{i - 1}^nb_i \ \% \ p = 0 i1nbi % p=0 必然成立.

AC代码

#include <bits/stdc++.h>
#define rep(i, n) for (int i = 1; i <= (n); ++i)
using namespace std;
typedef long long ll;
const int N = 1E5 + 10;
int n;
int a[N], b[N];
ll fact(ll p) {
	ll all = 0;
	rep(i, n) b[i] = a[i] % p, all += b[i];
	sort(b + 1, b + 1 + n);

	ll res = 0;
	for (int l = 1, r = n; l < r; ++l) {
		if (!b[l]) continue;

		while (b[l]) {
			ll qaq = min<ll>(b[l], p - b[r]);
			b[l] -= qaq, b[r] += qaq;
			res += qaq;

			if (b[r] == p) r--;
		}
	}
	return res;
}
int main()
{
	int t; cin >> t;
	while (t--) {
		scanf("%d", &n);
		ll all = 0;
		rep(i, n) scanf("%d", &a[i]), all += a[i];

		ll res = LLONG_MAX;
		for (int p = 2; p <= all / p; ++p) {
			if (all % p == 0) {
				while (all % p == 0) all /= p;

				res = min(res, fact(p));
			}
		}
		if (all > 1) res = min(res, fact(all));

		printf("%lld\n", res);
	}
	return 0;
}

END

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

逍遥Fau

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值