【jzoj4586】【Ned的难题】【单调栈】

题目大意

给出一个序列，求所有连续子序列gdc和。

题解

可以发现，所有质因子对答案的贡献都是独立的。考虑每个数所含当前质因子的个数，可以发现一段区间的gcd等于取min。可以用单调栈维护，保证单调不下降，栈里的每个位所代表的区间到当前区间可以形成区间，它们都可以贡献栈里所代表的区间的高度个当前质因子的答案，说不太清楚看一下code。这样一个一个搞就可以了。

观察一下本题的特殊性，扩大区间一定不会让gcd变大。设当前为第i位，f[j]表示j到i所代表的区间的gcd，显然不同的f[j]不会大于log(n)个，减小最少除以2。那我们可以用并查集维护相同的f[j]，新加入点时，之前的区间不会变多，只可能变少，暴力求gdc，用并查集合并就可以了。

code

#include<set>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
#define fo(i,j,k) for(int i=j;i<=k;i++)
#define fd(i,j,k) for(int i=j;i>=k;i--)
using namespace std;
int const maxn=50000,maxa=10000000,mod=1000000009;
int n,m,top,a[maxn+10],ss[maxn+10],stack[maxn+10],size[maxn+10],pre[maxn+10];
bool flag[maxn+10];
LL Pow(LL x,LL y){
    LL z=1;
    while(y){
        if(y&1)z=(z*x)%mod;
        x=(x*x)%mod;
        y/=2;
    }
    return z;
}
int main(){
    //freopen("ned.in","r",stdin);
    //freopen("ned.out","w",stdout);
    freopen("d.in","r",stdin);
    freopen("d.out","w",stdout);
    scanf("%d",&n);
    fo(i,1,n)
        scanf("%d",&a[i]);
    int sqa=sqrt(maxa);
    fo(i,2,sqa)
        if(!flag[i]){
            ss[++ss[0]]=i;
            fo(j,2,sqa/i)flag[i*j]=1;
        }
    LL ans=1;
    fo(i,1,ss[0]){
        long long sum=0;top=0;
        fo(j,1,n){
            int cnt=0;
            while(a[j]%ss[i]==0){
                cnt++;
                a[j]/=ss[i];
            }
            int tmp=1;
            while((top)&&(cnt<=stack[top]))
                tmp+=size[top--];
            if(cnt){
                stack[++top]=cnt;
                size[top]=tmp;
                pre[top]=pre[top-1]+cnt*tmp;
                sum+=pre[top];
            }
        }
        ans=(ans*Pow(ss[i],sum))%mod;
    }
    int i=1;
    while(i<=n){
        int j=i;while((j<n)&&(a[j+1]==a[i]))j++;
        ans=(ans*Pow(a[i],(j-i+1)*(j-i+2)/2))%mod;
        i=j+1;
    }
    printf("%lld\n",ans);
    return 0;
}