Ned 的难题

题目

分析

对于20%，\[ans=\Pi_{i=1}^{n}\Pi_{j=i}^{n}gcd(a_{i},a_{i+1},...,a_{j-1},a_{j})\]
显然这是会超时的，那么我们换个方法，
假设当前做到\(i\)，设\(b_{j}表示gcd(a_{j},a_{j+1},...,a_{i-1})\)
那么
\[ans=ans*a_{i}*\Pi_{j=1}^{i-1}gcd(b_{j},a_{i})\]
发现\(b\)一定是单调递增的，而且，其中有很多的\(b_{j}\)是相同的。接着，每次递增都至少是翻倍的。
那么，就可以加个链表。
时间复杂度\(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);
}