首页 > 代码库 > BZOJ 2186 SDOI 2008 沙拉公主的困惑 数论

BZOJ 2186 SDOI 2008 沙拉公主的困惑 数论

题目大意:求出1~N!中与M!互质的数的个数,保证N>M。


前言:最讨厌数学了。。。

思路:因为保证了n>m,所以n!一定是m!的倍数。如果找到了一个x使得gcd(x,m!)==1,那么gcd(x + m!,m!) == 1一定成立,gcd(x + k * m!,m!) == 1(k >= 1)也一定成立。x的个数就是φ(m!),那么总的个数就是φ(m!) * n! / m!。之后各种推公式,就可以用线性筛来解决这个问题。(我太弱,讲不明白。。。


CODE:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define MAX 10000100
using namespace std;

int T,P,n,m,primes;
long long prime[500100],ans[MAX],_link[MAX],rev[MAX];
bool notp[MAX];

void Pretreatment()
{
	for(int i = 2; i < MAX; ++i) {
		if(!notp[i])	prime[++primes] = i;
		for(int j = 1;j <= primes && prime[j] * i < MAX; ++j) {
			notp[prime[j] * i] = true;
			if(i % prime[j] == 0)
				break;
		}
	}
	_link[0] = 1;
	for(int i = 1; i < MAX; ++i)
		_link[i] = _link[i - 1] * i % P;
	rev[1] = 1;
	for(int i = 2; i < MAX && i < P; ++i)
		rev[i] = (P - P / i) * rev[P % i] % P;
	ans[1] = 1;
	for(int i = 2; i < MAX; ++i) {
		if(!notp[i])
			ans[i] = ans[i - 1] * (i - 1) % P * rev[i % P] % P;
		else	ans[i] = ans[i - 1];
	}
}

int main()
{
	cin >> T >> P;
	Pretreatment();
	for(int i = 1; i <= T; ++i) {
		scanf("%d%d",&n,&m);
		printf("%d\n",_link[n] * ans[m] % P);
	}
	return 0;
}


BZOJ 2186 SDOI 2008 沙拉公主的困惑 数论