L1-006 连续因子

L1-006 连续因子 (20 分)

一个正整数 N 的因子中可能存在若干连续的数字。例如 630 可以分解为 3×5×6×7,其中 5、6、7 就是 3 个连续的数字。给定任一正整数 N,要求编写程序求出最长连续因子的个数,并输出最小的连续因子序列。

输入格式:

输入在一行中给出一个正整数 N(1<N<231​​)。

输出格式:

首先在第 1 行输出最长连续因子的个数;然后在第 2 行中按 因子1*因子2*……*因子k 的格式输出最小的连续因子序列,其中因子按递增顺序输出,1 不算在内。

输入样例:

630

输出样例:

3
5*6*7

此题一般有两种解法:

第一种:

从1乘到N,每乘一次判定当前乘积prd是否为N的因子,即N%prd是否为0,若为0,比较上一次乘积因子序列的长度,若大于,则记录。 不断更新乘积因子序列长度,直到最后一个因子为N。当然,最需要注意的是最外层的循环变量i判定条件必须为i<=sqrt(n),如果是i*i<=n,则最后一个测试点会不过。本人猜测可能是平台的问题。虽然看上去很暴力,但此算法最后一个测试点耗时仅5ms,目前应该是最快的算法,而且比较好理解。

 

#include<iostream>
#include<cmath>
using namespace std;
typedef long long ll;
int main()
{
    ll n;
    cin>>n;
    ll prd=0;//定义乘积 
    int start=0,len=0;//定义最终得到序列开始的因子,序列的长度 
    for(int i=2;i<=sqrt(n);i++)//i从2到根号n 
    {
        prd=1;
        for(int j=i;prd*j<=n;j++)//从j开始一直乘到N为止,每次乘积判定是否小于等于N,若超过N,则结束循环 
        {
            prd*=j;//乘积迭代 
            if(n%prd==0&&j-i+1>len)//如果当前乘积为N的乘积因子且长度大于上一次 
            {//更新序列起始因子和长度 
                start=i;
                len=j-i+1;
            }
        }
    }
    if(start==0)//若起始因子为0,说明N为质数,因为质数=1*本身,而循环最多能表示1*本身的根号 
    {
        start=n;
        len=1;
    }
    cout<<len<<'\n'<<start;
    for(int i=start+1;i<start+len;i++)//输出,如果因子只有一个只输出一个 
    cout<<'*'<<i;
    return 0;
}

第二种:

思想:从大到小乘积,一次性乘完,再判定,这样当命中所求乘积时,因为是从大到小乘积,故此时的序列长度就是最大的,然后及时跳出循环,大大减少循环次数,乘积完一轮再减少乘积长度。根据题目N的上界,13!刚好可以覆盖这个上界。于是除去1的乘积长度不超过12,虽然没有了频繁更新序列长度的操作,但是其最坏复杂度则耗时较长,一方面是都需要乘积后判定,另一方面是还是有很多不必要的乘积操作,但是就此算法而言不可避免。因为拿N=36来说,根号N=6,于是按照此算法的思想,先从2乘到13,再从3乘到14...6乘到17,然后减少乘积长度,2乘到12...6乘到16,不断减少乘积长度至1。不难发现,第一次的乘积操作2乘到13已经覆盖N的上界了,再从3乘到14也大于第一次的乘积结果,每一次的乘积的中途可能就已经满足等于N或者是N的乘积因子的条件,即假设M满足条件,则3x4x5....xMx(M+1)x(M+2)x.....14。而此算法在当前一轮乘积中,只得到最终乘到14的结果,无法得到中间的M,如果要得到M,必须在每轮乘完后的序列减少操作不断执行后才可能得到。即必须乘完3到14,4到15,5到16,6到17,减少一次长度,再乘一轮,减少一次长度,总共需要进行14-M轮乘积才能得到M。这其中的操作也大大降低了算法执行的效率。需要注意的是,乘积结果和N都必须是long long 类型,而且循环变量的上界是小于等于根号N,如果是小于会有一个点不过。这是因为根号函数返回的是double,转换为int只能向下取整。此算法最后一个测试点耗时15ms,显然比第一种慢很多。

#include<iostream>
#include<cmath>
using namespace std;
typedef long long ll;
int main()
{
    ll n;
    cin>>n;
    ll prd;
    int rootn=sqrt(n);//得到根号N 
    int flag=0,start,len;//定义是否为乘积因子的标识,乘积序列开始的因子,序列长度 
    for(len=12;len>=1;len--)//序列最长为12,递减到1 
    {
        for(start=2;start<=rootn;start++)//从当前一轮乘积因子的上界从2开始到根号N,注意一定是小于等于,否则有一个点会不过 
        {
            prd=1;
            for(int i=start;i<start+len;i++)//从当前乘积因子开始乘积,乘积len个长度 
            prd*=i;
            if(n%prd==0)//如果找到乘积因子 
            {
                flag=1;
                break;//标识,及时退出 
            }
        }
        if(flag)
        break;
    }
    if(!flag)//如果未标识为1,说明是质数 
    cout<<1<<endl<<n;
    else
    {
        cout<<len<<endl<<start;
        for(int i=start+1;i<start+len;i++)//输出,如果只有1个输出一个 
        cout<<'*'<<i;
    }
    return 0;
}

 

 

转载于:https://www.cnblogs.com/yfr2zaz/p/10409076.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值