题目大意:
给你一个序列A,f(l,r) 表示 在[l,r]中 的Ai 对于每一个数Aj 都有 Ai%Aj!=0 的数目( i!=j )
卡了一段时间。。。。。
题解
简单题 定义两个数组L[i],R[i],表示第i数左侧和右侧最接近它且值是A[i]因子的数的位置,那么第i个数贡献的答案就是(R[i]-i ) * (i-L[i])
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <ctime>
#include <algorithm>
#include <iostream>
#include <sstream>
#include <string>
#define oo 0x13131313
const long long MOD=(1e9)+7;
using namespace std;
long long A[100005];
long long Q[100005];
long long H[100005];
long long C[100005];
long long P[100005];
long long pos[10005];
long long n;
void input()
{
memset(Q,0,sizeof(Q));
memset(H,0,sizeof(H));
memset(A,0,sizeof(A));
for(int i=1;i<=n;i++)
scanf("%d",&A[i]);
}
void YCL()
{
for(int i=1;i<=100000;i++)
{
P[i]=(P[i-1]+i)%MOD;
C[i]=(C[i-1]+P[i])%MOD;
}
}
void step1()
{
memset(pos,0,sizeof(pos));
for(int i=1;i<=n;i++)
{
long long MAX=0;
long long t=sqrt(A[i]);
for(int j=1;j<=t;j++)
{
if(A[i]%j==0)
{
MAX=max(MAX,pos[j]);
MAX=max(MAX,pos[A[i]/j]);
}
}
Q[i]=MAX;
pos[A[i]]=i;
}
for(int i=0;i<=10000;i++)
pos[i]=n+1;
for(int i=n;i>=1;i--)
{
long long MIN=n+1;
long long t=sqrt(A[i]);
for(int j=1;j<=t;j++)
{
if(A[i]%j==0)
{
MIN=min(MIN,pos[j]);
MIN=min(MIN,pos[A[i]/j]);
}
}
H[i]=MIN;
pos[A[i]]=i;
}
}
void solve()
{
long long ans=0;
for(int i=1;i<=n;i++)
{
long long len1=H[i]-Q[i];
long long len2=i-Q[i];
long long len3=H[i]-i;
ans=(ans+(len1*(len1+1)/2)-(len2*(len2+1)/2)-(len3*(len3+1)/2))%MOD;
}
cout<<ans<<endl;
}
int main()
{
// freopen("a.in","r",stdin);
YCL();
while(cin>>n)
{
input();
step1();
solve();
}
}