bzoj 3930 CQOI 选数

最新推荐文章于 2019-07-04 22:58:56 发布

我要吃熊猫

最新推荐文章于 2019-07-04 22:58:56 发布

阅读量256

点赞数

分类专栏：杜教筛莫比乌斯反演

本文链接：https://blog.csdn.net/chai_jing/article/details/72854393

版权

莫比乌斯反演同时被 2 个专栏收录

18 篇文章 0 订阅

订阅专栏

杜教筛

6 篇文章 0 订阅

订阅专栏

3930: [CQOI2015]选数
Time Limit: 10 Sec Memory Limit: 512 MB
Submit: 1200 Solved: 573
[Submit][Status][Discuss]
Description

我们知道，从区间[L,H]（L和H为整数）中选取N个整数，总共有(H-L+1)^N种方案。小z很好奇这样选出的数的最大公约数的规律，他决定对每种方案选出的N个整数都求一次最大公约数，以便进一步研究。然而他很快发现工作量太大了，于是向你寻求帮助。你的任务很简单，小z会告诉你一个整数K，你需要回答他最大公约数刚好为K的选取方案有多少个。由于方案数较大，你只需要输出其除以1000000007的余数即可。

Input

输入一行，包含4个空格分开的正整数，依次为N，K，L和H。

Output

输出一个整数，为所求方案数。

Sample Input

2 2 2 4

Sample Output

HINT

样例解释

所有可能的选择方案：(2, 2), (2, 3), (2, 4), (3, 2), (3, 3), (3, 4), (4, 2), (4, 3), (4, 4)

其中最大公约数等于2的只有3组：(2, 2), (2, 4), (4, 2)

对于100%的数据，1≤N,K≤10^9，1≤L≤H≤10^9，H-L≤10^5

Source

【分析】
不知道大家还记不记得Problem B…
异曲同工…繁衍还是当年的繁衍，加了个快速幂（平方变为n次方）
然后也不用当年傻傻的拆分询问了…QAQ
做四次反演++- -都是骗人的

【代码】

//bzoj 3930 CQOI 选数
#include<map>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define N 10000000
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mxn=10000005;
const int mod=1e9+7;
bool vis[mxn];
map <int,int> M;
int n,pri[mxn],mu[mxn];
inline void init()
{
    mu[1]=1;
    int i,j;
    fo(i,2,N)
    {
        if(!vis[i]) pri[++pri[0]]=i,mu[i]=-1;
        for(j=1;j<=N && (ll)i*pri[j]<=N;j++)
        {
            vis[i*pri[j]]=1;
            if(i%pri[j]==0) break;
            mu[i*pri[j]]=-mu[i];
        }
    }
    fo(i,1,N) mu[i]+=mu[i-1];
}
inline int res(int x)
{
    if(x<=N) return mu[x];
    if(M.count(x)) return M[x];
    ll sum=1;
    for(int i=2,last=0;i<=x;i=last+1)
    {
        last=x/(x/i);
        sum=(sum-(ll)(last-i+1)*res(x/i)%mod)%mod;
    }
    sum=(sum+mod)%mod;
    M[x]=(int)sum;
    return (int)sum;
}
inline ll ksm(ll x,int k)
{
    if(k==1) return x;
    ll tmp=ksm(x,k>>1);
    if(k&1) return tmp*tmp%mod*x%mod;
    return tmp*tmp%mod;
}
inline ll solve(int L,int H)
{
    ll sum=0;
    for(int i=1,last=0;i<=H;i=last+1)
    {
        if(L/i==0) last=H/(H/i);
        else last=min(L/(L/i),H/(H/i));
        sum=(sum+(ll)(res(last)-res(i-1))%mod*ksm(H/i-L/i,n)%mod)%mod;
    }
    sum=(sum+mod)%mod;
    return sum;
}
int main()
{
    init();
    int i,j,k,L,H;
    scanf("%d%d%d%d",&n,&k,&L,&H);
    L=(L-1)/k,H=H/k;
    printf("%lld\n",solve(L,H));
    return 0;
}
//2 1 
//2 1000000000