NP(np)
Time Limit:1000ms Memory Limit:64MB
题目描述
LYK 喜欢研究一些比较困难的问题,比如 np 问题。
这次它又遇到一个棘手的 np 问题。问题是这个样子的:有两个数 n 和 p,求 n 的阶乘
对 p 取模后的结果。
LYK 觉得所有 np 问题都是没有多项式复杂度的算法的, 所以它打算求助即将要参加 noip
的你,帮帮 LYK 吧!
输入格式(np.in)
输入一行两个整数 n,p。
输出格式(np.out)
输出一行一个整数表示答案。
输入样例
3 4
输出样例
2
数据范围
对于 20%的数据:n,p<=5。
对于 40%的数据:n,p<=1000。
对于 60%的数据:n,p<=10000000。
对于 80%的数据:n<=10^18,p<=10000000。
对于另外 20%的数据:n<=10^18,p=1000000007。
其中大致有 50%的数据满足 n>=p。
#include<cstdio>
#define name "np"
#define ll long long
using namespace std;
const int a[]={,,,,,,
,,,,,,,
,,,,,,,
,,,,,,,
,,,,,,,
,,,,,,,
,,,,,,,
,,,,,,,,
,,,,,,,,
,,,,,,,,
,,,,,,,,
,,,,,,,,
,,,,,,,,
,,,,};
const int mod=1e9+;
const int sz=1e7;
ll n,p;
ll ans=;
int main(){
scanf("%I64d%I64d",&n,&p);
if(n>=p){puts("");return ;}// 比较好想的一个暴力,但确实有用 并且用处很大
if(p==mod){// 分段打表
n<sz?ans=:ans=a[n/sz];
for(ll i=n/sz*sz+;i<=n;i++) ans=ans*i%p;
printf("%I64d",ans);
return ;
}
for(ll i=;i<=n;i++) ans=ans*i%p;
printf("%I64d",ans);
return ;
}
分段打表~~