bzoj5358

Problem A. 口算训练
Input file: stdin
Output file: stdout
Time limit: 5 seconds
Memory limit: 512 megabytes
小 Q 非常喜欢数学，但是他的口算能力非常弱。因此他找到了小 T，给了小 T 一个长度为 n 的正
整数序列 a1, a2, ..., an，要求小 T 抛出 m 个问题以训练他的口算能力。
每个问题给出三个正整数 l, r, d，小 Q 需要通过口算快速判断 al × al+1 × ... × ar−1 × ar 是不是 d
的倍数。
小 Q 迅速地回答了出来，但是小 T 并不知道正确答案是什么，请写一个程序帮助小 T 计算这些问
题的正确答案。
Input
第一行包含一个正整数 T(1 ≤ T ≤ 10)，表示测试数据的组数。
每组数据第一行包含两个正整数 n, m(1 ≤ n, m ≤ 100000)，分别表示序列长度以及问题个数。
第二行包含 n 个正整数 a1, a2, ..., an(1 ≤ ai ≤ 100000)，表示序列中的每个数。
接下来 m 行，每行三个正整数 l, r, d(1 ≤ l ≤ r ≤ n, 1 ≤ d ≤ 100000)，表示每个问题。
Output
对于每个问题输出一行，若是倍数，输出“Yes”，否则输出“No”。
Examples
stdin

1
5 4
6 4 7 2 5
1 2 24
1 3 18
2 5 17
3 5 35

stdout
Yes
No
No
Yes

 

sol：这题数据 ai 很小，所以XJB搞搞就能过，就把<=sqrtai的ai的质因数筛出来，然后记个数组，对于>=sqrtai的约数肯定只有一个，所以开个vector特判一下就能过了

#include <bits/stdc++.h>
using namespace std;
typedef int ll;
inline ll read()
{
    ll s=0;
    bool f=0;
    char ch=' ';
    while(!isdigit(ch))
    {
        f|=(ch=='-'); ch=getchar();
    }
    while(isdigit(ch))
    {
        s=(s<<3)+(s<<1)+(ch^48); ch=getchar();
    }
    return (f)?(-s):(s);
}
#define R(x) x=read()
inline void write(ll x)
{
    if(x<0)
    {
        putchar('-'); x=-x;
    }
    if(x<10)
    {
        putchar(x+'0');    return;
    }
    write(x/10);
    putchar((x%10)+'0');
    return;
}
#define W(x) write(x),putchar(' ')
#define Wl(x) write(x),putchar('\n')
const int N=200005,B=505;
int n,m;
int Prim[N];
bool Bo[N];
int S[N][B];
vector<int>Pos[N];
inline void Shai(int Up)
{
    int i,j;
    Bo[0]=Bo[1]=1;
    for(i=2;i<=Up;i++)
    {
        if(!Bo[i]) Prim[++*Prim]=i;
        for(j=1;j<=*Prim&&(Prim[j]*i<=Up);j++)
        {
            Bo[Prim[j]*i]=1;
            if(i%Prim[j]==0) break;
        }
    }
    return;
}
int main()
{
    int i,j;
    Shai(500);
    R(n); R(m);
    for(i=1;i<=n;i++)
    {
        int x=read();
        for(j=1;j<=*Prim;j++)
        {
            S[i][j]=S[i-1][j];
            if(x%Prim[j]==0)
            {
                while(x%Prim[j]==0)
                {
                    S[i][j]++; x/=Prim[j];
                }
            }
        }
        if(x>1) Pos[x].push_back(i);
    }
    while(m--)
    {
        int l=read(),r=read(),x=read();
        bool Bo=0;
        for(j=1;j<=*Prim;j++) if(x%Prim[j]==0)
        {
            int oo=0;
            while(x%Prim[j]==0)
            {
                oo++; x/=Prim[j];
            }
            if(S[r][j]-S[l-1][j]<oo)
            {
                Bo=1; break;
            }
        }
        if(x>1)
        {
            vector<int>::iterator it=lower_bound(Pos[x].begin(),Pos[x].end(),l);
            int Wz; 
            if(it!=Pos[x].end()) Wz=(*it);
            else Wz=n+1;
            if(Wz>r) Bo=1;
        }
        (Bo)?(puts("No")):puts("Yes");
    }
    return 0;
}
/*
input
5 4
6 4 7 2 5
1 2 24
1 3 18
2 5 17
3 5 35
output
Yes
No
No
Yes

input
5 1
2087 2 2 2 2
1 2 4174
output
Yes
*/