Advertisement
erfanul007

LA 2195

Jul 5th, 2021
1,120
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 4.76 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. // #include <iostream>
  3. // #include <cstdio>
  4. // #include <cstdlib>
  5. // #include <algorithm>
  6. // #include <cmath>
  7. // #include <vector>
  8. // #include <set>
  9. // #include <map>
  10. // #include <queue>
  11. // #include <stack>
  12. // #include <ctime>
  13. // #include <cassert>
  14. // #include <complex>
  15. // #include <string>
  16. // #include <cstring>
  17. // #include <bitset>
  18. using namespace std;
  19.  
  20. // #pragma GCC optimize("Ofast,no-stack-protector")
  21. // #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
  22. // #pragma GCC optimize("unroll-loops")
  23.  
  24. #define ll              long long int
  25. #define vi              vector< int >
  26. #define vll             vector< ll >
  27.  
  28. #define sc              scanf
  29. #define pf              printf
  30. #define cspf(i)         pf("Case %d:\n", i)
  31. #define spc             pf(" ")
  32. #define line            pf("\n")
  33.  
  34. #define ff              first
  35. #define ss              second
  36. #define mp              make_pair
  37. #define pb              push_back
  38. #define ppb             pop_back
  39. #define tp(v,j)         get<j>(v)
  40. #define Log(b,x)        (log(x)/log(b))
  41.  
  42. #define FOR(i,x,y)      for(int i = int(x); i < int(y); i++)
  43. #define ROF(i,x,y)      for(int i = int(x)-1; i >= int(y); i--)
  44. #define clr(arr,x)      memset(arr, x, sizeof arr)
  45. #define vout(v,sz)      for(int w=0;w<sz;w++){if(w) spc; cout<<v[w].ff << ' ' <<v[w].ss<<'\n';}
  46. #define all(v)          v.begin(), v.end()
  47. #define rall(v)         v.rbegin(), v.rend()
  48. #define unq(v)          sort(all(v)),(v).resize(unique(all(v))-v.begin())
  49. #define fastIO          ios::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr)
  50.  
  51. #define sc1(x)          sc("%d",&x)
  52. #define sc2(x,y)        sc("%d %d", &x, &y)
  53. #define sc3(x,y,z)      sc("%d %d %d", &x, &y, &z)
  54. #define scl1(x)         sc("%lld",&x)
  55. #define scl2(x,y)       sc("%lld %lld", &x, &y)
  56. #define scf1(x)         sc("%lf",&x);
  57. #define scf2(x,y)       sc("%lf %lf", &x, &y)
  58.  
  59. #define pf1(x)          pf("%d",x)
  60. #define pf2(x,y)        pf("%d %d", x, y)
  61. #define pfl1(x)         pf("%lld",x)
  62. #define pfl2(x,y)       pf("%lld %lld", x, y)
  63.  
  64. #define MOD             (int)(1e9)
  65. #define MaxN            4000000
  66. #define inf             0x3f3f3f3f
  67. #define PI              acos(-1.0)  // 3.1415926535897932
  68. #define eps             1e-6
  69.  
  70. #ifdef ERFANUL007
  71.     #define debug(...) __f(#__VA_ARGS__, __VA_ARGS__)
  72.     template < typename Arg1 >
  73.     void __f(const char* name, Arg1&& arg1){
  74.         cout << name << " = " << arg1 << std::endl;
  75.     }
  76.     template < typename Arg1, typename... Args>
  77.     void __f(const char* names, Arg1&& arg1, Args&&... args){
  78.         const char* comma = strchr(names, ',');
  79.         cout.write(names, comma - names) << " = " << arg1 <<" | ";
  80.         __f(comma+1, args...);
  81.     }
  82. #else
  83.     #define debug(...)
  84. #endif
  85.  
  86. template <class T> inline T bigMod(T p,T e,T M){T ret=1; for(;e>0;e>>=1){ if(e&1) ret=(ret*p)%M; p=(p*p)%M;} return (T)ret;}
  87. template <class T> inline T modInverse(T a,T M){return bigMod(a,M-2,M);}
  88. template <class T> inline T gcd(T a,T b){if(b==0)return a;return gcd(b,a%b);}
  89. template <class T> inline T lcm(T a,T b) {a=abs(a);b=abs(b); return (a/gcd(a,b))*b;}
  90.  
  91. int dx[] = { 1,-1, 0, 0};                //graph moves
  92. int dy[] = { 0, 0, 1,-1};               //graph moves
  93.  
  94. bool isprime[MaxN];
  95. vector< ll > prime;
  96.  
  97. void sieve(){
  98.     clr(isprime, 1);
  99.     for(int i=4; i<MaxN; i+=2) isprime[i] = 0;
  100.     for(ll i=3; i*i<MaxN; i+=2){
  101.         if(!isprime[i]) continue;
  102.         for(ll j = i*i; j<MaxN; j+=(i*2)){
  103.             isprime[j] = 0;
  104.         }
  105.     }
  106.     prime.pb(2);
  107.     for(int i=3; i<MaxN; i+=2){
  108.         if(isprime[i]) prime.pb(i);
  109.     }
  110. }
  111.  
  112. vector< pair< ll, ll > > fact;
  113. void primefact(ll n){
  114.     for(int i=0; prime[i] * prime[i] <= n; i++){
  115.         if(n % prime[i] == 0){
  116.             ll cnt = 0;
  117.             while(n % prime[i] == 0){
  118.                 n/=prime[i];
  119.                 cnt++;
  120.             }
  121.             fact.pb(mp(prime[i], cnt));
  122.         }
  123.     }
  124.     if(n > 1) fact.pb(mp(n, 1));
  125. }
  126.  
  127. vector< ll > divs;
  128. void BT(int i, ll val){
  129.     if(i == fact.size()){
  130.         divs.pb(val);
  131.         return;
  132.     }
  133.     BT(i+1, val);
  134.  
  135.     for(int j=0; j<fact[i].ss; j++){
  136.         val *= fact[i].ff;
  137.         BT(i+1, val);
  138.     }
  139. }
  140.  
  141. int solve(ll n){
  142.     fact.clear();
  143.     primefact(n);
  144.     divs.clear();
  145.     BT(0, 1);
  146.     int sum = 0;
  147.     for(int i=1; i<divs.size(); i++){
  148.         ll m = n;
  149.         while(m % divs[i] == 0){
  150.             m /= divs[i];
  151.             sum++;
  152.         }
  153.     }
  154.     return sum;
  155. }
  156.  
  157. int main()
  158. {
  159.     #ifdef ERFANUL007
  160.         clock_t tStart = clock();
  161.         freopen("input.txt", "r", stdin);
  162.         freopen("output.txt", "w", stdout);
  163.     #endif
  164.  
  165.     sieve();
  166.  
  167.     ll n;
  168.  
  169.     while(~scl1(n) && n){
  170.         pfl1(n); spc;
  171.         pf1(solve(n)); line;
  172.     }
  173.  
  174.  
  175.     #ifdef ERFANUL007
  176.         fprintf(stderr, "\n>> Runtime: %.10fs\n", (double) (clock() - tStart) / CLOCKS_PER_SEC);
  177.     #endif
  178.  
  179.     return 0;
  180. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement