济南学习 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。

 1 #include<cstdio>
 2 #define name "np"
 3 #define ll long long
 4 using namespace std;
 5 const int a[101]={0,682498929,491101308,76479948,723816384,67347853,
 6     27368307,625544428,199888908,888050723,927880474,281863274,661224977,
 7     623534362,970055531,261384175,195888993,66404266,547665832,109838563,
 8     933245637,724691727,368925948,268838846,136026497,112390913,135498044,
 9     217544623,419363534,500780548,668123525,128487469,30977140,522049725,
10     309058615,386027524,189239124,148528617,940567523,917084264,429277690,
11     996164327,358655417,568392357,780072518,462639908,275105629,909210595,
12     99199382,703397904,733333339,97830135,608823837,256141983,141827977,696628828,
13     637939935,811575797,848924691,131772368,724464507,272814771,326159309,456152084,
14     903466878,92255682,769795511,373745190,606241871,825871994,957939114,435887178,
15     852304035,663307737,375297772,217598709,624148346,671734977,624500515,748510389,
16     203191898,423951674,629786193,672850561,814362881,823845496,116667533,256473217,
17     627655552,245795606,586445753,172114298,193781724,778983779,83868974,315103615,
18     965785236,492741665,377329025,847549272,698611116};
19 const int mod=1e9+7;
20 const int sz=1e7;
21 ll n,p;
22 ll ans=1;
23 int main(){
24     scanf("%I64d%I64d",&n,&p);
25     if(n>=p){puts("0");return 0;}// 比较好想的一个暴力，但确实有用 并且用处很大 
26     if(p==mod){//    分段打表 
27         n<sz?ans=1:ans=a[n/sz];
28         for(ll i=n/sz*sz+1;i<=n;i++) ans=ans*i%p;
29         printf("%I64d",ans);
30         return 0;
31     }
32     for(ll i=1;i<=n;i++) ans=ans*i%p;
33     printf("%I64d",ans);
34     return 0;
35 }

分段打表~~