#include <cstdio> #include <cstdlib> #include <iostream> #include <fstream> #include <sstream> #include <set> #include <map> #include <vector> #include <list> #include <algorithm> #include <cstring> #include <cmath> #include <string> #include <queue> #include <bitset> //UWAGA - w czasie kompilacji musi byc znany rozmiar wektora - nie mozna go zmienic #include <cassert> #include <iomanip> //do setprecision #include <ctime> #include <complex> using namespace std; #define FOR(i,b,e) for(int i=(b);i<(e);++i) #define FORQ(i,b,e) for(int i=(b);i<=(e);++i) #define FORD(i,b,e) for(int i=(b)-1;i>=(e);--i) #define REP(x, n) for(int x = 0; x < (n); ++x) #define ALL(u) (u).begin(),(u).end() #define ST first #define ND second #define PB push_back #define MP make_pair #define LL long long #define ULL unsigned LL #define LD long double typedef pair<int, int> PII; const double pi = 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342; const int MR = 32000; bool noPrime[MR]; vector<PII> fact(int n, const vector<int> &primes) { vector<PII> res; for (int p : primes) { if (p > sqrt(n)) break; int cnt = 0; while (n % p == 0) { n /= p; cnt++; } if (cnt) res.push_back(MP(p, cnt)); } if (n > 1) res.push_back(MP(n, 1)); return res; } int main() { int n; scanf("%d", &n); vector<int> div; int lim = sqrt(n); vector<int> primes; { int lim1 = sqrt(lim); noPrime[0] = noPrime[1] = 1; FORQ(i, 2, lim1) if (!noPrime[i]) { for (int j = i * i; j <= lim; j += i) noPrime[j] = 1; } FORQ(i, 2, lim) if (!noPrime[i]) primes.push_back(i); } FORQ(i, 1, lim) if (n%i == 0) { div.push_back(i); if (n / i != i) div.push_back(n / i); } LL res = 0; for (int a : div) { int r = n / a; r--; if (r < 6) continue; auto pr = fact(r, primes); int s = 1; for (const auto &p : pr) s *= p.second + 1; // dodaj kazdy rozklad - wyliczony z liczby wszystkich dzielnikow // oprocz 1 i r s -= 2; // kazdy dzielnik daje dobra pare oprocz potencjalnie ostatniego // z+1 > 2 czyli jak r parzyste to ostatni dzielnik da nam z+1 == 2 // np dla 12, ostatni dzielnik = 6, to y = 6, a z = 12/6 = 2 res += s - (1 - r % 2); } printf("%lld\n", res); 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 121 122 123 124 125 126 127 128 129 | #include <cstdio> #include <cstdlib> #include <iostream> #include <fstream> #include <sstream> #include <set> #include <map> #include <vector> #include <list> #include <algorithm> #include <cstring> #include <cmath> #include <string> #include <queue> #include <bitset> //UWAGA - w czasie kompilacji musi byc znany rozmiar wektora - nie mozna go zmienic #include <cassert> #include <iomanip> //do setprecision #include <ctime> #include <complex> using namespace std; #define FOR(i,b,e) for(int i=(b);i<(e);++i) #define FORQ(i,b,e) for(int i=(b);i<=(e);++i) #define FORD(i,b,e) for(int i=(b)-1;i>=(e);--i) #define REP(x, n) for(int x = 0; x < (n); ++x) #define ALL(u) (u).begin(),(u).end() #define ST first #define ND second #define PB push_back #define MP make_pair #define LL long long #define ULL unsigned LL #define LD long double typedef pair<int, int> PII; const double pi = 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342; const int MR = 32000; bool noPrime[MR]; vector<PII> fact(int n, const vector<int> &primes) { vector<PII> res; for (int p : primes) { if (p > sqrt(n)) break; int cnt = 0; while (n % p == 0) { n /= p; cnt++; } if (cnt) res.push_back(MP(p, cnt)); } if (n > 1) res.push_back(MP(n, 1)); return res; } int main() { int n; scanf("%d", &n); vector<int> div; int lim = sqrt(n); vector<int> primes; { int lim1 = sqrt(lim); noPrime[0] = noPrime[1] = 1; FORQ(i, 2, lim1) if (!noPrime[i]) { for (int j = i * i; j <= lim; j += i) noPrime[j] = 1; } FORQ(i, 2, lim) if (!noPrime[i]) primes.push_back(i); } FORQ(i, 1, lim) if (n%i == 0) { div.push_back(i); if (n / i != i) div.push_back(n / i); } LL res = 0; for (int a : div) { int r = n / a; r--; if (r < 6) continue; auto pr = fact(r, primes); int s = 1; for (const auto &p : pr) s *= p.second + 1; // dodaj kazdy rozklad - wyliczony z liczby wszystkich dzielnikow // oprocz 1 i r s -= 2; // kazdy dzielnik daje dobra pare oprocz potencjalnie ostatniego // z+1 > 2 czyli jak r parzyste to ostatni dzielnik da nam z+1 == 2 // np dla 12, ostatni dzielnik = 6, to y = 6, a z = 12/6 = 2 res += s - (1 - r % 2); } printf("%lld\n", res); return 0; } |