POJ 3421 X-factor Chains(dp or 素因子分解)

X-factor Chains
Time Limit: 1000MS Memory Limit: 65536K
Total Submissions:7875 Accepted: 2505

Description

Given a positive integer X, an X-factor chain of length m is a sequence of integers,

1 = X0X1X2, …, Xm = X

satisfying

Xi < Xi+1 and Xi | Xi+1 where a | b means a perfectly divides into b.

Now we are interested in the maximum length of X-factor chains and the number of chains of such length.

Input

The input consists of several test cases. Each contains a positive integer X (X ≤ 220).

Output

For each test case, output the maximum length and the number of such X-factors chains.

Sample Input

2
3
4
10
100

Sample Output

1 1
1 1
2 1
2 2
4 6
dp:
1000000的因子不足50个,直接dp便可求得。
代码:


#include<cstdio>
#include<algorithm>
#include<string.h>
using namespace std;
int dp[3003];
long long sum[3003];
int main()
{
    int x;
    while(~scanf("%d",&x))
    {
        int y[3003],r=0;
        for(int i=2;i*i<=x;i++)
        {
            if(x%i==0)
            {
                if(i*i==x)
                    y[r++]=i;
                else y[r++]=i,y[r++]=x/i;
            }
        }


        y[r++]=x;
        sort(y,y+r);
        for(int i=0;i<=3000;i++)
            dp[i]=sum[i]=1;
        int ma=0;
        for(int i=0;i<r;i++)
        {
            for(int j=0;j<i;j++)
            {
                if(y[i]%y[j]==0)
                {
                    if(dp[j]+1>dp[i])
                        dp[i]=dp[j]+1,sum[i]=sum[j];
                    else if(dp[j]+1==dp[i])
                    {
                        sum[i]=sum[i]+sum[j];
                    }
                }
            }
            ma=max(ma,dp[i]);
        }
        long long s=0;
        for(int i=0;i<r;i++)
            if(dp[i]==ma)
            s+=sum[i];
        printf("%d %lld\n",ma,s);
    }
    return 0;
}
素因子分解:
那100来说,先让他除以最小素数直到不能增除,100/2=50;
50/2=25。不能整出了,100=2^2*25;
25除以下一个素数,25/3不能整除,25/5=5,5/5=1;除到1结束。
所以100=2^2*5^2;最长因子便是指数相加2+2=4;
因为因子要保证后一项整除前一项,
这四个因子便是
1、2 2*2 2*2*5 2*2*5*5
2、2 2*5 2*5*2 2*5*2*5
3、2 2*5 2*5*5 2*5*5*2
4、5 5*5 5*5*2 5*5*2*2
5、5 5*2 5*2*2 5*2*2*5
6、5 5*2 5*2*5 5*2*5*2
可以看出如果有长度为4 便是四个因子的随机组合,有4!种
但是其中有重复,2有2个,5有2个,去掉重复的情况便是
4!/(2!+2!)
代码:

#include<cstdio>
#include<string.h>
bool mp[1100005];
long long prim[300002],K[33];
int L;
using namespace std;
void sushu()
{
    L=0;
    memset(mp,1,sizeof(mp));
    for(int i=2;i<=1100000;i++)
    {
        if(mp[i])
        {
            prim[++L]=i;
            for(int j=i+i;j<=1100000;j=j+i)
                mp[j]=0;
        }
    }
    K[1]=1;
    for(long long i=2;i<=30;i++)
        K[i]=i*K[i-1];
}
int main()
{
    sushu();
    long long n;
    while(~scanf("%lld",&n))
    {
        int i,j;bool bb=0;
        long long a=0,s=0,r=1;
        for(i=1;i<=L;i++)
        {
            while(n%prim[i]==0)
            {
                bb=1;a++;
                n=n/prim[i];
            }
            if(bb)bb=0,s+=a,r=r*K[a],a=0;
            if(n==1)break;
        }
        printf("%lld %lld\n",s,K[s]/r);
    }
    return 0;
}


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
题目描述 给出一个$n\times m$的矩阵,每个位置上有一个非负整数,代表这个位置的海拔高度。一开始时,有一个人站在其中一个位置上。这个人可以向上、下、左、右四个方向移动,但是只能移动到海拔高度比当前位置低或者相等的位置上。一次移动只能移动一个单位长度。定义一个位置为“山顶”,当且仅当从这个位置开始移动,可以一直走到海拔高度比它低的位置上。请问,这个矩阵中最多有多少个“山顶”? 输入格式 第一行两个整数,分别表示$n$和$m$。 接下来$n$行,每行$m$个整数,表示整个矩阵。 输出格式 输出一个整数,表示最多有多少个“山顶”。 样例输入 4 4 3 2 1 4 2 3 4 3 5 6 7 8 4 5 6 7 样例输出 5 算法1 (递归dp) $O(nm)$ 对于这道题,我们可以使用递归DP来解决,用$f(i,j)$表示以$(i,j)$为起点的路径最大长度,那么最后的答案就是所有$f(i,j)$中的最大值。 状态转移方程如下: $$ f(i,j)=\max f(x,y)+1(x,y)是(i,j)的下一个满足条件的位置 $$ 注意:这里的状态转移方程中的$x,y$是在枚举四个方向时得到的下一个位置,即: - 向上:$(i-1,j)$ - 向下:$(i+1,j)$ - 向左:$(i,j-1)$ - 向右:$(i,j+1)$ 实现过程中需要注意以下几点: - 每个点都需要搜一遍,因此需要用双重for循环来枚举每个起点; - 对于已经搜索过的点,需要用一个数组$vis$来记录,防止重复搜索; - 在进行状态转移时,需要判断移动后的点是否满足条件。 时间复杂度 状态数为$O(nm)$,每个状态转移的时间复杂度为$O(1)$,因此总时间复杂度为$O(nm)$。 参考文献 C++ 代码 算法2 (动态规划) $O(nm)$ 动态规划的思路与递归DP类似,只不过转移方程和实现方式有所不同。 状态转移方程如下: $$ f(i,j)=\max f(x,y)+1(x,y)是(i,j)的下一个满足条件的位置 $$ 注意:这里的状态转移方程中的$x,y$是在枚举四个方向时得到的下一个位置,即: - 向上:$(i-1,j)$ - 向下:$(i+1,j)$ - 向左:$(i,j-1)$ - 向右:$(i,j+1)$ 实现过程中需要注意以下几点: - 每个点都需要搜一遍,因此需要用双重for循环来枚举每个起点; - 对于已经搜索过的点,需要用一个数组$vis$来记录,防止重复搜索; - 在进行状态转移时,需要判断移动后的点是否满足条件。 时间复杂度 状态数为$O(nm)$,每个状态转移的时间复杂度为$O(1)$,因此总时间复杂度为$O(nm)$。 参考文献 C++ 代码

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值