lightoj 1236 pairs of lcm

1236 - Pairs Forming LCM

PDF (English)

Statistics

Forum

Time Limit: 2 second(s)

Memory Limit: 32 MB

Find the result of the following code:

long long pairsFormLCM( int n ) {
    long long res = 0;
    for( int i = 1; i <= n; i++ )
        for( int j = i; j <= n; j++ )
           if( lcm(i, j) == n ) res++; // lcm means least common multiple
    return res;
}

A straight forward implementation of the code may time out.If you analyze the code, you will find that the code actually counts the numberof pairs (i, j) for which lcm(i, j) = n and (i ≤ j).

Input

Input starts with an integer T (≤ 200),denoting the number of test cases.

Each case starts with a line containing an integer n (1≤ n ≤ 10¹⁴).

Output

For each case, print the case number and the value returnedby the function 'pairsFormLCM(n)'.

Sample Input	Output for Sample Input
15 2 3 4 6 8 10 12 15 18 20 21 24 25 27 29	Case 1: 2 Case 2: 2 Case 3: 3 Case 4: 5 Case 5: 4 Case 6: 5 Case 7: 8 Case 8: 5 Case 9: 8 Case 10: 8 Case 11: 5 Case 12: 11 Case 13: 3 Case 14: 4 Case 15: 2

 

---

Problem Setter: Jane Alam Jan

题意：求1-n中不同的两个数的最小公倍数等于n的个数。

分析：遇到这种题，暴力是不可能的，只有考虑别的方法。考虑将n进行质因数分解。从质因数中找一些线索。

因数分解以后，n=p1^a1 * p2^a2 * p3^a3 * ......  如果两个数（a,b）的lcm为n，把n拆成p1^a1, p2^a2, p3^a3...若干部分，如果a与b质因数分解后与之对应的每一部分幂的最大值小于n的对应部分的幂，那么lcm肯定不是n，所以n的每一部分的幂等于a与b对应部分的幂的最大值。例如第一部分，如果a取了最大值a1，那么b可以取0-a1共a1+1个，如果b取了最大值a1，那么a可以取0-a1共a1+1个，因此就是每次质因数分解质因数的个数 cnt*2-1 个，因为有两次（a1,a1）组合，所有情况相乘，最后除以二（每个数重复了两次）就是答案。

#include<bitset>
#include<map>
#include<vector>
#include<cstdio>
#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<stack>
#include<queue>
#include<set>
#define inf 0x3f3f3f3f
#define mem(a,x) memset(a,x,sizeof(a))

using namespace std;

typedef long long ll;
typedef pair<int,int> pii;

inline ll in()
{
    ll res=0;char c;
    while((c=getchar())<'0' || c>'9');
    while(c>='0' && c<='9')res=res*10+c-'0',c=getchar();
    return res;
}
const int N=10000010;
bitset<N> vis;
int prime[680000];
//一千万以内的素数开68万
//一百万以内的素数开8万
//十万以内的素数开1万
//基本除以10就行
//一万以内的素数直接开一万就行- -
int p;
int main()
{
    int T=in(),ii=1;
    for(int i=2;i<N;i++)
    {
        if(!vis[i])
        {
            prime[p++]=i;
            for(ll j=1LL*i*i;j<N;j+=i) vis[j]=1; //i*i别溢出了，不要只定义j为ll
        }
    }
    while(T--)
    {
        ll n=in();
        ll ans=1;
        for(int i=0;i<p && 1LL*prime[i]*prime[i]<=n;i++)
        {
            int cnt=0;
            while(n%prime[i] == 0)
            {
                n/=prime[i];
                cnt++;
            }
            ans *= ((cnt+1)<<1)-1;
        }
        if(n>1) ans *= 3;
        printf("Case %d: %lld\n",ii++,(ans+1)>>1);
    }
    return 0;
}