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计算这些问题的正确答案。

 

Input

第一行包含一个正整数T(1≤T≤10)T(1≤T≤10)，表示测试数据的组数。 

每组数据第一行包含两个正整数n,m(1≤n,m≤100000)n,m(1≤n,m≤100000)，分别表示序列长度以及问题个数。 

第二行包含nn个正整数a1,a2,...,an(1≤ai≤100000)a1,a2,...,an(1≤ai≤100000)，表示序列中的每个数。 

接下来mm行，每行三个正整数l,r,d(1≤l≤r≤n,1≤d≤100000)l,r,d(1≤l≤r≤n,1≤d≤100000)，表示每个问题。

 

Output

对于每个问题输出一行，若是倍数，输出Yes，否则输出No。

 

Sample Input

1

5 4

6 4 7 2 5

1 2 24

1 3 18

2 5 17

3 5 35

 

Sample Output

Yes

No

No

Yes

这里直接gcd暴力肯定是会超时的。题目数据量最大是100000，他的质因数的个数肯定小于等于70。因为1-100000总共也才70个素数；）。

我们把输入的所有数分解后放入vector<int > g[xi]中       g[xi]中存下质因数xi的下标 。

然后我们对输入的值x也进行分解，并计数为cot。 在l-r的范围中 如果质因数的个数小于cot。那么这个数肯定不可能整除，一直将x分解完即可。

#include <queue>
#include <cstdio>
#include <set>
#include <string>
#include <stack>
#include <cmath>
#include <climits>
#include <map>
#include <cstdlib>
#include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
#include <stdio.h>
#include <ctype.h>
#include <bitset>
#include <time.h>
#define LL long long
#define ULL unsigned long long
#define mod 1000000007
#define INF 0x3ffffff
#define mem(a,b) memset(a,b,sizeof(a))
#define MODD(a,b) (((a%b)+b)%b)
#define PI 3.1415927
#define maxn 100005
using namespace std;
char a[maxn];
int n;
int r[maxn];
vector<int> g[maxn];
void di_prime(int x,int n)
{
    int p = 0;
    for(int i = 2; i * i <= n; i++){
        if(n % i == 0){
           while(n % i == 0){
            g[i].push_back(x);
             n/=i;
           }
        }

    }
    if(n > 1) g[n].push_back(x);

}
int findP(int x,int l,int r)
{
    return (upper_bound(g[x].begin(),g[x].end(),r) - lower_bound(g[x].begin(),g[x].end(),l));//在l-r中 质因数x的个数
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--){
        int n,m;
        scanf("%d%d",&n,&m);
        for(int i = 0; i <= maxn; i++) g[i].clear();
        for(int i = 1; i <= n; i++){
            int k;
            scanf("%d",&k);
            di_prime(i,k);
        }
        for(int i = 0; i < m; i++){
            int l,r,x;
            int flag = 0;
            scanf("%d%d%d",&l,&r,&x);
            for(int j = 2; j * j <= x; j++){
                int cot = 0;
                if(x % j == 0){
                    while(x % j == 0){
                      cot++;
                      x/=j;
                    }
                    if(cot > findP(j,l,r)){
                        flag = 1;
                        break;
                    }

                }

            }
            if(x > 1 && findP(x,l,r) == 0) flag = 1;
            if(!flag) printf("Yes\n");
            else printf("No\n");

        }



    }

    return 0;

}