Pairs Forming LCM 素数筛和素因子分解

本文详细解析了一个关于最小公倍数(LCM)的算法题目,通过素数筛和素因子分解的方法实现了题目的解答。介绍了两种素数筛算法,并强调了选择合适算法的重要性以避免超时。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

Pairs Forming LCM  :http://acm.hust.edu.cn/vjudge/contest/view.action?cid=109329#problem/B    密码:nefu

一个题目分了好几波才把它完全搞明白,公式还是很多的,直接偷懒上手写了。








这次写的也是比较啰嗦,但却是自己的思路,先来两个素数筛的代码,再贴上素因子分解的代码,最后才是本题的AC代码哟~

素数筛一:

void getprime1()
{
    memset(prime,0,sizeof(prime));
    for(int i=2;i<MAXN1;i++)
    {
        if(!prime[i])
        prime[++prime[0]]=i;
        for(int j=1;j<=prime[0]&&prime[j]<=MAXN1/i;j++)
        {
            prime[prime[j]*i]=1;
            if(i%prime[j]==0)
            break;
        }
    }
}

素数筛二:

void getprime()
{
    memset(isprime,true,sizeof(isprime));
    int nprime=0;
    isprime[1]=0;
    for(long long i=2;i<MAXN;i++)
    {
        if(isprime[i])
        {
            prime[++nprime]=i;  //从下标为1开始使用
            for(long long j=i*i;j<MAXN;j+=i)
            {
                isprime[j]=false;
            }
        }
    }
}

需要注意的是,方法二比方法一更优一些,此题必须用方法二,否则会超时。


素因子分解:

long long factor[10010];
int cnt;
int getfactors(long long x)
{
    cnt=0;
    long long tmp=x;
    for(int i=1;prime[i]<=tmp/prime[i];i++)  //从下标1开始使用素数数组
    {
        factor[cnt]=0; //从下标0开始存储每一个素数的指数
        if(tmp%prime[i]==0)
        {
            while(tmp%prime[i]==0)
            {
                factor[cnt]++;
                tmp/=prime[i];
            }
            cnt++;
        }
    }
    if(tmp!=1)
    {
        factor[cnt++]=1;
    }
    return cnt;
}

最终代码实现:

#include <iostream>
#include <cstdio>
#include <cstring>

using namespace std;
const int MAXN=10000100;
const int MAXN1=1000000;
bool isprime[MAXN];
int prime[MAXN/10];

void getprime()
{
    memset(isprime,true,sizeof(isprime));
    int nprime=0;
    isprime[1]=0;
    for(long long i=2;i<MAXN;i++)
    {
        if(isprime[i])
        {
            prime[++nprime]=i;  //从下标为1开始使用
            for(long long j=i*i;j<MAXN;j+=i)
            {
                isprime[j]=false;
            }
        }
    }
}

void getprime1()
{
    memset(prime,0,sizeof(prime));
    for(int i=2;i<MAXN1;i++)
    {
        if(!prime[i])
        prime[++prime[0]]=i;
        for(int j=1;j<=prime[0]&&prime[j]<=MAXN1/i;j++)
        {
            prime[prime[j]*i]=1;
            if(i%prime[j]==0)
            break;
        }
    }
}

long long factor[10010];
int cnt;
int getfactors(long long x)
{
    cnt=0;
    long long tmp=x;
    for(int i=1;prime[i]<=tmp/prime[i];i++)  //从下标1开始使用素数数组
    {
        factor[cnt]=0; //从下标0开始存储每一个素数的指数
        if(tmp%prime[i]==0)
        {
            while(tmp%prime[i]==0)
            {
                factor[cnt]++;
                tmp/=prime[i];
            }
            cnt++;
        }
    }
    if(tmp!=1)
    {
        factor[cnt++]=1;
    }
    return cnt;
}

int main()
{
    int t,casenum;
    casenum=0;
    long long n;
    long long sum;
    getprime();
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lld",&n);
        sum=1;
        int k=getfactors(n);
        for(int i=0;i<k;i++)
        {
            sum*=(2*factor[i]+1);
        }
        sum=(sum+1)/2;
        printf("Case %d: %lld\n",++casenum,sum);
    }
    return 0;
}






### Lua 中 `ipairs` `pairs` 的区别及用法 #### 使用场景差异 `ipairs` 主要用于遍历数组类型的表,即索引为连续正整数的表。这种情况下,`ipairs` 只会访问那些具有连续整数键的元素[^1]。 ```lua local array = { "apple", "banana", "orange" } for index, value in ipairs(array) do print(index, value) end ``` 相比之下,`pairs` 更加通用,可以用来遍历任何类型的表,无论其键是什么形式(字符串、数字或其他)。这使得 `pairs` 能够访问到所有的键值对,而不仅仅是带有整数键的部分。 ```lua local mixedTable = { fruit = "apple", [2] = "banana", ["third"] = "orange" } for key, value in pairs(mixedTable) do print(key, value) end ``` #### 返回迭代器的不同 `ipairs` 提供了一个简单的迭代接口来按顺序获取数值型键及其对应的值;当遇到第一个非数值型键或超出范围时停止迭代。因此,在处理可能含有间隙或者混合类型键的表时应谨慎使用 `ipairs`。 另一方面,`pairs` 则返回一个完整的迭代函数,它能够依次取出表中的每一个键值对直到结束。这意味着即使存在不规则分布的数据结构也能被完全枚举出来。 #### 性能考量 由于 `ipairs` 需要在每次调用时检查当前项是否满足条件(如是否为有效的整数下标),所以在某些特定条件下可能会稍微慢于直接通过 `pairs` 进行无差别扫描的方式。然而对于大多数实际应用场景而言,两者之间的性能差距通常是可以忽略不计的。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值