hdu 6287 口算训练

题意：

小Q非常喜欢数学，但是他的口算能力非常弱。因此他找到了小T，给了小T一个长度为nn的正整数序列a1,a2,...,ana1,a2,...,an，要求小T抛出mm个问题以训练他的口算能力。 

每个问题给出三个正整数l,r,dl,r,d，小Q需要通过口算快速判断al×al+1×...×ar−1×aral×al+1×...×ar−1×ar是不是dd的倍数。 

小Q迅速地回答了出来，但是小T并不知道正确答案是什么，请写一个程序帮助小T计算这些问题的正确答案。

思路：

一开始bit套map 用前缀和处理，tle了？

对于每一个输入的数字，分解质因子，对于每一个质因子，把这个数字的位置push到这个质因子的vector里面去。

那么每一个vector里面的数字就是非递减的了。

然后对于每一个输入的l，r，d，把d分解质因子，对于每一个质因子，在这个质因子的vector里面找第一个大于等于l的位置和最后一个小于等于r的位置，相减，比较是否有d中这个质因子的数量多。

找的时候当然利用前面的非递减，所以可以二分。

代码：

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <vector>
using namespace std;
const int N = 1e5 + 10;

vector<int> p;
vector<int> G[N],g[N];
bool pp[N];
int a[N];
int n,m;

int main()
{
    for (int i = 2;i < N;i++)
    {
        if (!pp[i])
        {
            p.push_back(i);
            for (int j = i * 2;j < N;j += i)
            {
                G[j].push_back(i);
                pp[j] = 1;
            }
        }
    }
    int T;
    scanf("%d",&T);
    while (T--)
    {
        scanf("%d%d",&n,&m);
        for (int i = 0;i < p.size();i++)
        {
            g[p[i]].clear();
        }
        for (int i = 1;i <= n;i++)
        {
            int x;
            scanf("%d",&x);
            if (x == 1) continue;
            if (!pp[x])
            {
                g[x].push_back(i);
                continue;
            }
            for (int j = 0;j < G[x].size();j++)
            {
                int cnt = 0;
                int tmp = x;
                int v = G[x][j];
                while (tmp % v == 0)
                {
                    tmp /= v;
                    g[v].push_back(i);
                }
            }
        }
        while (m--)
        {
            int l,r,d;
            scanf("%d%d%d",&l,&r,&d);
            if (d == 1)
            {
                puts("Yes");
                continue;
            }
            if (!pp[d])
            {
                int x = lower_bound(g[d].begin(),g[d].end(),l) - g[d].begin();
                int y = upper_bound(g[d].begin(),g[d].end(),r) - g[d].begin();
                y--;
                if (y - x + 1 < 1) puts("No");
                else puts("Yes");
            }
            else
            {
                bool f = 0;
                for (int i = 0;i < G[d].size();i++)
                {
                    int v = G[d][i];
                    int cnt = 0;
                    int tmp = d;
                    while (tmp % v == 0)
                    {
                        cnt++;
                        tmp /= v;
                    }
                    int x = lower_bound(g[v].begin(),g[v].end(),l) - g[v].begin();
                    int y = upper_bound(g[v].begin(),g[v].end(),r) - g[v].begin();
                    y--;
                    if (y - x + 1 < cnt) f = 1;
                    if (f) break;
                }
                if (f) puts("No");
                else puts("Yes");
            }
        }
    }
    return 0;
}