分析
对于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);
}