sacgajcvs

Untitled

Feb 9th, 2020
185
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.67 KB | None | 0 0
  1. /*
  2. _____ _ _ _ _
  3. |_ _| |__ ___ / \ _ __ ___| |__ _ _| |
  4. | | | '_ \ / _ \ / _ \ | '_ \/ __| '_ \| | | | |
  5. | | | | | | __// ___ \| | | \__ \ | | | |_| | |
  6. |_| |_| |_|\___/_/ \_\_| |_|___/_| |_|\__,_|_|
  7.  
  8. */
  9. #include<bits/stdc++.h>
  10. #include <ext/pb_ds/assoc_container.hpp>
  11. #include <ext/pb_ds/tree_policy.hpp>
  12. #define ll long long
  13. #define pb push_back
  14. #define ppb pop_back
  15. #define endl '\n'
  16. #define mii map<ll,ll>
  17. #define msi map<string,ll>
  18. #define mis map<ll, string>
  19. #define rep(i,a,b) for(ll i=a;i<b;i++)
  20. #define repr(i,a,b) for(ll i=b-1;i>=a;i--)
  21. #define trav(a, x) for(auto& a : x)
  22. #define pii pair<ll,ll>
  23. #define vi vector<ll>
  24. #define vii vector<pair<ll, ll>>
  25. #define vs vector<string>
  26. #define all(a) (a).begin(),(a).end()
  27. #define F first
  28. #define S second
  29. #define sz(x) (ll)x.size()
  30. #define hell 1000000007
  31. #define lbnd lower_bound
  32. #define ubnd upper_bound
  33. #define max(a,b) (a>b?a:b)
  34. #define min(a,b) (a<b?a:b)
  35.  
  36. /* For Debugging */
  37. #define DEBUG cerr<<"\n>>>I'm Here<<<\n"<<endl;
  38. #define display(x) trav(a,x) cout<<a<<" ";cout<<endl;
  39. #define what_is(x) cerr << #x << " is " << x << endl;
  40.  
  41. std::mt19937_64 rng(std::chrono::steady_clock::now().time_since_epoch().count());
  42. #define ordered_set tree<ll, null_type,less<ll>, rb_tree_tag,tree_order_statistics_node_update>
  43. #define TIME cerr << "\nTime elapsed: " << setprecision(5) <<1000.0 * clock() / CLOCKS_PER_SEC << "ms\n";
  44. #define DECIMAL(n) cout << fixed ; cout << setprecision(n);
  45. #define FAST ios_base::sync_with_stdio(false);cin.tie(0);cout.tie(0);
  46. using namespace __gnu_pbds;
  47. using namespace std;
  48. #define PI 3.141592653589793
  49. #define N 100005
  50. ll add(ll a,ll b,ll mod=hell)
  51. {
  52. return (a+b)%mod;
  53. }
  54. ll sub(ll a,ll b,ll mod=hell)
  55. {
  56. return (a-b+mod)%mod;
  57. }
  58. ll mul(ll a,ll b,ll mod=hell)
  59. {
  60. a=a%hell;
  61. b=b%hell;
  62. return (a*b)%mod;
  63. }
  64. ll expo(ll base, ll exponent, ll mod) { //return base^exponent modulo modulus
  65. ll ans = 1;
  66. while(exponent !=0 ) {
  67. if((exponent&1) == 1) {
  68. ans = ans*base ;
  69. ans = ans%mod;
  70. }
  71. base = base*base;
  72. base %= mod;
  73. exponent>>= 1;
  74. }
  75. return ans%mod;
  76. }
  77. void solve()
  78. {
  79. ll n;
  80. cin>>n;
  81. ll ans=0,vl;
  82. ll modi=expo(2,hell-2,hell);
  83. for (ll i = 1, la; i <= n; i = la + 1)
  84. {
  85. la = n / (n / i);
  86. vl=mul(sub(mul(la,(la+1)),mul(i,(i-1+hell))),modi);
  87. vl=mul(vl,(n/i));
  88. ans=add(ans,vl);
  89. }
  90. cout<<ans<<endl;
  91. return;
  92. }
  93. int main()
  94. {
  95. FAST
  96. int TESTS=1;
  97. // cin>>TESTS;
  98. while(TESTS--)
  99. {
  100. solve();
  101. }
  102. TIME
  103. return 0;
  104. }
Advertisement
Add Comment
Please, Sign In to add comment