题目描述

大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞票。房地产第一大户沙拉公主决定预测一下大富翁国现在所有真钞票的数量。现在,请你帮助沙拉公主解决这个问题,由于可能张数非常大,你只需计算出对R取模后的答案即可。R是一个质数。//codevs这里有坑,R是合数

输入输出格式

输入格式:

第一行为两个整数T,R。R<=10^9+10,T<=10000,表示该组中测试数据数目,R为模 后面T行,每行一对整数N,M,见题目描述 m<=n

输出格式:

共T行,对于每一对N,M,输出1至N!中与M!素质的数的数量对R取模后的值

输入输出样例

输入样例#1:

1 11
4 2
输出样例#1:

1

数据范围:
对于100%的数据,1 < = N , M < = 10000000

解题思路

  http://www.cnblogs.com/yangyaojia/p/6434611.html

  http://blog.csdn.net/loi_dqs/article/details/50520309

吐槽

  这题一题能抵三题做了。

  不想被卡常就离线吧,筛法筛到最大的n即可,否则卡常、卡内存太严重了,我的好不容易洛谷ac,跑了6s。

  求逆元可以线性递推,也可以扩欧。

  popoqqq大神博客的代码(https://www.luogu.org/record/show?rid=2565704)、黄学长博客(https://www.luogu.org/record/show?rid=2565713)里线性递推逆元的代码都在洛谷爆零,不知为啥。

  我的代码离线int会溢出,估计是筛法出问题,还不够优化,于是递推求逆元的数组开不下了,就用扩欧吧。

  卡一波常1.3s跑完(BZOJ不开O2跑了2.9s)。

源代码

6s的,求逆元用线性递推(BZOJ TLE)

#include<cstdio>
#include<cstring>
#include<algorithm>
int T,R,n,m;
long long jc[]={};//阶乘取模
long long inv[]={};//逆元
bool not_prime[]={};//筛子
long long ans[]={};//pi累乘 void init()
{
inv[]=;
for(int i=;i<=&&i<R;i++)
inv[i]=(R-R/i)%R*inv[R%i]%R;
ans[]=;jc[]=;
for(int i=;i<=;i++)
{
jc[i]=jc[i-]*i%R;
if(!not_prime[i])
{
ans[i]=ans[i-]*(i-)%R*inv[i%R]%R;
for(int j=i<<;j<=;j+=i)
not_prime[j]=;
}
else
ans[i]=ans[i-];
}
} int main()
{
scanf("%d%d",&T,&R);
init();
while(T--)
{
scanf("%d%d",&n,&m);
printf("%lld\n",(long long)(jc[n])*(long long)(ans[m])%R);
}
return ;
}

1.3s的,求逆元用扩欧

#include<cstdio>
#include<algorithm>
int T,R,maxn=-;
long long jc[]={};//阶乘取模
//int inv[10000010]={0};//逆元
bool not_prime[]={};//筛子
long long ans[]={};//pi累乘
int n[]={},m[]={}; void exgcd(int a,int b,int &x,int &y)
{
if(b==){x=;y=;return;}
exgcd(b,a%b,x,y);
int t=x;x=y;y=t-a/b*y;
}
inline int getine(int t)
{
int x,y;
exgcd(t,R,x,y);
return (x%R+R)%R;
}//扩欧求逆元 inline void Read(int &in){
static char ch;
for(ch=getchar();ch>''||ch<'';ch=getchar()) ;
for(in=;ch>=''&&ch<='';ch=getchar()) in=in*+ch-'';
}
void init()
{
/*inv[1]=1;
for(int i=2;i<=maxn&&i<R;i++)
inv[i]=(R-R/i)%R*inv[R%i]%R;*///线性递推求逆元
ans[]=;jc[]=;
for(int i=;i<=maxn;i++)
{
jc[i]=jc[i-]*i%R;
if(!not_prime[i])
{
ans[i]=ans[i-]*(i-)%R*getine(i%R)%R;
for(int j=i<<;j<=maxn;j+=i)
not_prime[j]=;
}
else
ans[i]=ans[i-];
}
} int main()
{
Read(T);Read(R);
for(register int i=;i<=T;i++)
Read(n[i]),Read(m[i]),maxn=std::max(maxn,n[i]);
init();
for(register int i=;i<=T;i++)
printf("%lld\n",jc[n[i]]*ans[m[i]]%R);
return ;
}

洛谷记录里最优的代码,用户ID小小莫。Orz,洛谷跑了0.9s,BZOJ不开O2跑了2.6s,素数筛法值得学习。

#include <cstdio>
#include <algorithm>
using namespace std;
static const int maxm = + ;
static const int maxx = + ;
typedef long long LL;
bool Notprime[maxx];
int inv[maxx],phi[maxx],fac[maxx],prime[maxx];
int T,n[maxm],m[maxm],p,cnt,tot; void Read(int &); inline void Linear_shaker(){
int maxn = tot+;
fac[]=;inv[]=;phi[]=;
for(register int i=;i<maxn;i++){
if(!Notprime[i])
prime[++cnt]=i;
for(register int j=;j<=cnt && i*prime[j]<maxn;j++){
Notprime[i*prime[j]]=;
if(i%prime[j]==) break;
}
}
for(register int i=;i<=maxn;i++){
fac[i] = ((LL)fac[i-]*i)%p;
if(i<p) inv[i] = (LL)(p-p/i)*inv[p%i]%p;
if(!Notprime[i]) phi[i] = (LL)phi[i-]*(i-)%p*inv[i%p]%p;
else phi[i] = phi[i-];
}
}
int main(){
Read(T),Read(p);
for(register int i=;i<=T;i++) Read(n[i]),Read(m[i]),tot = max(tot,n[i]);
Linear_shaker();
for(register int i=;i<=T;i++)
printf("%lld\n",(LL)fac[n[i]]*phi[m[i]]%p);
return ;
} void Read(int &in){
static char ch;
for(ch=getchar();ch>''||ch<'';ch=getchar()) ;
for(in=;ch>=''&&ch<='';ch=getchar()) in=in*+ch-'';
}

洛谷 P2155 BZOJ 2186 codevs 2301 [SDOI2008]沙拉公主的困惑的更多相关文章

  1. 洛谷 P2155 [SDOI2008]沙拉公主的困惑 解题报告

    P2155 [SDOI2008]沙拉公主的困惑 题目描述 大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为\(1\)到\(N\)的阶乘,但是,政府只发行编号与\(M!\ ...

  2. 洛咕 P2155 [SDOI2008]沙拉公主的困惑

    洛咕 P2155 [SDOI2008]沙拉公主的困惑 有个结论,就是如果\(gcd(a,b)=1\),那么\(gcd(a+kb,b)=1\).证明比较显然. 所以这个题目要问的\(n!\)就可以分成\ ...

  3. Bzoj 2186: [Sdoi2008]沙拉公主的困惑 乘法逆元,线性筛,欧拉函数,数论

    2186: [Sdoi2008]沙拉公主的困惑 Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 2560  Solved: 857[Submit][St ...

  4. 数学(逆元):BZOJ 2186: [Sdoi2008]沙拉公主的困惑

    2186: [Sdoi2008]沙拉公主的困惑 Description 大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞 ...

  5. 【BZOJ 2186】 2186: [Sdoi2008]沙拉公主的困惑 (欧拉筛,线性求逆元)

    2186: [Sdoi2008]沙拉公主的困惑 Description 大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M!互质的钞 ...

  6. 【bzoj2186】[Sdoi2008]沙拉公主的困惑

    2186: [Sdoi2008]沙拉公主的困惑 Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 3303  Solved: 1129[Submit][S ...

  7. BZOJ2186: [Sdoi2008]沙拉公主的困惑(求[1,N!]与M!互素的个数)(线性筛)

    2186: [Sdoi2008]沙拉公主的困惑 Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 6103  Solved: 2060[Submit][S ...

  8. BZOJ2186 [Sdoi2008]沙拉公主的困惑 【数论,欧拉函数,线性筛,乘法逆元】

    2186: [Sdoi2008]沙拉公主的困惑 Time Limit: 10 Sec  Memory Limit: 259 MB Submit: 5003  Solved: 1725 [Submit] ...

  9. 【BZOJ2186】[Sdoi2008]沙拉公主的困惑 线性筛素数

    [BZOJ2186][Sdoi2008]沙拉公主的困惑 Description 大富翁国因为通货膨胀,以及假钞泛滥,政府决定推出一项新的政策:现有钞票编号范围为1到N的阶乘,但是,政府只发行编号与M! ...

随机推荐

  1. KVO中你所不知道的"坑"

      一.什么是 KVO 首先让我们了解一下什么KVO,全称为Key-Value Observing,是iOS中的一种设计模式,用于检测对象的某些属性的实时变化情况并作出响应.键值观察Key-Value ...

  2. [Swift通天遁地]七、数据与安全-(5)高效操作SQLite数据库

    ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★➤微信公众号:山青咏芝(shanqingyongzhi)➤博客园地址:山青咏芝(https://www.cnblogs. ...

  3. Boost.Build特点(译)

    Boost.Build Boost.Build makes it easy to build C++ projects, everywhere. Boost.Build让构建C++项目在任何地方都很容 ...

  4. 安卓5.0新特性之Palette

    根据图片来决定标题的颜色和标题栏的背景色,这样视觉上更具有冲击力和新鲜感,而不像统一色调那样呆板. Palette这个类能提取以下突出的颜色: Vibrant(充满活力的) Vibrant dark( ...

  5. BZOJ 3779 LCT 线段树 DFS序 坑

    hhhh抄了半天lty代码最后T了  对拍也没事.. 药丸 mine #pragma GCC optimize("O3") //By SiriusRen #include < ...

  6. cmd bat 相对命令

    "%~dp0",在BAT中,是不是“相对路径”的意思 (2013-08-21 12:19:32) 转载▼ 标签: 杂谈 分类: C# 0念 零 ,代表你的批处理本身. d p是FO ...

  7. WPF PasswordBox MVVM 实现

    由于PasswordBox.Password属性非依赖属性,所以不能作为绑定的目标,以下是本人的MVVM实现方法. PasswordBox.Password与TextBox.Text同步,TextBo ...

  8. 再战primer——decltype 和引用

    刷primer看到原文讲到"引用从来都作为其所指对象的同义词出现,只有用在decltype处是一个例外.",我很是好奇. 这个“引用”当然是指引用类型,like this: ; i ...

  9. Java系列学习(三)-基础语法

    1.关键字 特点:全部小写 2.标识符 (1)就是给类,接口,方法等起名字的字符序列 (2)组成规则: A:英文大小写字母 B:数字 C:$和_ (3)注意事项: A:不能以数字开头 B:不能是jav ...

  10. Android 使用MySQL直接访问数据库

    在实际项目中,一般很少直接访问MySQL数据库,一般情况下会通过http请求将数据传送到服务端,然后在服务端连接mysql数据库. 在android 中,会通过使用Jdbc 连接MySQL 服务器 p ...