版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/chimchim04/article/details/90067104
口算训练
小Q非常喜欢数学,但是他的口算能力非常弱。因此他找到了小T,给了小T一个长度为nn的正整数序列a1,a2,...,an要求小T抛出mm个问题以训练他的口算能力。
每个问题给出三个正整数l,r,dl,r,d,小Q需要通过口算快速判断al×al+1×...×ar−1×ar是不是d的倍数。
小Q迅速地回答了出来,但是小T并不知道正确答案是什么,请写一个程序帮助小T计算这些问题的正确答案。
Input
第一行包含一个正整数T(1≤T≤10)T(1≤T≤10),表示测试数据的组数。
每组数据第一行包含两个正整数n,m(1≤n,m≤100000)分别表示序列长度以及问题个数。
第二行包含nn个正整数a1,a2,...,an(1≤ai≤100000)表示序列中的每个数。
接下来mm行,每行三个正整数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
思路:分解每个数的质因子并记录,二分查找d的质因子在l,r之间是否够,具体实现看代码
代码:
#include<stdio.h>
#include<string.h>
#include<vector>
#include<algorithm>
using namespace std;
vector<int>e[100000];
void init(int x,int d)
{
for(int i=2;i*i<=x;i++)
{
if(x%i==0)
{
while(x%i==0)
{
x/=i;
e[i].push_back(d);
}
}
}
if(x>1) e[x].push_back(d);
}
int fun(int l,int r,int d)
{
for(int i=2;i*i<=d;i++)
{
if(d%i==0)
{
int k=0;
while(d%i==0)
{
d/=i;
k++;
}
int cn=upper_bound(e[i].begin(),e[i].end(),r)-lower_bound(e[i].begin(),e[i].end(),l); //l-r中的质因子i的个数
if(cn<k) return 0; //如果不够,那就不是d的倍数
}
}
if(d>1)
{
int cn=cn=upper_bound(e[d].begin(),e[d].end(),r)-lower_bound(e[d].begin(),e[d].end(),l);
if(cn<1) return 0;
}
return 1;
}
int main()
{
int t;
scanf("%d",&t);
while(t--)
{
for(int i=1;i<=100000;i++) //初始化
e[i].clear();
int n,m,x;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%d",&x);
init(x,i); //记录质因子
}
int l,r,d;
for(int i=0;i<m;i++)
{
scanf("%d%d%d",&l,&r,&d);
int f=fun(l,r,d);
if(f) printf("Yes\n");
else printf("No\n");
}
}
return 0;
}