#include<cstdio>
#include<iostream>
#include<algorithm>
#include<string>
#include<vector>
#include<set>
#include<cmath>
using namespace std;
typedef vector<int> VI;
typedef long long LL;
#define FOR(x, b, e) for(int x=b; x<=(e); ++x)
#define FORD(x, b, e) for(int x=b; x>=(e); --x)
#define REP(x, n) for(int x=0; x<(n); ++x)
#define VAR(v, n) __typeof(n) v = (n)
#define ALL(c) (c).begin(), (c).end()
#define SIZE(x) ((int)(x).size())
#define FOREACH(i, c) for(VAR(i, (c).begin()); i != (c).end(); ++i)
#define PB push_back
#define ST first
#define ND second
class Primes {
int n;
VI primes;
public:
VI primeList(int n) {
VI w(1, 2);
int s = 0, i = 2;
FOR(l, 3, n - 1) {
i = 0;
while (w[s] * w[s] <= l) s++;
while (i < s && l % w[i]) i++;
if (i == s) w.PB(l);
}
return w;
}
Primes(int n): n(n) {
primes = primeList(n);
}
vector<pair<int, int> > getPrimesFactorial(int a) {
vector<pair<int, int> > result;
if(a==1) {
result.PB(make_pair(0,0));
return result;
}
int i = 0;
while(a > 1 && i < primes.size()) {
int d = 0;
while(a % primes[i] == 0) {
a /= primes[i];
d++;
}
if(d>0) result.PB(make_pair(primes[i], d));
i++;
}
if(a>1) result.PB(make_pair(a, 1));
return result;
}
set<int> getDivisors(int a) {
set<int> result;
for(int i=1;i<=sqrt(a);i++) {
if(a%i==0) {
result.insert(i);
result.insert(a/i);
}
}
return result;
}
int getDivisorsCount(int a) {
vector<pair<int, int> > prime_fact = getPrimesFactorial(a);
int result = 1;
for(int i=0;i<prime_fact.size();i++) {
result *= prime_fact[i].second + 1;
}
return result;
}
};
class Solver {
int n;
public:
void readInput() {
cin >> n;
}
void solve() {
if(n < 6) {
cout << 0 << '\n';
return;
}
Primes p((int)ceil(sqrt(n))+1);
LL result = 0;
for(auto d:p.getDivisors(n)) {
if(n/d > 1) {
int tmp = p.getDivisorsCount(n/d - 1) - 2;
if((n/d-1)%2 == 0) tmp--;
result += max(0,tmp);
}
}
cout << result << '\n';
}
};
int main() {
ios_base::sync_with_stdio(0);
Solver solver;
solver.readInput();
solver.solve();
return 0;
}
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 | #include<cstdio> #include<iostream> #include<algorithm> #include<string> #include<vector> #include<set> #include<cmath> using namespace std; typedef vector<int> VI; typedef long long LL; #define FOR(x, b, e) for(int x=b; x<=(e); ++x) #define FORD(x, b, e) for(int x=b; x>=(e); --x) #define REP(x, n) for(int x=0; x<(n); ++x) #define VAR(v, n) __typeof(n) v = (n) #define ALL(c) (c).begin(), (c).end() #define SIZE(x) ((int)(x).size()) #define FOREACH(i, c) for(VAR(i, (c).begin()); i != (c).end(); ++i) #define PB push_back #define ST first #define ND second class Primes { int n; VI primes; public: VI primeList(int n) { VI w(1, 2); int s = 0, i = 2; FOR(l, 3, n - 1) { i = 0; while (w[s] * w[s] <= l) s++; while (i < s && l % w[i]) i++; if (i == s) w.PB(l); } return w; } Primes(int n): n(n) { primes = primeList(n); } vector<pair<int, int> > getPrimesFactorial(int a) { vector<pair<int, int> > result; if(a==1) { result.PB(make_pair(0,0)); return result; } int i = 0; while(a > 1 && i < primes.size()) { int d = 0; while(a % primes[i] == 0) { a /= primes[i]; d++; } if(d>0) result.PB(make_pair(primes[i], d)); i++; } if(a>1) result.PB(make_pair(a, 1)); return result; } set<int> getDivisors(int a) { set<int> result; for(int i=1;i<=sqrt(a);i++) { if(a%i==0) { result.insert(i); result.insert(a/i); } } return result; } int getDivisorsCount(int a) { vector<pair<int, int> > prime_fact = getPrimesFactorial(a); int result = 1; for(int i=0;i<prime_fact.size();i++) { result *= prime_fact[i].second + 1; } return result; } }; class Solver { int n; public: void readInput() { cin >> n; } void solve() { if(n < 6) { cout << 0 << '\n'; return; } Primes p((int)ceil(sqrt(n))+1); LL result = 0; for(auto d:p.getDivisors(n)) { if(n/d > 1) { int tmp = p.getDivisorsCount(n/d - 1) - 2; if((n/d-1)%2 == 0) tmp--; result += max(0,tmp); } } cout << result << '\n'; } }; int main() { ios_base::sync_with_stdio(0); Solver solver; solver.readInput(); solver.solve(); return 0; } |
English