LightOJ1236 Pairs Forming LCM 水題？。。

---

LightOJ1236 Pairs Forming LCM

---

標籤

---

前言

• 我的csdn和博客園是同步的，歡迎來訪danzh-博客園~

---

簡明題意

$\sum_{i=1}^n\sum_{j=1}^n[lcm(i,j)=n]$

---

思路

• 對n質因數分解，$n=p_1^{c1}p_2^{c_2}...p_k^{c_k}$，由lcm的定義可以知道，ij中，對於每一個$p_i$，都應該有$max(c_1,c_2)=c$，因此，需要讓ij中的一個爲c，另一個任選，這樣就有$2c+1$種選法。所以答案就應該是$\prod (2c_i+1)$。
• 答案需要a<=b，所以答案需要/2（且要向上取整）

---

注意事項

• 無

---

總結

• 無

---

AC代碼

#include<cstdio>

const int maxn = 1e7 + 10;

bool no_prime[maxn];
int prime[(int)7e5];
int shai(int n)
{
	int cnt = 0;
	no_prime[1] = 1;

	for (int i = 2; i <= n; i++)
	{
		if (!no_prime[i])
			prime[++cnt] = i;

		for (int j = 1; j <= cnt && prime[j] * i <= n; j++)
		{
			no_prime[prime[j] * i] = 1;
			if (i % prime[j] == 0) break;
		}
	}
	return cnt;
}

void solve()
{
	int cnt = shai(maxn - 10);

	int t;
	scanf("%d", &t);
	for (int i = 1; i <= t; i++)
	{
		long long n, r;
		scanf("%lld", &n);
		r = n;

		long long ans = 1;
		for (int i = 1; i <= cnt && 1ll * prime[i] * prime[i] <= r && n != 1; i++)
		{
			int cnt = 0;
			while (n % prime[i] == 0)
				n /= prime[i], cnt++;
			ans *= (2 * cnt + 1);
		}
		if (n != 1)
			ans *= 3;

		printf("Case %d: %lld\n", i, (ans + 1) / 2);
	}
}

int main()
{
	freopen("Testin.txt", "r", stdin);
	solve();
	return 0;
}