HDU 6287 口算训练

口算训练

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

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中国大学生程序设计竞赛-女生专场

题目大意:

       给出n个数，在这n个数上进行m次查询，每个查询给出三个正整数l,r,d，判断a(l)×a(l+1)×...×a(r−1)×a(r)是不是d的倍数。

解题思路：

       该题由于每一个a(i)达到了100000，如果简单模拟过程的话可能会溢出，并且查询都乘一遍会超时。故而采取分解质因数的方法。每次输入n个数，将每个数分解质因数，并且用vector数组，来存每一个质因子出现在哪个下标的数中。然后每次将d也分解质因子，将它需要的质因子的个数记录下来，然后到l，r范围中找，若个数够，找下一个，不够直接break输出No。

（一个数n的质因子大于sqrt（n）的最多只有一个）

AC代码：

#include <iostream>
#include <bits/stdc++.h>
using namespace std;
const int maxn=1e+5;
vector<int> num[maxn];//下标代表因子
int query(int ll,int rr,int xx)//该函数很巧妙，每次返回该质因子在ll-rr中出现的个数
{
    return upper_bound(num[xx].begin(),num[xx].end(),rr)-lower_bound(num[xx].begin(),num[xx].end(),ll);
}
int main()
{
    freopen("in.txt","r",stdin);
    int t;
    scanf("%d",&t);
    while(t--)
    {
        for(int i=0; i<maxn; ++i)//这里注意每次都要清空一次num数组，因为我们在这里定义的是全局数组
        {
            num[i].clear();
        }
        long long n,m;
        scanf("%lld%lld",&n,&m);
        long long temp;
        for(int i=1; i<=n; i++)
        {
            scanf("%lld",&temp);
            //将a[i]分解因式
            /*
            long long t=temp;
            for(int j=2; j*j<t; j++)此处如果这么写就会超时
            {
                while(temp%j==0)
                {
                    num[j].push_back(i);
                    temp/=j;    
                }
            }
            */
            for(int j=2; j*j<=temp; j++)//这里temp的确会发生变化,
            {                           //但是下一次进入循环时相当于是对一个新的数进行因式分解
                while(temp%j==0)        
                {
                    num[j].push_back(i);
                    temp/=j;
                }
            }
            if(temp>0)
                num[temp].push_back(i);
        }
        long long l,r,d;
        int cnt;
        int flag;
        while(m--)
        {
            scanf("%lld%lld%lld",&l,&r,&d);
            flag=1;
            for(int j=2; j*j<=d; j++)
            {
                if(d%j==0)//此处也优化了一下，因为只有在d%j==0时
                {         //才需要进行query函数的调用
                    cnt=0;
                    while(d%j==0)
                    {
                        cnt++;
                        d/=j;
                    }

//                cout<<j<<" "<<query(l,r,j)<<" "<<cnt<<endl;
                    if(query(l,r,j)<cnt)
                    {
                        flag=0;
                        break;
                    }
                }
            }
            if(d>0)
                if(query(l,r,d)>=1)
                    d/=d;
            if(d!=1)
                printf("No\n");
            else
            {
                if(flag)
                    printf("Yes\n");
                else
                    printf("No\n");
            }
        }
    }
    return 0;
}

如有说的不对的地方，还请大佬们指出。