码迷,mamicode.com
首页 > 其他好文 > 详细

hdu 5288 OO’s Sequence(2015多校第一场第1题)枚举因子

时间:2015-07-23 01:02:36      阅读:176      评论:0      收藏:0      [点我收藏+]

标签:因子枚举

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5288


题意:在闭区间[l,r]内有一个数a[i],a[i]不能整除 除去自身以外的其他的数,f(l,r)表示在这区间内a[i]这样的数的个数,,现给你n个数,求所有区间的f(l,r)的和。


思路:对于每个数a[i]求出他的左右侧最靠近他的且是他的因子的位置L、R,并记录,那么对于每个数a[i]都有了他的L,R,而对于每个a[i]在f(l,r)有价值的次数之和就是(i-L+1)*(R-i+1)


代码:

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <iostream>
#include <queue>
#include <algorithm>
#include <vector>
using namespace std;
#define LL __int64
#define INF 0x3f3f3f3f
const int MAXN=100005;
#define mod 1000000007
int l[MAXN];
int r[MAXN];
int a[MAXN];
int visl[MAXN];
int visr[MAXN];
LL ans;
int main()
{
    int n,i,j;
    while(~scanf("%d",&n))
    {
        memset(visl,0,sizeof(visl));
        memset(visr,0,sizeof(visr));
        for(i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
            l[i]=1;
            r[i]=n;
        }
        for(i=1;i<=n;i++)
        {
            for(j=a[i];j<=10000;j+=a[i])
            {
                if(visr[j]&&r[visr[j]]==n)
                    r[visr[j]]=i-1;
            }
            visr[a[i]]=i;
        }
        for(i=n;i>0;i--)
        {
            for(j=a[i];j<=10000;j+=a[i])
            {
                if(visl[j]&&l[visl[j]]==1)
                    l[visl[j]]=i+1;
            }
            visl[a[i]]=i;
        }
        ans=0;
        for(i=1;i<=n;i++)
        {
            //printf("%d %d\n",l[i],r[i]);
            ans+=(LL)(i+1-l[i])*(r[i]-i+1);
            ans%=mod;
        }
        printf("%I64d\n",ans%mod);
    }
    return 0;
}


本文出自 “AcFun” 博客,请务必保留此出处http://9712315.blog.51cto.com/9702315/1677136

hdu 5288 OO’s Sequence(2015多校第一场第1题)枚举因子

标签:因子枚举

原文地址:http://9712315.blog.51cto.com/9702315/1677136

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!