问题描述:

林记在做数学习题的时候,经常遇到这种情况:苦思冥想了很久终于把问题解出来,结果发现答案是0,久而久之林记在得到习题答案是0的时候就没有了做出一道难题的成就感。于是林记决定:以后出题,答案一定不能是0,例如求n!最低位非零数这样的习题就很不错了。

现在林记提出了一个更难一点的问题:求n!在K进制下的最低位非零数。其中K符合一些特殊的条件:K是由若干个互不相同的质数相乘得出来的,例如K=2,3,5,6,7,10……

输入格式:

首先输入的第一行是一个整数Q,表示询问的个数。

接下来是Q个询问,每个询问有两行组成,第一行首先是一个整数nk,然后紧跟着nk个正整数:m1,m2,……mnk,则K为所有mi的乘积,输入保证mi是有若干个互不相同的质数相乘得出来的。接下来一行给出一个整数n。

输出格式:

对于每个询问,如果K不符合题目限制,则输出“I hate zero.”(不用加双引号)。否则输出K进制下n!的最低位非零数。

输入样例:

输出样例:

3

1 17

9

1 6

49

1 93

16

15

2

78

数据范围:

对于10%的数据满足:n≤10,K≤100

对于50%的数据满足:n≤100000,K≤100

对于100%的数据满足:n≤1015,K≤1015,mi≤106,Q≤20

【分析】

  快速求阶乘大法!!

  下面某步还用到 威尔逊定理 (p-1)!%p=p-1...... (也可以预处理出来的)

  把m分解质因数,对于质数p单独求n!=p^a+b的a的值和b%p的值(logn递归搞定),然后用中国剩余定理合并起来。

  

  long long的乘积会爆,所以...要用那个不断乘2的方法:

  

LL mul(LL x,LL y,LL p)
{
if(x==0||y==0) return 0;
LL now=x,yy=y,add=0;
while(yy>1)
{
if(yy%2!=0) add=(add+now)%p;
now=(now*2)%p;
yy/=2;
}
return (now+add)%p;
}

代码如下:

 #include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cmath>
using namespace std;
#define Maxn 1000010
#define INF 0xfffffff
#define LL long long LL nk,m[Maxn],ml,K;
LL n;
LL prime[Maxn],pl;
bool q[Maxn];
LL mn; LL mymin(LL x,LL y) {return x<y?x:y;} void get_pri(LL mx)
{
pl=;
memset(q,,sizeof(q));
for(LL i=;i<=mx;i++)
{
if(q[i]) prime[++pl]=i;
for(int j=;j<=pl;j++)
{
if(i*prime[j]>mx) break;
q[i*prime[j]]=;
if(i%prime[j]==) break;
}
}
} bool div(LL mm)
{
for(LL i=;prime[i]*prime[i]<=mm;i++) if(mm%prime[i]==)
{
m[++ml]=prime[i];
mm/=prime[i];
while(mm%prime[i]==) return ;
}
if(mm!=) m[++ml]=mm;
return ;
} bool init()
{
ml=;
scanf("%lld",&nk);
K=;
bool ok=;
for(int i=;i<=nk;i++)
{
LL mm;
scanf("%lld",&mm);K*=mm;
if(!div(mm)) ok=;
}
sort(m+,m++ml);
for(LL i=;i<=ml;i++) if(m[i]==m[i-]) {ok=;break;}
scanf("%lld",&n);
if(!ok) return ;
return ;
} LL f[Maxn]; struct node{LL x,y;};
node c[Maxn]; node ffind(LL x,LL p)
{
node ans;
if(x<p)
{
ans.x=;
ans.y=f[x];
return ans;
}
ans.x=;ans.y=;
node t=ffind(x/p,p);
ans.x+=t.x; ans.y*=t.y;
ans.y=(ans.y*f[x%p])%p;
ans.x+=x/p; LL y=(((x/p)%)==)?:p-;
ans.y=(ans.y*y)%p; return ans;
} LL mul(LL x,LL y,LL p)
{
if(x==||y==) return ;
LL now=x,yy=y,add=;
while(yy>)
{
if(yy%!=) add=(add+now)%p;
now=(now*)%p;
yy/=;
}
return (now+add)%p;
} LL qpow(LL x,LL b,LL p)
{
LL ans=;
while(b)
{
if(b&) ans=mul(ans,x,p);
x=mul(x,x,p);
b>>=;
}
return ans;
} LL d[Maxn]; LL get_ans()
{
mn=c[].x;
for(LL i=;i<=ml;i++) mn=mymin(mn,c[i].x);
for(LL i=;i<=ml;i++)
{
if(c[i].x>mn) {d[i]=;continue;}
LL bm=qpow(K/m[i],mn,m[i]);
bm=qpow(bm,m[i]-,m[i]);
d[i]=mul(bm,c[i].y,m[i]);//(bm*c[i].y)%m[i];
}
LL ans=;
for(LL i=;i<=ml;i++)
{
LL y=qpow(K/m[i],m[i]-,m[i]);
//ans=(ans+((K/m[i]*y)%K)*d[i])%K;
ans=(ans+mul(mul(K/m[i],y,K),d[i],K))%K;
}
return ans;
} int main()
{
get_pri();
int T;
scanf("%d",&T);
while(T--)
{
if(!init()) {printf("I hate zero.\n");continue;}
for(LL i=;i<=ml;i++)
{
f[]=;
for(LL j=;j<m[i];j++) f[j]=(f[j-]*j)%m[i];
c[i]=ffind(n,(LL)m[i]);
}
printf("%lld\n",get_ans());
}
return ;
}

2016-09-04 19:37:26

【GDOI 2011 DAY2 T3】零什么的最讨厌了 (快速求阶乘、中国剩余定理)的更多相关文章

  1. 【NOIP 2017】Day2 T3 列队

    Problem Description \(Sylvia\) 是一个热爱学习的女孩子. 前段时间,\(Sylvia\) 参加了学校的军训.众所周知,军训的时候需要站方阵. \(Sylvia\) 所在的 ...

  2. 数论-质数 poj2689,阶乘分解,求阶乘的尾零hdu1124, 求尾零为x的最小阶乘

    /* 要求出[1,R]之间的质数会超时,但是要判断[L,R]之间的数是否是素数却不用筛到R 因为要一个合数n的最大质因子不会超过sqrt(n) 所以只要将[2,sqrt(R)]之间的素数筛出来,再用这 ...

  3. 【BZOJ 4518】【SDOI 2016 Round1 Day2 T3】征途

    比较明显的斜率优化DP,省选时因为时间太紧张和斜率DP写得不熟等原因只写了60分的暴力DP,其实当时完全可以对拍来检验标算的正确,但是我当时too naive- 很快打完了,调了将近一晚上QAQ,因为 ...

  4. 【NOIP 2013 DAY2 T3】 华容道(spfa)

    题目描述 [问题描述] 小 B 最近迷上了华容道,可是他总是要花很长的时间才能完成一次.于是,他想到用编程来完成华容道:给定一种局面, 华容道是否根本就无法完成,如果能完成, 最少需要多少时间. 小 ...

  5. 【NOIP 2015 DAY2 T3】 运输计划 (树链剖分-LCA)

    题目背景 公元 2044 年,人类进入了宇宙纪元. 题目描述 L 国有 n 个星球,还有 n-1 条双向航道,每条航道建立在两个星球之间,这 n-1 条航道连通了 L 国的所有星球. 小 P 掌管一家 ...

  6. 【NOIP2015提高组】 Day2 T3 运输计划

    题目描述 L 国有 n 个星球,还有 n-1 条双向航道,每条航道建立在两个星球之间,这 n-1 条航道连通了 L 国的所有星球. 小 P 掌管一家物流公司,该公司有很多个运输计划,每个运输计划形如: ...

  7. 一道搜索题【2013 noip提高组 DAY2 t3】华容道

    这篇不多说,具体的解释都在程序里 题目描述 [问题描述] 小 B 最近迷上了华容道,可是他总是要花很长的时间才能完成一次.于是,他想到用编程来完成华容道:给定一种局面, 华容道是否根本就无法完成,如果 ...

  8. 【NOIP 2015】Day2 T3 运输计划

    Problem Background 公元 \(2044\) 年,人类进入了宇宙纪元. Description 公元\(2044\) 年,人类进入了宇宙纪元. $L $国有 \(n\) 个星球,还有 ...

  9. 【NOIP 2016】Day2 T3 愤怒的小鸟

    Problem Description \(Kiana\) 最近沉迷于一款神奇的游戏无法自拔. 简单来说,这款游戏是在一个平面上进行的. 有一架弹弓位于 \((0,0)\) 处,每次 \(Kiana\ ...

随机推荐

  1. 使用ROW_NUMBER进行的快速分页

    DECLARE @pageSize INT ; DECLARE @pageIndex INT ; SET @pageSize = 5 SET @pageIndex =2 ; --第二页,每页显示5条数 ...

  2. 经典SQL语句大全(绝对的经典)

    ”,start为起始位置,length为字符串长度,实际应用中以len(expression)取得其长度3,right(char_expr,int_expr) 返回字符串右边第int_expr个字符, ...

  3. primefaces 带参数的组件

    分析了一下primefaces extension showcase的源码,总结一下组件传参的方式. 下面定义了一个组件,名为fourTabsDecoratorCustom.xhtml 1: < ...

  4. 不同浏览器下css 透明度的写法

    filter:alpha(opacity=90); /* IE transparent */ -moz-opacity:0.9; /* Moz + FF transparent */ opacity: ...

  5. python-MySQLdb-练习

    看完视频,自己练习一遍. 还是遇到问题,不过最终还是解决了.贴上完成的代码. CREATE TABLE `NewTable` ( `acctid` int(11) NOT NULL AUTO_INCR ...

  6. 03_MySQL中文乱码处理_02_确保MySQL插入数据不乱码的解决方法

    [要统一MySQL数据库客户端 和 服务端的字符集] 如下,就是不统一的情况: 通常,MySQL数据库的下面几个字符集(客户端 和 服务端)统一成一个字符集(即" show variable ...

  7. How to change Visual Studio default environment setting

    如何改变 Visual Studio 的默认环境设置: 1. 工具栏 Tools --> Import and Export Settings... 2. 选择 Reset All Settin ...

  8. 九度OJ 1352 和为S的两个数字

    题目地址:http://ac.jobdu.com/problem.php?pid=1352 题目描述: 输入一个递增排序的数组和一个数字S,在数组中查找两个数,是的他们的和正好是S,如果有多对数字的和 ...

  9. 百度UEditor(富文本编辑器)的基础用法

    百度的这个编辑器挺强大的,这里只是用他的文本功能,没有介绍上传图片视频的. 我用是的SSH来写的项目. 1. 把下载的UEditor(ueditor1_4_3_1-utf8-jsp)解压后全部复制到W ...

  10. IOS 学习笔记 2015-04-15 手势密码(原)

    // // WPSignPasswordView.h // 网投网 // // Created by wangtouwang on 15/4/9. // Copyright (c) 2015年 wan ...