sacgajcvs

Untitled

Dec 23rd, 2019
207
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.66 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 l,r,n;
  51. ll dp[100];
  52. ll fun(ll num)
  53. {
  54. return (r/num-(l-1)/num);
  55. }
  56. void solve(ll test)
  57. {
  58. cin>>l>>r>>n;
  59. ll tmp=n;
  60. vi vec;
  61. for(ll i=2;i*i<=n;i++)
  62. {
  63. if(tmp%i==0)
  64. {
  65. vec.pb(i);
  66. while(tmp%i==0)
  67. tmp/=i;
  68. }
  69. }
  70. if(tmp>1)
  71. vec.pb(tmp);
  72. rep(i,0,sz(vec))
  73. dp[i]=0;
  74. rep(i,1,(1<<sz(vec)))
  75. {
  76. ll num=1;
  77. rep(j,0,sz(vec))
  78. {
  79. if(i&(1<<j))
  80. num*=vec[j];
  81. }
  82. dp[__builtin_popcount(i)-1]+=fun(num);
  83. }
  84. ll ans=0;
  85. rep(i,0,sz(vec))
  86. {
  87. if(!(i&1))
  88. ans+=dp[i];
  89. else
  90. ans-=dp[i];
  91. }
  92. cout<<"Case #"<<test<<": "<<r-l+1-ans<<endl;
  93. return;
  94. }
  95. int main()
  96. {
  97. FAST
  98. // #ifndef ONLINE_JUDGE
  99. // // freopen ("1","r",stdin);
  100. // // freopen("output.txt","w",stdout);
  101. // #endif
  102. int TESTS=1;
  103. cin>>TESTS;
  104. rep(i,1,TESTS+1)
  105. {
  106. solve(i);
  107. }
  108. // TIME
  109. return 0;
  110. }
Advertisement
Add Comment
Please, Sign In to add comment