求逆元~

最新推荐文章于 2023-09-21 21:20:20 发布

lentty

最新推荐文章于 2023-09-21 21:20:20 发布

阅读量619

点赞数

分类专栏：大学ACM记忆

本文链接：https://blog.csdn.net/lentty1452/article/details/8685082

版权

大学ACM记忆专栏收录该内容

85 篇文章 0 订阅

订阅专栏

Problem G: 数学

Time Limit: 1 Sec Memory Limit: 128 MB
Submit: 707 Solved: 103

Submit

Status

Web Board

Description

给出数列A1,A2 ,...,AN，并设

现要求把所有的Bi 算出来。

Input

输入包含多组测试数据。对于每组数据，第 1 行，1 个整数N(1 ≤N≤100,000), 表示数列的长度。第 2 行，N 个整数A1,A2,...,AN(1 ≤Ai≤109)，表示给出的数列。输入以一个 0 表示结尾。

Output

对于每组数据，输出一行，N 个整数用空格分隔，表示算出的B1,B2,...,BN。

Sample Input

31 2 30

Sample Output

6 3 2

HINT

用扩展欧几里得求逆元，没啥好说的~

#include<stdio.h>
#include<stdlib.h>
#include<memory.h>
#define maxn 100100
#define mod 1000000007
typedef long long ll;
int a[maxn];
ll d,x,y;
void extend_euclid(ll a,ll b)
{
    if(b==0)
    {
    d=a;x=1;y=0;
    }
    else
    {
        extend_euclid(b,a%b);
        ll tmp=x;
        x=y;
        y=tmp-a/b*y;
    }
}

ll inverse(ll a,ll n)
{
    extend_euclid(a,n);
    if(d==1)
    return (x%n+n)%n;
    else
    return -1;
}
int main()
{
    freopen("in.txt","r",stdin);
    int n;
    while(scanf("%d",&n)==1&&n)
    {
        int i;
        ll sum=1;
        for(i=0;i<n;i++)
        {
          scanf("%d",&a[i]);
          sum=(sum*a[i])%mod;
        }
        for(i=0;i<n;i++)
        {
            if(i>0) printf(" ");
            ll ans=inverse(a[i],mod);
            printf("%lld",(sum*ans)%mod);
        }
        printf("\n");
    }
    return 0;
}