hdu6169 Senior PanⅡ【2017多校第九场】

转载请注明出处，谢谢http://blog.csdn.net/bigtiao097?viewmode=contents
参考文章http://blog.csdn.net/wubaizhe/article/details/77484454
题目传送门

题意：

给定区间[L,R]和整数K，求区间内所有满足其最小因子(1除外)为K的数的和
$1\leq L \leq R \leq 10^{11}$、$2\leq K \leq 10^{11}$

思路：

问题可以转化为求$1~n$中最小因子（1除外）为K的数的和
首先，$k$不是素数时，和为0
其次，$k>\sqrt n 时$，和为k，由于$n<10^{11}$,故只需算$k<320000$的情况

$dp[i][j]$表示1~i 这 i个数被前j个素数筛过以后剩下的数的和。
则：

$$dp[i][j]=dp[i][j−1]−prime[j]∗dp[i/prime[j]][j−1]$$
解释：对于 $dp[i][j]$，只需要考虑在 $dp[i][j−1]$的基础上减去第j个素数造成的影响。
对于 $[1,i]$中会被第 $j个素数prime[j]$筛去的数来说，其可以表示成 $x∗prime[j]$
其中 $x\leq i/prime[j]$，而且x不会被前j−1个素数筛去，故 $dp[i/prime[j]][j−1]$中的每一个数，乘上 $prime[j]$之后一定是会被第j个素数筛去的数，故递推式即为所求。

然后最终答案就等于：(若K为第cnt+1个素数)

$$K∗dp[R/K][cnt]−K∗dp[(L−1)/K][cnt]$$
因为L,R的范围很大，考虑小规模记忆化，大规模搜索的DP策略。

---

具体代码如下：
Result：Accepted$\ \ \ $Memory: 12480K$\ \ \ $Time : 670MS

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod = 1e9+7;
const int maxn = 3.2e5+5;
const ll inv2 = 500000004;
const int maxa = 1e4+5;
const int maxb = 1e2+5;
ll dp[maxa][maxb];
ll l,r,k;
int T,cases;
const int maxp = 3.2e5+5;
int prime[maxp] = {0},cnt = 0;
bool notPri[maxp] = {1,1};
void linearGetPrimes()
{
    for(int i = 2; i < maxp; i++)
    {
        if(!notPri[i]) prime[++cnt] = i;
        for(int j = 1; j <= cnt && (long long)i*prime[j] < maxp; j++)
        {
            notPri[i*prime[j]] = 1;
            if(!(i%prime[j])) break;
        }
    }
}
void init()
{
    linearGetPrimes();
    for(int i=1;i<maxa;i++)
    {
        dp[i][0] = i*(i+1)/2;
        for(int j=1;j<maxb;j++)
            dp[i][j] = dp[i][j-1]-prime[j]*dp[i/prime[j]][j-1];
    }

}

ll calc(ll n,ll m)
{
    if(n<=1)return n;
    if(!m) return n%mod*(n%mod+1)%mod*inv2%mod;
    if(n<maxa&&m<maxb)return dp[n][m];
    while(m&&prime[m]>n) m--;
    return (calc(n,m-1)-1LL*prime[m]*calc(n/prime[m],m-1)%mod+mod)%mod;
}
bool isprime(ll x)
{
    for(ll i=2;i*i<=x;i++)
        if(x%i==0)return 0;
    return 1;
}
int main()
{
    init();
    ios::sync_with_stdio(false);
    cin>>T;
    while(T--)
    {
        cout<<"Case #"<<++cases<<": ";
        cin>>l>>r>>k;
        if(!isprime(k))
        {
            cout<<0<<endl;
            continue;
        }
        if(k>maxn)
        {
            if(l<=k&&r>=k)
                cout<<k%mod<<endl;
            else
                cout<<0<<endl;
            continue;
        }
        int x = 0;
        while(prime[x+1]<k) x++;
        cout<<(calc(r/k,x)*k%mod-calc((l-1)/k,x)*k%mod+mod)%mod<<endl;

    }
}