当前位置: 代码迷 >> 综合 >> LightOJ1236 - Pairs Forming LCM (LCM·唯一分解)
  详细解决方案

LightOJ1236 - Pairs Forming LCM (LCM·唯一分解)

热度:10   发布时间:2023-12-06 19:47:19.0

题意 :给你一个数n  求满足lcm(a, b) == n, a <= b 的 (a,b) 的个数

唯一分解定理内容:对于任意一个大于1的数都可以唯一分解为若干个素数的乘积,即n=a1^b1*a2^b2*......an^bn;

容易知道 n 是由a和b的所有素因子构成的,n的中的素因子指数等于a,b中相同素因子指数较大的指数。

先将n分解为素数指数积的形式  n = π(pi^ei)    那么对于每个素因子pi  pi在a,b中的指数ai, bi 至少有一个等于pi, 另一个小于等于pi

先不考虑a, b的大小  对于每个素因子pi

1. 在a中的指数 ai == ei   那么 pi 在 b 中的指数可取 [0, ei] 中的所有数  有 ei + 1 种情况

2. 在a中的指数 ai < ei  即 ai 在 [0, ei) 中   那么 pi 在 b 中的指数只能取 ei  有 ei 种情况

那么对与每个素因子都有 2*ei + 1种情况  也就是满足条件的 (a, b) 有 π(2*ei + 1)个  考虑大小时 除了 (n, n) 所有的情况都出现了两次  那么满足a<=b的有 (π(2*ei + 1)) / 2 + 1  个

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 1e7+5;
int prim[N/10], m;
bool vis[N];
void initPrime()
{m = 0;for(ll i = 2; i < N; i++){if(!vis[i])prim[m++] = i;for(ll j = i*i; j < N; j += i)vis[j] = true;}
}
int main()
{initPrime();int T;ll n, ans, c;scanf("%d", &T);for(int kase = 1; kase <= T; kase++){scanf("%lld", &n);ans = 1;for(int i = 0; i < m; i++){if(ll(prim[i])*prim[i] > n)break;c = 0;while(n % prim[i] == 0){n /= prim[i];c++;}if(c)ans *= c * 2 + 1;}if(n > 1)ans *= 3;printf("Case %d: %lld\n", kase, ans/2+1);}return 0;
}