Superhero’s Job

It’s tough to be a superhero. And it’s twice as tough to resist the supervillain who is cool at math. Suppose that you’re an ordinary Batman in an ordinary city of Gotham. Your enemy Joker mined the building of the city administration and you only have several minutes to neutralize the charge. To do that you should enter the cancel code on the bomb control panel.

However, that mad man decided to give you a hint. This morning the mayor found a playing card under his pillow. There was a line written on the card:

The bomb has a note saying ” J(x) = A“, where A is some positive integer. You suspect that the cancel code is some integer x that meets the equation J(x) = A. Now in order to decide whether you should neutralize the bomb or run for your life, you’ve got to count how many distinct positive integers x meet this equation.

Input

The single line of the input contains a single integer A (1 ≤ A ≤ 1012).

Output

Print the number of solutions of the equation J(x) = A.

Examples

input

3

output

1

input

24

output

3

Note

Record x|n means that number n divides number x.

 is defined as the largest positive integer that divides both a and b.

In the first sample test the only suitable value of x is 2. Then J(2) = 1 + 2.

In the second sample test the following values of x match:

  • x = 14, J(14) = 1 + 2 + 7 + 14 = 24
  • x = 15, J(15) = 1 + 3 + 5 + 15 = 24
  • x = 23, J(23) = 1 + 23 = 24

Solution:

#include <bits/stdc++.h>

using namespace std;

int const P = 1234567;
int primes[P], np[P];

int const N = 1234567;
long long d[N];
pair<long long, long long> zd[N];
long long dp[N], ndp[N];

int main() {
long long a;
cin >> a;
int cn = 0;
for (long long i = 1; i * i <= a; i++) {
        if (a % i != 0) continue;
        d[cn++] = i;
        if (i * i != a) d[cn++] = a / i;
    }
    for (int i = 2; i * i < P; i++) {
        if (!np[i]) {
            for (int j = i * i; j < P; j += i) np[j] = 1;
        }
    }
    int pc = 0;
    for (int i = 2; i < P; i++) {
        if (!np[i])
            primes[pc++] = i;
    }
    sort(d, d + cn);
    fill(dp, dp + cn, 0);
    dp[0] = 1;
    int edgecount = 0;
    for (int i = 0; i < cn; i++) {
        if (d[i] <= 2) continue;
        long long f = d[i] - 1;
        int bad = 0;
        long long prime = -1;
        for (int j = 0; j < pc; j++) {
            if (f % primes[j] == 0) {
                while (f % primes[j] == 0) {
                    f /= primes[j];
                }
                bad = f != 1;
                prime = primes[j];
                break;
            }
        }
        if (prime < 0) prime = f;
        if (!bad) {
            zd[edgecount++] = {prime, d[i]};
        }
    }
    sort(zd, zd + edgecount);
    for (int i = 0; i < edgecount;) {
        int j = i;
        while (j < edgecount && zd[i].first == zd[j].first) ++j;
        for (int e = 0; e < cn; e++) ndp[e] = dp[e];
        while (i < j) {
            long long t = zd[i].second;
            for (int old = 0, where = 0; old < cn; old++) {
                if (1. * (double) d[old] * (double) t > 1e13) break;
long long newValue = d[old] * t;
while (where < cn && d[where] < newValue) ++where;
                if (where >= cn) break;
if (d[where] == newValue) {
ndp[where] += dp[old];
}
}
i++;
}
for (int e = 0; e < cn; e++) dp[e] = ndp[e];
    }
    cout << dp[cn - 1] << endl;
}