HDU 5288 乱搞

最新推荐文章于 2018-07-16 20:24:05 发布

Dan__ge

最新推荐文章于 2018-07-16 20:24:05 发布

阅读量632

点赞数

分类专栏：杂文章标签： hdu ACM

本文链接：https://blog.csdn.net/Dan__ge/article/details/51801715

版权

杂专栏收录该内容

24 篇文章 0 订阅

订阅专栏

点击打开链接

题意：问你给定的方程式最后的值，表达式的含义是在区间L到R中有多少个元素不是其它任意一个元素的倍数，如2 3 4 5 ，结果就是3，因为2，3，5不是剩下的元素的倍数

思路：区间太多了，所以不用考虑区间的做法，我们可以考虑每个元素它可以成功的次数，而它可以成功的条件是这个区间没有它的素因子，那么我们可以找到它左侧第一个是它素因子的位置以及右侧第一个是它素因子的位置，那么它可以提供的次数就是（R-i）*（i-L），最后加起来就行了，而处理的时候复杂度为N*素因子个数，因为数字最大10000；所以素因子不会很大满足时限

#include <math.h>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3fll;
const int maxn=100010;
const ll mod=1000000007;
int num[maxn],vis[maxn],L[maxn],R[maxn];
int main(){
    int n;
    while(scanf("%d",&n)!=-1){
        for(int i=0;i<n;i++) scanf("%d",&num[i]);
        memset(vis,-1,sizeof(vis));
        for(int i=0;i<n;i++){
            int t=(int)sqrt((double)num[i]),l=-1;
            for(int j=1;j<=t;j++){
                if(num[i]%j==0){
                    if(vis[j]!=-1) l=max(l,vis[j]);
                    if(vis[num[i]/j]!=-1) l=max(l,vis[num[i]/j]);
                }
            }
            vis[num[i]]=i;
            if(l==-1) L[i]=-1;
            else L[i]=l;
        }
        memset(vis,-1,sizeof(vis));
        for(int i=n-1;i>=0;i--){
            int t=(int)sqrt((double)num[i]),r=inf;
            for(int j=1;j<=t;j++){
                if(num[i]%j==0){
                    if(vis[j]!=-1) r=min(r,vis[j]);
                    if(vis[num[i]/j]!=-1) r=min(r,vis[num[i]/j]);
                }
            }
            vis[num[i]]=i;
            if(r==inf) R[i]=n;
            else R[i]=r;
        }
        ll ans=0;
        for(int i=0;i<n;i++) ans=ans+(ll)(i-L[i])*(ll)(R[i]-i)%mod;
        printf("%I64d\n",ans%mod);
    }
    return 0;
}