Counting Divisors of a Number in [tutorial]

№	Пользователь	Рейтинг
1	tourist	4009
2	jiangly	3823
3	Benq	3738
4	Radewoosh	3633
5	jqdai0815	3620
6	orzdevinwang	3529
7	ecnerwala	3446
8	Um_nik	3396
9	ksun48	3390
10	gamegame	3386

№	Пользователь	Вклад
1	cry	167
2	Um_nik	163
3	maomao90	162
3	atcoder_official	162
5	adamant	159
6	-is-this-fft-	158
7	awoo	157
8	TheScrasse	154
9	Dominater069	153
9	nor	153

Link to PDF (Latex Formatted)

Topic : Counting Divisors of a Number

Pre Requisites : Basic Maths , $\text{[math]}$ Factorisation , Primality testing

Motivation Problem :

There are T test cases. Each test case contains a number N. For each test case , output the number of factors of N.

1 < = T < = 10

1 < = N < = 10¹⁸

Note : 1 and N are also treated as factors of the number N.

Solution :

Naive Solution — $\text{[math]}$ Factorisation

The most naive solution here is to do the $\text{[math]}$ factorisation. But as we can see, it will surely receive the Time Limit Exceeded verdict. You may try to compute the prime numbers till required range and loop over them , but that too exceeds the usual limit of 10⁸ operations. Some optimisations and heuristics may allow you to squeeze through your solution, but usually at the cost of a few TLE's.

Supposing you fit in the above description and cannot think of anything better than the naive solution, I would like to present to you a very simple algorithm which helps you count the number of divisors in $\text{[math]}$ , which would be useful for this kind of questions. This algorithm only gives us the count of factors, and not the factors itself.

Firstly, let's do a quick recap of how we obtain the $\text{[math]}$ algorithm for any number N :

We write N as product of two numbers P and Q.

$\text{[math]}$

Looping over all values of P gives us the count of factors of N.

Before you move forward to the next section, it would be useful if you try to come up with an algorithm that finds the count of factors in $\text{[math]}$ . As a hint, the way of thinking is same as that of getting to the $\text{[math]}$ algorithm.

Counting factors in $\text{[math]}$

Let's start off with some maths to reduce our $\text{[math]}$ factorisation to $\text{[math]}$ for counting factors :

We write N as product of three numbers P, Q and R.

$\text{[math]}$

We can loop over all prime numbers in range $\text{[math]}$ and try to reduce N to it's prime factorisation, which would help us count the number of factors of N.

To know how to calculate divisors using prime factorisation, click here.

We will split our number N into two numbers X and Y such that X * Y = N. Further, X contains only prime factors in range $\text{[math]}$ and Y deals with higher prime factors ( $\text{[math]}$ ). Thus, gcd(X , Y) = 1. Let the count of divisors of a number N be denoted by the function F(N). It is easy to prove that this function is multiplicative in nature, i.e., F(m * n) = F(m) * F(n), if gcd(M,N) = 1. So, if we can find F(X) and F(Y), we can also find F(X * Y) or F(N) which is the required quantity.

For finding F(X), we use the naive trial division to prime factorise X and calculate the number of factors. Once this is done, we have Y = N / X remaining to be factorised. This may look tough, but we can see that there are only three cases which will cover all possibilities of Y :

$\text{[math]}$ is a prime number : F(Y) = 2.
$\text{[math]}$ is square of a prime number : F(Y) = 3.
$\text{[math]}$ is product of two distinct prime numbers : F(Y) = 4.

We have only these three cases since there can be at max two prime factors of Y. If it would have had more than two prime factors, one of them would surely have been $\text{[math]}$ , and hence it would be included in X and not in Y.

So once we are done with finding F(X) and F(Y), we are also done with finding F(X * Y) or F(N).

Pseudo Code :

N = input()
primes = array containing primes till 10^6
ans = 1
for all p in primes :
            if p*p*p > N:
                  break
            count = 1
            while N divisible by p:
                  N = N/p
                  count = count + 1
            ans = ans * count
if N is prime:
            ans = ans * 2
else if N is square of a prime:
            ans = ans * 3
else if N != 1:
            ans = ans * 4

Checking for primality can be done quickly using Miller Rabin. Thus, the time complexity is $\text{[math]}$ for every test case and hence we can solve our problem efficiently.

At this point, you may think that in a similar way, we can reduce this to $\text{[math]}$ by handling some cases. I have not thought much on it, but the number of cases to be handled are high since after trial division N could be factorised into one, two or three primes. This is easy enough to code in contest environment, which is our prime objective.

This trick is not quite commonly known and people tend to make bugs in handling the three cases. A problem in regionals which uses this trick directly :

Problem F | Codeforces Gym

You can try also this technique on problems requiring $\text{[math]}$ factorisation for practice purposes.

Hope you found this useful! Please suggest more problems to be added as well as any edits, if required.

Happy Coding!

/* this function calculates (a*b)%c taking into account that a*b might overflow */ long long mulmod(long long a,long long b,long long c){ long long x = 0,y=a%c; while(b > 0){ if(b%2 == 1){ x = (x+y)%c; } y = (y*2)%c; b /= 2; } return x%c; }

#include <algorithm> #include <iostream> #include <cstdlib> #include <vector> #include <ctime> #include <cmath> using namespace std; typedef long long ll; vector<int> prime, lpf; void sieve(int n) { prime.assign(1, 2); lpf.assign(n + 1, 2); for (int x = 3; x <= n; x += 2) { if (lpf[x] == 2) prime.push_back(lpf[x] = x); for (int i = 0; i < prime.size() && prime[i] <= lpf[x] && prime[i] * x <= n; ++i) lpf[prime[i] * x] = prime[i]; } } ll addMOD(ll a, ll b, const ll MOD) { return (a + b) % MOD; } ll mulMOD(ll a, ll b, const ll MOD) { // if (double(a) * b <= 1e18) return (a * b) % MOD; ll res = 0; for (a %= MOD; b > 0; a <<= 1, b >>= 1) { if (a >= MOD) a -= MOD; if (b & 1) { res += a; if (res >= MOD) res -= MOD; } } return res; } ll powMOD(ll x, ll n, const ll MOD) { ll res = 1; for (; n > 0; n >>= 1) { if (n & 1) res = mulMOD(res, x, MOD); x = mulMOD(x, x, MOD); } return res; } bool mr_test(ll n, int a, ll d, int s) { ll x = powMOD(a, d, n); if (x == 1 || x == n - 1) return false; for (int r = 1; r < s; ++r) { x = mulMOD(x, x, n); if (x == n - 1) return false; } return true; } bool miller_rabin(ll n, int k = 5) { if (n < 4) return n == 2 || n == 3; if (n % 2 == 0 || n % 3 == 0) return false; int s = __builtin_ctz(n - 1); ll d = (n - 1) >> s; for (int it = 1; it <= 5; ++it) { int a = 2 + rand() % (n - 3); if (mr_test(n, a, d, s)) return false; } return true; } bool isPrime(ll n) { if (n < 2) return false; if (n < lpf.size()) return lpf[n] == n; return miller_rabin(n); } ll any_divisor_of(ll n) { if (n <= 1) return 1; if (n < lpf.size()) return lpf[n]; if (isPrime(n)) return n; ll d = n; for (ll c = 2; d == n; ++c) { ll x = 2, y = 2, i = 2, k = 2; while (true) { x = (mulMOD(x, x, n) + c); if (x >= n) x -= n; d = __gcd(abs(x - y), n); if (d > 1 && n % d == 0) break; if (i++ == k) y = x, k <<= 1; } } return d; } bool isSquarePrime(ll n) { ll t = round(sqrt(n)); return (t * t == n && isPrime(t)); } bool isCubePrime(ll n) { ll t = round(cbrt(n)); return (t * t * t == n && isPrime(t)); } bool isSemiPrime(ll n) { ll d = any_divisor_of(n); return isPrime(d) && isPrime(n / d); } bool isSphenicPrime(ll n) { ll p = any_divisor_of(n); n /= p; if (isPrime(n)) swap(p, n); ll q = any_divisor_of(n); n /= q; return (p != q) && (q != n) && (n != p) && isPrime(p) && isPrime(q) && isPrime(n); } ll d(ll n) { int t = sqrt(sqrt(n)) + 1; sieve(t); ll res = 1; for (int p : prime) { if (1LL * p * p * p * p > n) break; if (n % p == 0) { int f = 0; do ++f; while ((n /= p) % p == 0); res *= (f + 1); } } if (n == 1) return res; if (isPrime(n)) return res * 2; if (isSquarePrime(n)) return res * 3; if (isSemiPrime(n)) return res * 4; if (isCubePrime(n)) return res * 4; if (isSphenicPrime(n)) return res * 8; return res * 6; } int main() { srand(time(NULL)); ll n; cin >> n; cout << d(n); return 0; } /// N = X * Y /// K = N^(1/4) /// X <= K /// /// for a, b, c distint prime /// > Y = a /// /// > Y = a^2 /// > Y = a * b /// /// > Y = a^3 /// > Y = a^2 * b /// > Y = a * b * c

#include <algorithm> #include <iostream> #include <cstdlib> #include <vector> #include <ctime> #include <cmath> #include <map> using namespace std; typedef long long ll; vector<int> prime, lpf; void sieve(int n) { prime.assign(1, 2); lpf.assign(n + 1, 2); for (int x = 3; x <= n; x += 2) { if (lpf[x] == 2) prime.push_back(lpf[x] = x); for (int i = 0; i < prime.size() && prime[i] <= lpf[x] && prime[i] * x <= n; ++i) lpf[prime[i] * x] = prime[i]; } } ll addMOD(ll a, ll b, const ll MOD) { return (a + b) % MOD; } ll mulMOD(ll a, ll b, const ll MOD) { // if (double(a) * b <= 1e18) return (a * b) % MOD; ll res = 0; for (a %= MOD; b > 0; a <<= 1, b >>= 1) { if (a >= MOD) a -= MOD; if (b & 1) { res += a; if (res >= MOD) res -= MOD; } } return res; } ll powMOD(ll x, ll n, const ll MOD) { ll res = 1; for (; n > 0; n >>= 1) { if (n & 1) res = mulMOD(res, x, MOD); x = mulMOD(x, x, MOD); } return res; } bool mr_test(ll n, int a, ll d, int s) { ll x = powMOD(a, d, n); if (x == 1 || x == n - 1) return false; for (int r = 1; r < s; ++r) { x = mulMOD(x, x, n); if (x == n - 1) return false; } return true; } bool miller_rabin(ll n, int k = 5) { if (n < 4) return n == 2 || n == 3; if (n % 2 == 0 || n % 3 == 0) return false; int s = __builtin_ctz(n - 1); ll d = (n - 1) >> s; for (int it = 1; it <= 5; ++it) { int a = 2 + rand() % (n - 3); if (mr_test(n, a, d, s)) return false; } return true; } bool isPrime(ll n) { if (n < 2) return false; if (n < lpf.size()) return lpf[n] == n; return miller_rabin(n); } ll any_divisor_of(ll n) { if (n <= 1) return 1; if (n < lpf.size()) return lpf[n]; if (isPrime(n)) return n; ll d = n; for (ll c = 2; d == n; ++c) { ll x = 2, y = 2, i = 2, k = 2; while (true) { x = (mulMOD(x, x, n) + c); if (x >= n) x -= n; d = __gcd(abs(x - y), n); if (d > 1 && n % d == 0) break; if (i++ == k) y = x, k <<= 1; } } return d; } ll d(ll n) { if (n < 0) n = -n; if (n == 0) return 1LL << 63 - 1; if (n == 1) return 1; vector<ll> S; if (n != 1) S.push_back(n); map<int, int> M; vector<ll> factor; while (S.size()) { ll x = S.back(); S.pop_back(); while (x > 1) { ll d = any_divisor_of(x); x /= d; if (isPrime(d)) ++M[d]; else S.push_back(d); } } ll res = 1; for (const pair<int, int> &e : M) res *= (e.second + 1); return res; } int main() { srand(time(NULL)); ll n; cin >> n; cout << d(n); return 0; } /// N = X * Y /// K = N^(1/4) /// X <= K /// /// for a, b, c distint prime /// > Y = a /// /// > Y = a^2 /// > Y = a * b /// /// > Y = a^3 /// > Y = a^2 * b /// > Y = a * b * c

Комментарии (53)

Показать архивные | Написать комментарий?

BLANKRK

9 лет назад, # |

Good one :) (Y)

→ Ответить

klamathix

You can't really call this " $\text{[math]}$ factorization" though, you only get the exponents, not the factors.

himanshujaju

9 лет назад, # ^ |

← Rev. 2 →

Agreed, that's why I put the title as "Counting Divisors". For the sub heading, I couldn't think of anything short and simple, so I went ahead with that. What could it be changed to?

Edit : I edited the post so that it isn't misleading anymore!

Looks good now, good work :).

VastoLorde95

← Rev. 5 →

No wait I missed some cases, ignore :P

$\text{[math]}$ isnt right, take 101 * 101 * 97 as an example, you just loop till 40!

What is N here?

N = 101 * 101 * 97.

Ah, I see, we don't necessarily have a prime factor less than the 4th root of the number

Auto comment: topic has been updated by himanshujaju (previous revision, new revision, compare).

Batman

-13

To check if N is prime or not at the end of your algorithm you have to use O(sqrt(n)) algorithm...

+11

Read up on Miller Rabin, I mentioned that explicitly in the post.

Ok...tnx...I'll read about it

kushalsingh007

You mentioned the following part -- We have only these three cases since there can be at max two prime factors of Y. If it would have had more than two prime factors, one of them would surely have been , and hence it would be included in X and not in Y.

Could you please explain this part a bit in detail. I am still now sure about why there can be at max two prime factors for Y and not 3 or more.

Suppose N has 3 primes factors left after the trial division is done. So, N = p₁ * p₂ * p₃.

Now , p₁ , p₂ and p₃ all have to be greater than $\text{[math]}$ or else they would have been found in the trial division. Their product will be greater than N , which is not possible and hence we contradict what we started with.

Try a few examples in pen and paper and read this proof again in case you still don't understand. Don't hesitate to reply for any clarifications.

fofao_funk

What is a good way for checking whether a given number is a square?

I'm using this:

bool is_square(ll n) {
	ll sqr = sqrt(n);
	return sqr * sqr == n;
}

bhishma

← Rev. 8 →

-20

Kostroma

It checks whether n is a power of two, not a square.

-8

Sorry guys I accept my mistake, is there a way to remove the comment

vikky_codder

Someone please provide code for Miller-Rabin primality test. Thanks!

mitsuki

8 лет назад, # ^ |

This provides a nice explanation as well

iabhishek15

4 года назад, # ^ |

Explanation was good but the code was not fully correct, plz can someone provide the code.

← Rev. 3 →

I tried solving the problem using Miller-Rabin as you said. The thing is: when we are limited to ~ 64 bits, doing an exponentiation a^d (mod n) may give overflow. This happens because, if n is too large (like 10^18), the value of a^d will overflow before we can apply the modulus operator (applying it when it is less than n won't help). Is there any way to overcome this or do we need to work with bigger types?

I use this to avoid overflow (taken from topcoder tutorial on primality) :

Illidan

+10

2015:


#define LL long long


LL mulmod(LL a, LL b, LL c) {
  return (__int128)a * b % c
}

CountZero

__int128 is a non-standart GCC extension. For example, Microsoft has no plans to implement it in MSVC

masterwayne

can u provide links to some spoj problem that uses this technique?

Unfortunately, I know only one question yet. You can use this on any question where $\text{[math]}$ logic for counting factors also passes. Havnt been solving on SPOJ recently, so no idea.

NaiveNaive

+17

The Pollard-rho algorithm can give the factorization of n in O(n^1 / 4) time.

uttom

8 лет назад, # |

When Y consist of two prime or Y is a prime how can i check it. Thanks in Advance.

Arunnsit

suppose Y<=10^18.

iterate over all primes from 2 to 10^6 and keep dividing Y by them meanwhile keep a count of factors.

After this step either you have found the solution or the Y left contains prime factors > 10^6 .

NOTE : now notice that Y cant have more than two prime factors at this step because factors left are already greater than 10^6 .

case 1: Y is prime , using miller rabin test check it .
else case 2 : Y has two prime factors and both of them are >= 10^6 .

Jori

Am I allowed to translate it to Chinese with your name on it?

Sure! With/without name doesn't matter much :)

XR_Maala

please can you explain why F(y) ={ 2 , 3 , 4 } in the three cases !! I understand why we have only three cases but I can't understand who we get F(y) value .. thanks in advance

sam29

Can you please share your accepted code for the above Problem.

Hossam

Can someone elaborate on how we can check if Y is any of the three cases? I can't see this through!

Vicennial

Check if Y is a prime number by using Miller Rabin test.
If it fails, check if the sqrt of Y is prime by using the same test.
If that fails too, then the number must belong to the third case.

Protap_Ghose

Did u solve the above problem Vicennial?-

upd: solved by removing the sqrt() fun

sagarthecoder2

7 лет назад, # |

ok. suppose, we want to check nod of a prime number and this number is n>10^6; so,in your pseudo code, the number failed to go through in 1st while loop; after it comes [ if N is prime ] this condition. but we stored prime only up to 10^6. since the number can't go through while loop and the value of the number (n) cannot be changed, how I check the number is prime or not???? how it passes this [ if N is prime ] condition???

Example, if a number is 50761746622, what is it's output??

Aemon

7 лет назад, # ^ |

You don't use the sieve to check if the result is prime. You use Rabin-miller probabilistic primality testing to check if it is prime. Look up "Rabin-Miller" and go to wikipedia to read about it. It is a way to check with high probability that a number is prime or not.

Surana

If N is a prime number greater than 1e17 then Y=N. So, how to check whether Y is prime or Y is a square of prime or Y is product of two prime? (Initially it is not known that N is prime greater than 1e17)

amanbol

Fermat's primality test

drinkWater.

6 лет назад, # |

I have written the for this topic, You can check it here : http://codeforces.me/blog/entry/62105

Amoo_Safar

Pollard's rho algorithm can do the factorization in $\text{[math]}$ complexity

Snapper_001

10 месяцев назад, # ^ |

-18

implementation any?

prodipdatta7

This one is the same as described one..

WillAchieveTheImpossible

5 лет назад, # |

← Rev. 6 →

i submitted this code on this problem and i got WA do you know what the wrong ?
my code : 69619960

SPyofgame

3 года назад, # |

At this point, you may think that in a similar way, we can reduce this to $$$O(N^{\frac{1}{4}})$$$ by handling some cases.

I know this is necroposting but we can reduce the complexity downto $$$O(N^{\frac{1}{4}})$$$ or generalized to $$$O(N^{\frac{1}{k}})$$$ (though I am not sure about the time complexity)
Let $$$N = X \times Y$$$ for $$$1 \leq X \leq N^{\frac{1}{4}} < Y \leq N$$$
Yes, it is possible, for $$$A, B, C$$$ distint primes

1) $$$Y = A$$$ is a prime $$$\rightarrow F(Y) = 2$$$

2) $$$Y = A^2$$$ is a square of prime $$$\rightarrow F(Y) = 3$$$

3) $$$Y = A^3$$$ is a cube of prime $$$\rightarrow F(Y) = 4$$$

4) $$$Y = A \times B$$$ is a squarefree semiprime $$$\rightarrow F(Y) = 2 \times 2$$$

5) $$$Y = A \times B \times C$$$ is a sphenic number $$$\rightarrow F(Y) = 2 \times 2 \times 2$$$

6) $$$Y = A^2 \times B$$$ is a cubefree product of square of a prime and a prime $$$\rightarrow F(Y) = 3 \times 2$$$

Here are the steps to do it

First let make a sieve

/// Linear Sieve
vector<int> prime, lpf;
void sieve(int n)
{
    prime.assign(1, 2);
    lpf.assign(n + 1, 2);

    for (int x = 3; x <= n; x += 2)
    {
        if (lpf[x] == 2) prime.push_back(lpf[x] = x);
        for (int i = 0; i < prime.size() && prime[i] <= lpf[x] && prime[i] * x <= n; ++i)
            lpf[prime[i] * x] = prime[i];
    }
}

Then we extract X from N

/// divisor count of (n)
ll d(ll n)
{   
    int t = sqrt(sqrt(n)) + 1;
    sieve(t);

    ll res = 1;
    for (int p : prime)
    {
        if (1LL * p * p * p * p > n) break;
        if (n % p == 0)
        {
            int f = 0;
            do ++f; while ((n /= p) % p == 0);
            res *= (f + 1);
        }
    }

    .... /// extract Y
}

Let check if (n) is a prime

/// Calculate a * b (mod MOD)
ll mulMOD(ll a, ll b, const ll MOD) /// O(log n)
{
    // if (double(a) * b <= 1e18) return (a * b) % MOD; /// for small number | should not to call when there are many large multiplications

    ll res = 0;
    for (a %= MOD; b > 0; a <<= 1, b >>= 1)
    {
        if (a >= MOD) a -= MOD;
        if (b & 1) 
        {
            res += a;
            if (res >= MOD) res -= MOD;
        }
    }

    return res;
}

/// Calculate x ^ n (mod MOD)
ll powMOD(ll x, ll n, const ll MOD) /// O(log n) * O(log n)
{
    ll res = 1;
    for (; n > 0; n >>= 1)
    {
        if (n & 1) res = mulMOD(res, x, MOD);
        x = mulMOD(x, x, MOD);
    }

    return res;
}

/// miller rabin primality test for specific value | Accuracy: 75%+
bool mr_test(ll n, int a, ll d, int s) /// O(s * log(n) * log(n))
{
    ll x = powMOD(a, d, n);
    if (x == 1 || x == n - 1) return false;

    for (int r = 1; r < s; ++r)
    {
        x = mulMOD(x, x, n);
        if (x == n - 1) return false;
    }

    return true;
}

/// miller rabin function to generate primality tests
bool miller_rabin(ll n, int k = 5) /// O(k * log n * log n * log n)
{
    if (n < 4) return n == 2 || n == 3;
    if (n % 2 == 0 || n % 3 == 0) return false;

    int s = 0;
    ll d = n - 1;
    while (!(d & 1)) d >>= 1, ++s;

    for (int it = 1; it <= k; ++it)
    {
        int a = 2 + rand() % (n - 3);
        if (mr_test(n, a, d, s)) return false;
    }

    return true;
}

/// check if (n) is a prime
bool isPrime(ll n) /// O(log^3) or O(1) with precalculation
{
    if (n < 2) return false;    
    if (n < lpf.size()) return lpf[n] == n; /// small number - O(1)
    return miller_rabin(n);                 ///  big  number - O(log^3)
}

Let check if (n) is a square of a prime

/// check if (n) is a square of a prime
bool isSquarePrime(ll n) /// O(log^3) or O(1) with precalculation
{
    ll t = round(sqrt(n));
    return (t * t == n && isPrime(t));
}

Let check if (n) is a cube of a prime

/// check if (n) is a cube of a prime
bool isCubePrime(ll n) /// O(log^3) or O(1) with precalculation
{
    ll t = round(cbrt(n));
    return (t * t * t == n && isPrime(t));
}

We can find a divisor of integer by pollard rho and miller rabin

/// find any divisor of (n), better to not be (-n), (-1), (1), (n)
ll any_divisor_of(ll n) /// ˜O(s(N)1/2) or O(1) with precalculation
{
    if (n <= 1) return 1;
    if (isPrime(n)) return n;
    if (n < lpf.size()) return lpf[n];

    ll d = n;
    for (ll c = 2; d == n; ++c)
    {    
        ll x = 2, y = 2, i = 2, k = 2;
        while (true)
        {
            x = (mulMOD(x, x, n) + c);
            if (x >= n)	x -= n;
            d = __gcd(abs(x - y), n);
            if (d > 1 && n % d == 0) break;
            if (i++ == k) y = x, k <<= 1;
        }
    }

    return d;
}

Let check if (n) is a semiprime

/// check if (n) is a semiprime
/// both 2 divisors are primes
bool isSemiPrime(ll n) /// ˜O(s(N)1/2) or O(1) with precalculation
{
    ll d = any_divisor_of(n);
    return isPrime(d) && isPrime(n / d);
}

Let check if (n) is a sphenic number

/// check if (n) is a sphenic prime
/// all 3 divisors are distinct primes
bool isSphenicPrime(ll n) /// ˜O(s(N)1/2) or O(1) with precalculation
{
    ll p = any_divisor_of(n);       n /= p;
    if (isPrime(n)) swap(p, n);

    ll q = any_divisor_of(n);       n /= q;
    return (p != q) && (q != n) && (n != p) && isPrime(p) && isPrime(q) && isPrime(n);
}

Then we extract Y from N

/// divisor count of (n)
ll d(ll n)
{   
    ... /// extract X

    if (n == 1) return res;
    if (isPrime(n)) return res * 2;
    if (isSquarePrime(n)) return res * 3;
    if (isSemiPrime(n)) return res * 4;
    if (isCubePrime(n)) return res * 4;
    if (isSphenicPrime(n)) return res * 8;
    return res * 6;
}

Finally the main function



int main()
{
    srand(time(NULL));

    ll n;
    cin >> n;

    cout << d(n);
    return 0;
}

You can submit at GCPC15 — F. It is as fast as $$$O(N^{\frac{1}{3}})$$$ solution in this case

Full Code ~ Asleast asfast as O(N^(1/3)) solution

Simpler code but run slower

8 месяцев назад, # ^ |

What is the complexity of the the function (anyDivisor of n) , Also what you do in it

← Rev. 4 →

It is $$$Õ\left(\sqrt{_\sqrt{n}^{\text{ }}}\right) = Õ\left(n^{1/4}\right)$$$ for which $$$Õ(n) = O\left(n \text{ polylog } n\right)$$$

In my implementation, you can optimize them in certain ways to optimize further for general cases. (For examples: deterministic miller rabin, $$$O(1)$$$ modulo multiplication, reduce amount of modulo operations).

I was explaining that if there is an algorithm better than pollard rho exists, let it be $$$O(f(n))$$$ then we can count divisors by making $$$N = X_1 \times X_2 \times ... \times X_k$$$, for $$$1 \leq X_1 \leq N^{1/k} \leq \frac{N}{X_1} \leq N$$$ then check for each case whether $$$X_i = X_j$$$ are the same prime or not (each yelling a different $$$F(Y)$$$ value mentioned earlier in my comment)

Ohh it is pollard rho I never have its implementation and now i have

Thanks

Bashca

21 месяц назад, # |

We don't really need to use Miller-Rabin test, because there are no Carmichael numbers with two prime factors. we can use Fermat test instead.

_lbw_

Why not use Pollard's rho algorithm?

Блог пользователя himanshujaju

Counting factors in

Counting factors in $\text{[math]}$