Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- /*
- _____ _ _ _ _
- |_ _| |__ ___ / \ _ __ ___| |__ _ _| |
- | | | '_ \ / _ \ / _ \ | '_ \/ __| '_ \| | | | |
- | | | | | | __// ___ \| | | \__ \ | | | |_| | |
- |_| |_| |_|\___/_/ \_\_| |_|___/_| |_|\__,_|_|
- */
- #include<bits/stdc++.h>
- #include <ext/pb_ds/assoc_container.hpp>
- #include <ext/pb_ds/tree_policy.hpp>
- #define ll long long
- #define pb push_back
- #define ppb pop_back
- #define endl '\n'
- #define mii map<ll,ll>
- #define msi map<string,ll>
- #define mis map<ll, string>
- #define rep(i,a,b) for(ll i=a;i<b;i++)
- #define repr(i,a,b) for(ll i=b-1;i>=a;i--)
- #define trav(a, x) for(auto& a : x)
- #define pii pair<ll,ll>
- #define vi vector<ll>
- #define vii vector<pair<ll, ll>>
- #define vs vector<string>
- #define all(a) (a).begin(),(a).end()
- #define F first
- #define S second
- #define sz(x) (ll)x.size()
- #define hell 1000000007
- #define lbnd lower_bound
- #define ubnd upper_bound
- #define max(a,b) (a>b?a:b)
- #define min(a,b) (a<b?a:b)
- /* For Debugging */
- #define DEBUG cerr<<"/n>>>I'm Here<<</n"<<endl;
- #define display(x) trav(a,x) cout<<a<<" ";cout<<endl;
- #define what_is(x) cerr << #x << " is " << x << endl;
- std::mt19937_64 rng(std::chrono::steady_clock::now().time_since_epoch().count());
- #define ordered_set tree<ll, null_type,less<ll>, rb_tree_tag,tree_order_statistics_node_update>
- #define TIME cerr << "\nTime elapsed: " << setprecision(5) <<1000.0 * clock() / CLOCKS_PER_SEC << "ms\n";
- #define DECIMAL(n) cout << fixed ; cout << setprecision(n);
- #define FAST ios_base::sync_with_stdio(false);cin.tie(0);cout.tie(0);
- using namespace __gnu_pbds;
- using namespace std;
- #define PI 3.141592653589793
- #define N 100005
- ll l,r,n;
- ll dp[100];
- ll fun(ll num)
- {
- return (r/num-(l-1)/num);
- }
- void solve(ll test)
- {
- cin>>l>>r>>n;
- ll tmp=n;
- vi vec;
- for(ll i=2;i*i<=n;i++)
- {
- if(tmp%i==0)
- {
- vec.pb(i);
- while(tmp%i==0)
- tmp/=i;
- }
- }
- if(tmp>1)
- vec.pb(tmp);
- rep(i,0,sz(vec))
- dp[i]=0;
- rep(i,1,(1<<sz(vec)))
- {
- ll num=1;
- rep(j,0,sz(vec))
- {
- if(i&(1<<j))
- num*=vec[j];
- }
- dp[__builtin_popcount(i)-1]+=fun(num);
- }
- ll ans=0;
- rep(i,0,sz(vec))
- {
- if(!(i&1))
- ans+=dp[i];
- else
- ans-=dp[i];
- }
- cout<<"Case #"<<test<<": "<<r-l+1-ans<<endl;
- return;
- }
- int main()
- {
- FAST
- // #ifndef ONLINE_JUDGE
- // // freopen ("1","r",stdin);
- // // freopen("output.txt","w",stdout);
- // #endif
- int TESTS=1;
- cin>>TESTS;
- rep(i,1,TESTS+1)
- {
- solve(i);
- }
- // TIME
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment