【JZOJ 4586】 Ned 的难题

Description

1<=n<=50000
1<=a[i]<=10^7

Analysis

突然发现漏写了一题==
其实不难==
设$b[j]=gcd(a[j],a[j+1]...a[i])(j\leq i)$
往后一位，对于i+1的答案，为
$a[i+1]*\prod_{j=1}^igcd(b[j],a[i+1])$
再更新b数组。
但是忽略了b的一个性质：b[i]~b[1]单调不升。
而且b[i]的取值最多只有log种，因为每次gcd一遍若其发生改变必然会除以二。
所以我们相同的b一起计算，多开一个数组记录做完当前往前跳到哪里即可。

Code

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
typedef long long ll;
const int N=50010,mo=int(1e9)+9;
int n;
ll a[N],b[N];
ll qmi(ll x,int n)
{
    ll t=1;
    for(;n;n>>=1)
    {
        if(n&1) t=t*x%mo;
        x=x*x%mo;
    }
    return t;
}
ll gcd(ll x,ll y)
{
    if(x%y==0) return y;
    else return gcd(y,x%y);
}
int main()
{
    freopen("ned.in","r",stdin);
    freopen("ned.out","w",stdout);
    ll ans=1;
    scanf("%d",&n);
    fo(i,1,n) scanf("%lld",&a[i]),ans=ans*a[i]%mo,b[i]=i;
    fo(i,1,n)
    {
        for(int j=i-1,k=i;j;k=j,j=b[j]-1)
        {
            a[j]=gcd(a[j],a[i]);
            if(a[j]==a[k]) b[k]=b[j];
            ans=ans*qmi(a[j],j-b[j]+1)%mo;
        }
    }
    printf("%lld",ans);
    return 0;
}