P3467【2015多校联训1】统计

时间限制 : 20000 MS 空间限制 : 165536 KB

问题描述

现在有一个数组 a，数组中有 n 个元素。定义一个函数 f(l,r)表示 i（l<=i<=r）的 个数,其中 i 符合条件：不存在 j (l<=j<=r 且 j≠i)满足 ai mod aj = 0 求
∑i=1n∑j=inf(i,j) mod （109+7）.
即所有区间中包含的符合条件的 i 的个数。

输入格式

第一行一个整数 n（n<=100000）
第二行有 n 个数，表示数组中的元素 ai，0

题解

详情见训练赛七第二题
ps：快速求1~n的因子：
for(i=1;i<=maxn-2;i++)
{
for(j=i;j<=maxn-2;j+=i)
{
m[j][0]++;
m[j][m[j][0]]=i;
}
}

代码

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;
#define ll long long
#define maxn 10002
#define maxn1 100005
#define inf 1e9
const int mod=1e9+7;
ll m[maxn][105],a[maxn1],ls[maxn1],rs[maxn1],Last[maxn1];
ll n;
int main()
{
    int i,j,k;
    scanf("%lld",&n);
    for(i=1;i<=n;i++) scanf("%lld",&a[i]);
    for(i=1;i<=maxn-2;i++)
    {
        for(j=i;j<=maxn-2;j+=i)
        {
            m[j][0]++;
            m[j][m[j][0]]=i;
        }
    }
    for(i=1;i<=n;i++)
    {
        ll s=-inf;
        for(j=1;j<=m[a[i]][0];j++) s=max(s,Last[m[a[i]][j]]);
        ls[i]=s;
        Last[a[i]]=i; 
    }
    for(i=1;i<=100000;i++) Last[i]=n+1;
    for(i=n;i>=1;i--)
    {
        ll s=inf;
        for(j=1;j<=m[a[i]][0];j++) s=min(s,Last[m[a[i]][j]]);
        rs[i]=s;
        Last[a[i]]=i; 
    }
    ll ans=0;
    for(i=1;i<=n;i++) ans=(ans+(i-ls[i])*(rs[i]-i)%mod)%mod;
    cout<<ans;
}