LightOJ1236の LCMを形成するペア
タグ
序文
- 私のブログのCSDNと庭園がようこそ、同期されdanzh-ブログパーク〜
質問の簡潔な意味
\ [\ sum_ {i = 1} ^ n個の\ sum_ {J = 1} ^ N [LCM(I、J)= N] \]
考え
- Nの素因数分解の\(P_1 ^ N = {C1} ^ {P_2 P_KのC_2} ... {C_K} ^ \) LCMによって定義され、それぞれについて、のIJを知っている\(P_I \) 、いるはず\(MAX(C_1と、C_2)= C \) 、従って、IJを加える必要がCであり、別のオプションは、これが\(2C + 1 \)が選択される方法を化合物。だから、答えはする必要があります\(\のProd(2c_i + 1)\) 。
- 答えは<= 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;
}