POJ 3421(X-factor Chains)素数

最新推荐文章于 2022-01-15 21:20:59 发布

Eknight123

最新推荐文章于 2022-01-15 21:20:59 发布

阅读量341

点赞数

分类专栏： ACM 文章标签： prime

本文链接：https://blog.csdn.net/Eknight123/article/details/81260819

版权

该博客探讨了POJ 3421问题，涉及寻找序列1, a0, a1, ..., am (am = X)，其中满足任意i时，ai % ai-1 == 0。最大m值等于X的质因子个数，而序列个数为所有质因子的全排列数。举例说明，当X=4时，质因子为{2,2}，m为2，全排列只有1种。" 107653372,9949672,AHB总线Slave设备的hready_in与hready_out交互解析,"['IC设计', '总线协议']

摘要由CSDN通过智能技术生成

题目链接：https://vjudge.net/problem/POJ-3421

题意：输入1个X，求解序列1,a0,a1,a2......am(am=X)，其中满足任意i(ai%ai-1==0)，求解最大的m和满足的序列个数

求解X的质因子个数即为m，再算出所有质因子的全排列个数。

n=4时，质因子{2,2}，m为2，全排列只有一种

#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<cmath>
#include<vector>
#include<cstring>
#define LL long long
using namespace std;
#define N 1100000
int a[1100000];
void is_sushu()
{
    memset(a,0,sizeof(a));//对数组a进行初始化为0，不是素数的标记为1，剩下为0的就是素数了
    a[1]=1;//1既不是素数也不是合数，先标记为0
    for(int i=2; i<=sqrt(N); i++)
    {
        if(a[i]==0)//如果i是素数
        {
            for(int j=2; j*i<=N; j++) //循环标记的范围是i*j<N
            {
                a[i*j]=1;//如果i是素数，那么i*j肯定不是素数