hdu6287 口算训练

口算训练

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


 

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

2018CCPC女生赛题目。

题意是给你一个n个数的序列,然后q次查询,问你l到r区间里,区间内所有数相乘能否是d的倍数

 

 

思路:首先,我们必须要先知道,如何判断相乘结果是否是d的倍数,比如,我们得到区间内的这些数后,可以分解质因子,然后记录下来,同时也把d分解质因子,然后如果对于d的质因子每一个都小于等于区间内所有数分解出的质因子的个数,则证明区间所有数相乘的结果是d的倍数。

所以首先我们可以处理一遍输入的n个数,分解每个数的质因子,然后得到的因子,把这个数的下标在一个存在一个vector中,(注意,二维vector的下标v[a].push_back(i),代表的就是在第i个位置分解出a这个因子,所以把它放进去,方便后面我们进行二分查询)。接着重点来了,我们对于这个d也进行分解质因子,然后去判断这个因子的vector中,是否存在一个l到r区间中可以满足大于等于d当前的质因子的个数的,(要利用l二分在vector中的位置,r也是要二分出vector的位置,然后相减,就可以得到在l到r区间里当前质因子的个数,所以就可以去判断和当前因子的个数大小,如果全部都满足,即可以满足l到r区间,d是这些数相乘的倍数)

 

思路就是上面这样,下面具体可以看下代码 , 应该就会懂了。

 

#include <bits/stdc++.h>
using namespace std;

const int maxn = 1e5 + 10;
int n,t,m,d;
int a[maxn];
vector<int>v[maxn];
bool vis[maxn * 10];
int b[maxn];
void fj(int x,int id) {
    for (int i = 2; i * i <= x; i++)
        while (x % i == 0)
            v[i].push_back(id),  x /= i;
    if (x > 1) v[x].push_back(id);
}


bool check(int x,int l,int r,int cnt1,int d) {
    int id = x;
    vector<int>::iterator it1,it2;
    it1 = lower_bound(v[id].begin(),v[id].end(),l);
    it2 = upper_bound(v[id].begin(),v[id].end(),r);
    int num;
    num = it2 - it1;
    if(num < cnt1) {
        return false ;
    }
    return true;
}
int main() {
    scanf("%d",&t);
    while(t--) {
        for(int i = 0; i < maxn; i++) v[i].clear();
        scanf("%d%d",&n,&m);
        for(int i = 1; i <= n; i++)
            scanf("%d",&a[i]),fj(a[i],i);
        int l,r;
        while(m--) {
            int g[maxn];
            set<int>s;
            scanf("%d%d%d",&l,&r,&d);
            int x = d ;
            int flag = 0;
            for(int i = 2; i * i <= x; i++) {
                int cnt1 = 0;
                while(x % i == 0)
                    cnt1++, x /= i;

                if(!check(i,l,r,cnt1,d)) {
                    flag = 1;
                    break;
                }
            }
            if(x > 1 && !check(x,l,r,1,d)) flag = 1;
            if(!flag)printf("Yes\n");
            else printf("No\n");
        }
    } 
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值