沙拉公主的困惑（bzoj2186 luogu2155）

这是一道很恶心的题目，老实说我第一次对着标程调了2个多小时

题目不说了，相信肯定看完了

显然N!是M!的倍数，答案为

预处理出1~10^7的质数及逆元还有阶乘

递推求

而答案为

该怎么递推求解f呢？

            

当m为质数 

       否则

这道题就这么简单的ac了

但是细节很多，如果你都开long long（P党不管了），你会爆空间，如果你用高精，你会贼烦。

所以记得%p，并且强制转long long

#include<cstdio>
#define ll long long
using namespace std;
const int N=1e7+5;
const int n=1e7;
int p,cnt;
int pri[N],jc[N],inv[N],s[N];
bool fl[N];
inline int read()
{
	char ch=getchar();
	int ret=0;
	while(ch<'0'||ch>'9') ch=getchar();
	while(ch>='0'&&ch<='9') 
		ret=(ret<<3)+(ret<<1)+ch-'0',ch=getchar();
	return ret;
}
inline void get_pri()
{
	for(int i=1;i<=n;i++) fl[i]=0;
	for(int i=2;i<=n;i++)
	{
		if(!fl[i]) pri[++cnt]=i;
		for(int j=1;j<=cnt;j++)
		{
			if(i*pri[j]>n) break;
			fl[i*pri[j]]=1;
			if(i%pri[j]==0)	break;
		}
	}
}
inline void get_jc()
{
	jc[1]=1;
	for(int i=2;i<=n;i++) jc[i]=(ll)jc[i-1]*i%p;
}
inline void get_inv()
{
	inv[1]=1;
	for(int i=2;i<=n;i++) 
		inv[i]=1LL*(p-p/i)*inv[p%i]%p;
}
inline void get_s()
{
	s[1]=1;
	for(int i=2;i<=n;i++)
		s[i]=(fl[i])?s[i-1]:1LL*s[i-1]*(i-1)%p*inv[i]%p;
}
int main()
{
	int T=read(); p=read();
	get_pri(); get_jc(); get_inv(); get_s();
	while(T--)
	{
		int a=read(),b=read();
		printf("%d\n",(ll)jc[a]*s[b]%p);
	}
	return 0;
}