luogu3172 && bzoj3930 [CQOI2015]选数

$题目传送门luogu3172$
$题目传送门bzoj3930$

---

题目

$3930: [CQOI2015]选数$

$Time Limit: 10 Sec$
$Memory Limit: 512 MB$

$Description$

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

$Input$

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

$Output$

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

$Sample Input$

2 2 2 4

$Sample Output$

3

$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\leq N,K\leq 10^{9}，1\leq L\leq H\leq 10^{9}，H-L\leq 10^{5}$

---

题解

看到$gcd$，第一反应这应该是莫反。
然而，这可以用数学方法。
没想到吧，我也没想到。
设$l=\lceil \frac{L}{K} \rceil$，$h=\lfloor \frac{H}{K} \rfloor$。
首先，我们可以把题目转化为从区间$[ l,h ]$中选取$N$个整数，求$N$个整数互质的选取方案有多少个。（因为仅当是$K$的倍数是才可以被选）
我们假设$N$个整数不完全相同。
设$f(i)$表示$i|gcd$的方案数，$g(i)$表示$i=gcd$的方案数。
设$ll=\lceil \frac{l}{i} \rceil$，$hh=\lfloor \frac{h}{i} \rfloor$。
则$f(i)=(hh-ll+1)^{n}-(hh-ll+1)$（因为当且仅当是$i$的倍数时可以被选）
$g(i)=f(i)-\sum\limits_{i|j}{g(j)}$
那么从$H-L$逆推到$1$即可。
$f(i)$用快速幂，$ll=1$时$g(1)$记得$+1$。
时间复杂度$O((H-L)logN)$。

---

总结

这道题其实可以莫反，也可以数学方法。但是我的第一反应仍是数学方法。这种题目如果不会莫反，就不要去惹它。想一想骗分的排列组合，说不定就过了。
不想新建分类，就当做是莫反。

---

标程

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int MOD=1000000007;
const int N=100005;
int n, k, l, h, ans;
int ll, hh, f[N];
int power(int a, int b)
{
    int c=1;
    while (b)
    {
        if (b&1) c=1LL*c*a%MOD;
        b>>=1; a=1LL*a*a%MOD;
    }
    return c;
}
void work()
{
    scanf("%d%d%d%d", &n, &k, &l, &h);
    l=(l-1)/k+1; h/=k;
    for (int i=1; i<=h-l; i++)
    {
        ll=(l-1)/i+1; hh=h/i;
        f[i]=(power(hh-ll+1, n)-(hh-ll+1)+MOD)%MOD;
    }
    for (int i=h-l; i; i--)
        for (int j=(i<<1); j<=h-l; j+=i)
            f[i]=(f[i]-f[j]+MOD)%MOD;
    if (l==1) f[1]=(f[1]+1)%MOD;
    printf("%d\n", f[1]);
}
int main()
{
    work();
    return 0;
}

---