【多校第一场】【筛法改】HDU 5288 OO’s Sequence

这道题的数据范围需要让我们有些反应才是，虽然n的组数是1e5的大小，但是数字的范围是1e4。

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

思路是找到每一个元素的左端首个因数的位置和右端因数的位置，这是显而易见的。但是这题的问题却不是在这里，而是在于如何操作。学长介绍的方法是将他们从小到大排列，然后去寻找当前值的倍数。。。有点说不太清不过对着代码模拟一遍一个就好，注释才是正文系列。

记住，long long不要乱开，sort时记得你的数组是从哪开始的。（因为这些而WA了一天TAT）。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#define LL long long
using namespace std;

struct n
{
    int num;
    int id;
}num[1000100];
LL ans;
const int mod=1e9+7;
vector <int>lr[10002];

bool cmp(n a,n b)
{
    return a.num<b.num;
}

int main()
{
    int n;
	while(~scanf("%d",&n))
	{
	    for(int i=0;i<=10001;i++)
            lr[i].clear();
        ans=0;
	    for(int i=1;i<=n;i++)
        {
            scanf("%d",&num[i].num);
            num[i].id=i;
        }

        sort(num+1,num+n+1,cmp);

        int k=1;
        for(int i=1;i<=n;i++)
        {
            for(int j=num[i].num;j<=10000;j+=num[i].num)
            {
                lr[j].push_back(num[i].id);
            }
            if(i!=n && num[i].num==num[i+1].num) continue;

            while(k<=i)
            {
                int l=0,r=n+1;
                for(int j=0;j<lr[num[k].num].size();j++)
                {
                    if(lr[num[k].num][j]<num[k].id && lr[num[k].num][j]>l)
                    {
                        l=lr[num[k].num][j];
                    }
                    else if(lr[num[k].num][j]>num[k].id && lr[num[k].num][j]<r)
                    {
                        r=lr[num[k].num][j];
                    }
                }
                ans=(ans+(num[k].id-l)*(r-num[k].id))%mod;
                k++;
//                cout<<l<<" "<<r<<" "<<ans<<endl;
//                for(int i=1;i<=n;i++)
//                {
//                    printf("%d  |",lr[i].size());
//                    for(int j=0;j<lr[i].size();j++) printf(" %d",lr[i][j]);
//                    printf("\n");
//                }
            }
        }

//        for(int i=1;i<=n;i++)
//        {
//            printf("%d  |",lr[i].size());
//            for(int j=0;j<lr[i].size();j++) printf(" %d",lr[i][j]);
//            printf("\n");
//        }

        printf("%lld\n",ans);
	}
	return 0;
}