济南学习 Day 3 T1 am

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 ;
 }

分段打表~~