zoukankan      html  css  js  c++  java
  • 济南学习 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 }

    分段打表~~

  • 相关阅读:
    (8)FastDFS
    (7)文件上传
    (6)品牌新增
    数据仓库_MySQL(2)
    数据仓库_Linux(5)&MySQL(1)
    J哥说生产事故之僵尸进程
    J哥说生产事故之CPU爆表
    idea classpath
    (五)返回两个数组之间的差异
    (四)数组扁平化
  • 原文地址:https://www.cnblogs.com/suishiguang/p/6037695.html
Copyright © 2011-2022 走看看