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;
}