2023牛寒--阿宁的倍数

阿宁有一个长度为n的数组a，下标从1开始，有q次操作。
修改操作：数组末尾增加一个数x，数组长度加1。
询问操作：有多少个i(i>x)，满足ai是ax​的倍数？

输入描述：
第一行两个整数n,q。
第二行n个整数ai。
接下来q行，每行两个数op,x，代表一次操作。
如果op是1代表是修改操作；如果是2代表是询问操作。
1≤n,q,ai≤2×10^5
1≤op≤2
保证至少有一次询问操作。
如果是修改操作，1≤x≤2×10^5。
如果是询问操作，x小于等于当前数组长度。

输出描述：

对于每次询问操作，输出一个整数，占一行。

输入：

5 5
2 4 1 5 6
2 1
1 5
1 1
2 3
2 4

输出：

2
4
1

解析：对于每个ai，设他的因子是j，那么我们存v[ j ].push( i )，想当初村下了 j 的倍数的下标，然后每次查询就是在在他倍数所有下标中大于x的个数，此时我们可以二分找到第一个大于等于x+1的位置，用v[ j ]的size减去该位置就是个数。

注意：a[ ]数组要开到4e5，因为操作可以使得a[ ]长度增加，防止段错误。

/*
   Author:  XiaoWhy
   Created: 2023.2.3
*/
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=4e5+5,M=2e5+5;
int a[N];
vector<int> v[N];
void solve()
{
	int n,q;
	scanf("%d%d",&n,&q);
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&a[i]);
		for(int j=1;j<=a[i]/j;j++)
		{
			if(a[i]%j==0)//因子
			{
				v[j].push_back(i);
				if(a[i]/j!=j) v[a[i]/j].push_back(i);
			}
		}
	}
	int m=n;//当前数组长度
	while(q--)
	{
		int op,x;
		scanf("%d%d",&op,&x);
		if(op==1)
		{
			a[++m]=x;
			for(int j=1;j<=x/j;j++)
			{
				if(x%j==0){
					v[j].push_back(m);
					if(x/j!=j) v[x/j].push_back(m);
				}
			}
		}else
		{
			int w=a[x];//所查询a[x]
			int l=0,r=v[w].size()-1,pos=0,k=x+1;
			if(r==-1||k>v[w][r]||k<v[w][0])//不存在倍数情况或者k不存在
			{
				printf("0\n");
				continue;
			}
			while(l<=r)
			{
				int mid=l+r>>1;
				if(v[w][mid]>=k) pos=mid,r=mid-1;
				else l=mid+1;
			}
			printf("%d\n",v[w].size()-pos);
		}
	}
}
int main()
{
	solve();
	return 0;
}