口算训练
Time Limit: 8000/4000 MS (Java/Others) Memory Limit: 512000/512000 K (Java/Others)
Total Submission(s): 1963 Accepted Submission(s): 419
Problem Description
小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。
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
Source
与倍数有关的,应该联想到因子,不然会超时,而因子数是很少的,运用唯一分解定理求出素因子的个数,再二分查找,判断【l,r】是否能完全覆盖其素因子数
#include<bits/stdc++.h>
using namespace std;
const int maxn=1e+5;
vector<int>G[maxn];
bool check(int l,int r,int d)
{
int ans;
for(int j=2; j*j<=d; j++)
{
if(d%j==0)
{
int cnt=0;
while(d%j==0)
cnt++,d/=j;
ans=upper_bound(G[j].begin(),G[j].end(),r)-lower_bound(G[j].begin(),G[j].end(),l);
if(cnt>ans) return false;
}
}
if(d>1)
{
ans=upper_bound(G[d].begin(),G[d].end(),r)-lower_bound(G[d].begin(),G[d].end(),l);
if(ans<1) return false;
}
return true;
}
int main()
{
int t,n,m,l,r,d,x;
scanf("%d",&t);
while(t--)
{
scanf("%d%d",&n,&m);
for(int i=1;i<=maxn;i++) G[i].clear();
for(int i=1; i<=n; i++)
{
scanf("%d",&x);
for(int j=2; j*j<=x; j++)
{
while(x%j==0)
{
G[j].push_back(i);
x/=j;
}
}
if(x>1)
G[x].push_back(i);
}
while(m--)
{
scanf("%d%d%d",&l,&r,&d);
bool flag=check(l,r,d);
printf("%s\n",flag?"Yes":"No");
}
}
}