CodeForces 83D Numbers （数论+递归）*

题目链接：http://codeforces.com/problemset/problem/83/D

#include<bits/stdc++.h>
using namespace std;

#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long

#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
const int  maxn =1e3+5;
const int mod=9999991;
const int ub=1e6;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}

/*
题目大意：给定一个范围和一个整数k，
问在范围内，是k的倍数但不能被2到k-1中任意一个数整除
的个数是多少。

令二维状态为（n，k），
我们可以看到在同一个n下，对于不同的k其状态是没有交集的。
那么构造出递归式子：(n,k)=n/k-sigma i=2~k-1 (n/k,i),

对于终止条件，不难看到如果k*k大于n的时候答案为1，
如果n小于k的时候为0，当n大于等于k的时候，如果k不为质数答案明显为0，
注意条件判定顺序。

这个二维状态的收敛速度应该是log级别的，所以时间复杂度应该是够的。
*/

ll a,b,k;

bool check(ll x)
{
    for(int i=2;1LL*i*i<=x;i++)
        if(x%i==0) return false;
    return true;
}

ll solve(ll n,ll k)
{
    ///if(k==1||n==0) return 0LL;
    if(n<k||!check(k)) return 0LL;///在n>=k的情况下话还要再判断一下
    if(k*k>n) return min(n,1LL);
    ll ans=n/k;
    for(int i=2;i<k;i++)
        ans-=solve(n/k,i);
    return ans;
}

int main()
{
    scanf("%lld%lld%lld",&a,&b,&k);
    ///cout<<solve(b,k)<<endl;
    ll ans=solve(b,k)-solve(a-1,k);
    printf("%lld\n",ans);
    return 0;
}