Ned 的难题

最新推荐文章于 2023-01-13 10:52:13 发布

无尽的蓝黄

最新推荐文章于 2023-01-13 10:52:13 发布

阅读量453

点赞数 1

分类专栏：题解数论

本文链接：https://blog.csdn.net/chen1352/article/details/51852561

版权

题解同时被 2 个专栏收录

177 篇文章 0 订阅

订阅专栏

数论

31 篇文章 0 订阅

订阅专栏

题目

这里写图片描述

分析

对于20%，

a n s = Π n i = 1 Π n j = i g c d (a i, a i + 1, . . ., a j - 1, a j)

$ans=\Pi_{i=1}^{n}\Pi_{j=i}^{n}gcd(a_{i},a_{i+1},...,a_{j-1},a_{j})$
显然这是会超时的，那么我们换个方法，
假设当前做到

i $i$ ，设

bj表示gcd(aj,aj+1,...,ai−1) $b_{j}表示gcd(a_{j},a_{j+1},...,a_{i-1})$
那么

a n s = a n s * a i * Π i - 1 j = 1 g c d (b j, a i)

$ans=ans*a_{i}*\Pi_{j=1}^{i-1}gcd(b_{j},a_{i})$
发现

b $b$ 一定是单调递增的，而且，其中有很多的

bj $b_{j}$ 是相同的。接着，每次递增都至少是翻倍的。
那么，就可以加个链表。
时间复杂度

O(NlogN) $O(NlogN)$ 。

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const long long maxlongint=2147483647;
const long long mo=1000000009;
const long long N=50005;
using namespace std;
long long a[N],ans=1,n,m,b[N],tot,v[N];
long long gcd(long long x,long long y)
{
    return x==0?y:gcd(y%x,x);
}
long long mi(long long x,long long y)
{
    long long sum=1;
    while(y)
    {
        if(y&1) sum=sum*x%mo;
        x=x*x%mo;
        y/=2;
    }
    return sum;
}
int main()
{
    freopen("ned.in","r",stdin);
    freopen("ned.out","w",stdout);
    scanf("%lld",&n);
    for(long long i=1;i<=n;i++)
    {
        scanf("%lld",&a[i]);
        b[i]=a[i];
        v[i]=i;
    }
    ans=a[1];
    for(long long i=2;i<=n;i++)
    {
        ans=ans*a[i]%mo;
        long long k=i;
        for(long long j=i-1;j;j=v[j]-1)
        {
            b[j]=gcd(b[j],a[i]);
            if(b[j]==b[k]) 
                v[k]=v[j];
            ans=ans*mi(b[j],j-v[j]+1)%mo;
            k=j;
        }
    }
    printf("%d",ans);
}