hdu 6287 口算训练(二分+质因数分解+思维)

口算训练

Time Limit: 8000/4000 MS (Java/Others)    Memory Limit: 512000/512000 K (Java/Others)
Total Submission(s): 1968    Accepted Submission(s): 423

 

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

"字节跳动杯"2018中国大学生程序设计竞赛-女生专场

我也是看了题解才会做的，说到底还是数学太差了...

题意：很明显了...不多解释

思路：有很多样例，每组有m个数，n个要求的式子，r为左边界,l为右边界,r乘到l的值为需要求的数是否为其p倍数（好像讲的不是特别好...大概就是这样啦）

1.我们可以先把m个数分别对其求质因数，存入G[i]中,i为质因数，G[i]为下标。

2.再对p求质因数。

3.关键！！然后对询问的数分解质因式，对每一个分解得到的质因数，如果这个质因数的个数小于等于序列中该质因数的个数就可以整除，否则不能整除。这里用二分查找法查找在序列中质因数的个数。

好像讲的还是不太清楚...还是上代码吧！！

这个讲的挺好的,详见：https://blog.csdn.net/wjmwsgj/article/details/80519183

#include<iostream>
#include<algorithm>
#include<vector>
#include<stdio.h>
using namespace std;
vector<int>G[100001];
void fenjie(int id,int t)//求质因数
{
	for(int i=2;i*i<=t;i++)
	{
		while(t%i==0)
		{
			G[i].push_back(id);//加入到其尾部,G[x]含有质因数x的下标
			t=t/i;
		}
	}
	if(t>1)//注意！！
	{
		G[t].push_back(id);
	}
}
int main()
{
	int m,n,t,flag;
	int count;
	scanf("%d",&count);
	while(count--)
	{
		scanf("%d%d",&m,&n);
		for(int i=0;i<100001;i++)清空容器
		{
			G[i].clear();
		}
		for(int i=0;i<m;i++)//求质因数
		{
			scanf("%d",&t);
			fenjie(i,t);
		}
		while(n--)
		{
			int l,r,p;
			flag=1;
			scanf("%d%d%d",&l,&r,&p);
			l--;
			r--;
			for(int i=2;i*i<=p;i++)//对其求解质因数
			{
				int ans=0;
				while(p%i==0)
				{
					ans++;
					p=p/i;
				}
				if(ans!=0)
				{
					int sum=upper_bound(G[i].begin(),G[i].end(),r)-lower_bound(G[i].begin(),G[i].end(),l);//二分查找
					if(sum<ans)
					{
						flag=0;
						break;
					}
				}
			}
			if(p>1)//注意！！
			{
				int sum=upper_bound(G[p].begin(),G[p].end(),r)-lower_bound(G[p].begin(),G[p].end(),l);//二分查找
				if(sum==0)
				{
					flag=0;
				}
			}
			if(flag==0)
			{
				printf("No\n");
			}
			else
			{
				printf("Yes\n");
			}
		}
	}
	return 0;
}