hdu 5288 OO’s Sequence(15多校第一场1001)

最新推荐文章于 2017-06-26 12:22:00 发布

原创最新推荐文章于 2017-06-26 12:22:00 发布 · 381 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#hdu5288 #OOs Sequence #2015多校

15多校专栏收录该内容

17 篇文章

订阅专栏

本文介绍了一种针对小范围因子查询的高效算法——暴力求解因子的方法。通过预处理和优化查找过程，实现快速计算特定数值范围内因子的数量，适用于数的范围较小的场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

//15多校第一场 1001 暴力
//若有多数需要求因子，且数的范围较小，可以考虑预处理
#include<cstdio>
#include<algorithm>
#include<vector>
using namespace std;
const int N=100100;
int a[N];
vector<int>f[10010],vec[10010];
const long long M=1000000007;
int n;
#define MAX(a,b) a<b?b:a
#define MIN(a,b) a<b?a:b
typedef long long LL;
void init(){
    for(int i=1;i<=10000;i++)
        for(int j=1;j*i<=10000;j++)
            f[i*j].push_back(i);
}
LL find(int pos,int ai){
    int l=0,r=n+1;
    for(int i=0;i<f[ai].size();i++){
        int fac=f[ai][i];
        for(int j=0;j<vec[fac].size();j++){
            if(vec[fac][j]<pos){
                if(l<vec[fac][j])l=vec[fac][j];
            }
            if(vec[fac][j]>pos){
                if(vec[fac][j]<r)r=vec[fac][j];
            }
        }
        
    }
    return (LL)(r-pos)*(pos-l)%M;
}
int main(){
    init();
    while(scanf("%d",&n)!=EOF){
        for(int i=0;i<=10000;i++)
            vec[i].clear();//从1到10000清零
        for(int i=1;i<=n;i++){
            scanf("%d",&a[i]);
            vec[a[i]].push_back(i);//每一个ai的位置
        }
        long long sum=0;
        for(int i=1;i<=n;i++){
            sum=(sum+find(i,a[i])%M)%M;
        }
        printf("%lld\n",sum);
    }
    return 0;
}