【LightOJ 1236】 Pairs Forming LCM 【算术基本定理+排列组合】

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

给n求, 上述代码的结果。

分析:从算术基本定理来考虑。

代码

#include<bits/stdc++.h>
using namespace std;
#define LL long long

const int N = 1e7+11;
const int M = 1e6+11;
const int mod = 1e9+7;
const int inf = 0x3f3f3f3;

bool su[N+1]={1,1,0};int prm[M],sz=0;
void init(){
    for(int i=2;i<N;i++){
        if(!su[i]) prm[sz++]=i;

        for(int j=0;j<sz;j++){
            LL t=prm[j]*i;
            if(t>N) break;
            su[t]=1;
            if(i%prm[j]==0) break;
        }
    }
   // cout<<sz;

}

LL solve(LL n){
    if(n==1) return 1;
    LL ans=1; int cnt=0;
    for(int i=0;i<sz && prm[i]*1ll*prm[i]<=n;i++){
        if(n%prm[i]==0){
            int ge=0; cnt++;
            while(n%prm[i]==0) { ge++; n/=prm[i];}
            ans=ans*(2*ge+1);
        }
        if(n==1) break;
    }

    if(n!=1) ans=ans*(2*1+1);
    ans++;
    return ans/2;
}

int main(){
    init();
    int t;scanf("%d",&t);
    int cas=1;
    while(t--){
        LL n;scanf("%lld",&n);
        printf("Case %d: %lld\n",cas++,solve(n));
    }
return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值