Description
1<=n<=50000
1<=a[i]<=10^7
Analysis
突然发现漏写了一题==
其实不难==
设
b[j]=gcd(a[j],a[j+1]...a[i])(j≤i)
往后一位,对于i+1的答案,为
a[i+1]∗∏ij=1gcd(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;
}