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

标签：bzoj   数学   数论   

题目大意：求出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 沙拉公主的困惑 数论

标签：bzoj   数学   数论   

原文地址：http://blog.csdn.net/jiangyuze831/article/details/40426615