求一个数的所有因数+质因数分解【数论】

先附上所有因数的求法:

我的做法:是今天误打误撞写出来的;

http://exam.upc.edu.cn/problem.php?id=5062

然后,我上网找居然没有人写一个高效一点的,我这个做法其实就是.

不一定要会比根号N快,但是

模拟求所有因子个数的做法:

大家知道为什么所有因子的个数为:

设P1,P2……Pn都是数的质因子,

设C1,C2……Cn是数的质因子的个数:

Ans=(C1+1)*(C2+1)*……*(Cn+1)

大家想知道模拟一下就知道为什么了。

比如120=2*2*2*3*5;

那么Ans=(3+1)*(1+1)*(1+1)

要是求它的因子那么就是:

   (1+2+4+8)*(1+3)*(1+5)

=(1+2+4+8)+(3+6+12+24)+(5+10+20+40)+(15+30+60+120)

其中所有因数为:

1 , 2 , 3 , 4 , 5 , 6 , 8 , 10 , 12  ,  15  ,  20  ,  24  ,  30  ,  40  ,  60  ,  120

因数之和就是

=(1+2+4+8)+(3+6+12+24)+(5+10+20+40)+(15+30+60+120)

大家琢磨琢磨为什么是这样。

#include
using namespace std;
const int N=5e6+10;
typedef long long ll;
ll qpow(ll a,ll b){
    ll ans=1;
    while(b){
        if(b&1){
            ans=ans*a;
        }
        a=a*a;
        b>>=1;
    }
    return ans;
}
/*
 (2*3*5*7*11)*(13*17*19*23*29)*(31*37*41*43*47)
=614,88978,25884,91410
    前15个质数相乘得到一个 6 e18的数字
*/
vectorv[20];
mapVis;
ll a[N],cnt,Type;
void dfs(ll x,int step){
    if(Vis[x]==0&&step==Type){
        Vis[x]=1;
        a[cnt++]=x;
        return ;
    }
    for(int i=0;i

 归纳总结一下:质因数分解更快的做法:


1589: 题目名称:分解质因数

时间限制: 1 Sec 内存限制: 128 MB

题目描述

求出区间[a,b]中所有整数的质因数分解。

 

输入

输入两个整数ab

 

输出

每行输出一个数的分解,形如k=a1*a2*a3...(a1<=a2<=a3...k也是从小到大的)(具体可看样例)

 

样例输入

3 10

样例输出

3=3

4=2*2

5=5

6=2*3

7=7

8=2*2*2

9=3*3

10=2*5

 

提示

先筛出所有素数,然后再分解。

数据规模和约定

2<=a<=b<=10000

#include
using namespace std;
const int N=100;
typedef long long ll;
ll a[N],cnt;
void solve(ll n){
    ll tn=n;
    cnt=0;
    for(ll i=2;i*i<=n;i++){
        if(n%i==0){
            int D=0;
            while(n%i==0){
                D++;
                n/=i;
                a[cnt++]=i;
            }
        }
    }
    if(n!=1){
        a[cnt++]=n;
    }
    sort(a,a+cnt);
    printf("%lld=",tn);
    for(int i=0;i

 

你可能感兴趣的:(数论)