liyunlonghere

你必须非常努力,才能看起来毫不费力。

lightoj 1236 pairs of lcm



1236 - Pairs Forming LCM
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 ≤ 1014).

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;
}






阅读更多
版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/liyunlong41/article/details/49932157
文章标签: 数论
个人分类: 数论相关
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

不良信息举报

lightoj 1236 pairs of lcm

最多只允许输入30个字

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭