csuoj1163寒衣调(乘法逆元)

Description

男从戎,女守家。一夜,狼烟四起,男战死沙场。从此一道黄泉,两地离别。最后,女终于在等待中老去逝去。逝去的最后是换尽一生等到的相逢和团圆。
某日两人至奈何桥前,服下孟婆汤。
每滴孟婆汤都有强度不一的药效,设一碗孟婆汤共N滴(0<N<100000),其中第i滴(0≤i<N)用b[i]表示。
孟婆汤的药效与原料有关,设熬制前同样有N滴原料,第i滴原料用a[i]表示,0≤a[i]<2^32。
药效b[i]的计算方法为b[i]=(a[0]*a[1]*...*a[N-1]/a[i])%m(假设0/0=1),0<b[i]<2^32。

Input

每行开头给出原料数量N,取模数m,紧接着的一行按顺序给出原料a[i]。求出熬制所成孟婆汤的药效b[i],每次输完一碗孟婆汤的药效后以换行结尾。

Output

 求出熬制所成孟婆汤的药效b[i],每碗孟婆汤后以换行结尾。

Sample Input

5 11
2 7 5 3 9
3 7
9 8 5

Sample Output

10 6 4 3 1
5 3 2

如果不用乘法逆元,数字爆炸
乘法逆元是什么?emmm  推荐一篇博客  http://www.cnblogs.com/dupengcheng/p/5487362.html  不过他的求法是用扩展欧几里得,这里直接用他的结论逆元x=a的p-2次方
代码
#include<iostream>
#include<cstdio>
#define ll long long
int a[100005];
ll b[100005];
ll pow(int a,int n,int p)//快速幂
{
  ll x=a;
  ll res=1;
  while(n)
  {
    if(n&1)
       res=((ll)res*(ll)x)%p;
    n>>=1;
    x=((ll)x*(ll)x)%p;
  }
  return res;
}
int main()
{
	int N,M;
	while (~scanf("%d%d",&N,&M))
	{
    long long sum=1;
		for(int i=0;i<N;i++)
    {
      scanf("%d",&a[i]);
      if(a[i]==0)
        continue;
      b[i]=pow(a[i],M-2,M);
      sum=(sum*a[i])%M;
    }
    for(int i=0;i<N;i++)
    printf("%lld ",a[i]==0?1:(sum*b[i])%M);
    printf("\n");
	}
	return 0;
}





评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值