CF 112D. Petya and Divisors 枚举(标记倍数最后出现位置)

题意:n次询问[x[i],y[i]],每次问x[i]有多少个因子不能整除 x[i-y[i]]...x[i-1] 中的任意一个?
n,x[i]<=1e5,y[i]<=i-1,5秒的时限.


暴力的做法:算出x[i]的因子 然后遍历[ x[i-y[i],..x[i-1]]判定 最坏(D*n^2).
因为因子的个数不会很多 && x[i]<=1e5 我们遍历因子x的倍数kx 需要知道那些倍数在范围内?

如果last[kx]最后一次位置>=i-y[i] 则淘汰该因子.查询完一个数以后更新last即可.O(n*Dlogn)

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+2,M=1e5;
int n,x[N],y[N],fac[N],vis[N],last[N];
int main()
{
    int n,mx=0;
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d%d",&x[i],&y[i]);
    for(int k=1;k<=n;k++)
    {
        int a=x[k],pn=0;
        for(int i=1;i*i<=a;i++)
        {
            if(a%i==0)
            {
                fac[++pn]=i;
                if(i*i!=a)
                    fac[++pn]=a/i;
            }
        }

        for(int i=1;i<=pn;i++)
            vis[i]=0;

        for(int i=1;i<=pn;i++)
        {
            if(i==1)
            {
                if(y[k])
                    vis[i]=1;
                continue;
            }
            int b=fac[i];
            for(int j=b;j<=mx;j+=b)
            {
                if(last[j]>=k-y[k])
                {
                    vis[i]=1;
                    break;
                }
            }
        }
        int ans=0;
        for(int i=1;i<=pn;i++)
        {
            if(!vis[i])
                ans++;
            last[fac[i]]=k;
        }
        printf("%d\n",ans);
        last[a]=k,mx=max(mx,a);
    }
    return 0;
}

use[j] 最后一个j的倍数的位置

#include <cstdio>
const int maxn = 100000 + 5;
int used[maxn];
int main(void) {
	int n;
	scanf("%d", &n);
	for (int i = 1; i <= n; i++) {
		int x, y;
		scanf("%d%d", &x, &y);
		int res = 0;
		for (int j = 1; j * j <= x; j++) {
			if (x % j == 0) {
				if (used[j] < i - y) {
					res++;
				}
				used[j] = i;
				if (j * j == x) {
					continue;
				}
				if (used[x / j] < i - y) {
					res++;
				}
				used[x / j] = i;
 			}
		}
		printf("%d\n", res);
	}
	return 0;

}