Advertisement
Saleh127

UVA 12396

May 11th, 2021
110
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.30 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define ll long long
  4. #define test int t; cin>>t; for(int cs=1;cs<=t;cs++)
  5. #define maX 10000002
  6. vector<ll>p;
  7. ll mod=1e9+7;
  8. bool marked[maX+2];
  9. void sieve()
  10. {
  11. marked[0]=1;
  12. marked[1]=1;
  13. for(ll i=4; i<=maX; i+=2)
  14. {
  15. marked[i]=1;
  16. }
  17. p.push_back(2);
  18. for(ll i=3; i<=maX; i+=2)
  19. {
  20. if(marked[i]==0)
  21. {
  22. p.push_back(i);
  23. for(ll j=i*i; j<=maX; j+=i+i)
  24. {
  25. marked[j]=1;
  26. }
  27. }
  28. }
  29. }
  30.  
  31. ll bigmod(ll a,ll c,ll d)
  32. {
  33. if(c==0) return 1LL;
  34. ll x=bigmod(a,c/2,d);
  35. x=(x*x)%d;
  36. if(c%2==1LL)
  37. {
  38. x=(x*a)%d;
  39. }
  40. return x;
  41. }
  42.  
  43. int main()
  44. {
  45. ios_base::sync_with_stdio(0);
  46. cin.tie(0);
  47. cout.tie(0);
  48.  
  49. sieve();
  50.  
  51. ll n;
  52.  
  53. while(cin>>n)
  54. {
  55. if(n==0) break;
  56.  
  57. ll i,j,k,c,ans=1;
  58.  
  59. for(i=0; i<p.size() && p[i]<=n/2; i++)
  60. {
  61. c=0;
  62. k=p[i];
  63. while(k<=n)
  64. {
  65. c+=(n/k);
  66. k*=p[i];
  67. }
  68. if(c<=1) break;
  69.  
  70. if(c%2) c--;
  71.  
  72. ans*=bigmod(p[i],c,mod);
  73.  
  74. ans%=mod;
  75. }
  76. cout<<ans<<endl;
  77. }
  78.  
  79. return 0;
  80. }
  81.  
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement