2021年中国大学生程序设计竞赛 女生专场 - 热身赛 Problem C. 口算训练(质因子分解)
Posted Frozen_Guardian
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了2021年中国大学生程序设计竞赛 女生专场 - 热身赛 Problem C. 口算训练(质因子分解)相关的知识,希望对你有一定的参考价值。
题目分析:判断
x
x
x 是
y
y
y 的倍数,等价于质因子分解后,
y
y
y 中的每个质因子的出现次数都小于等于其在
x
x
x 中的出现次数。
那么对于每次询问 [ l , r , d ] [l,r,d] [l,r,d],实质上就是将 d d d 质因子分解后,判断一下每个质因子在区间 [ l , r ] [l,r] [l,r] 中的出现次数。
问题转换为了如何快速求解区间内某个质因子的出现次数。
因为本题的值域特别小,所以考虑按值域分块, ⌊ 100000 ⌋ = 316 \\lfloor \\sqrt100000 \\rfloor=316 ⌊100000⌋=316:
- 小于等于 316 316 316 的质因子只有 65 65 65 个,预处理出 65 65 65 个前缀和暴力维护即可。
- 大于 316 316 316 的质因子在 d d d 中至多出现一次,所以问题又转换为了判断区间存在问题。
对于上面的第一种情况,暴力 O ( n ∗ 65 ) O(n*65) O(n∗65) 暴力维护即可;第二种情况的话也可以用二分轻松解决,时间复杂度 O ( n log n ) O(n\\log n) O(nlogn),空间复杂度均摊下来是 O ( n ) O(n) O(n) 的。
然后 d d d 的质因子至多有 log d \\log d logd 个,所以对于每次询问的复杂度就是 O ( log d + log n ) O(\\log d+\\log n) O(logd+logn)
其实可以将本题带修,前缀和改成线段树,vector 改成 set,复杂度多一个 log n \\log n logn 而已
代码:
// #pragma GCC optimize(2)
// #pragma GCC optimize("Ofast","inline","-ffast-math")
// #pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include<iostream>
#include<cstdio>
#include<string>
#include<ctime>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<stack>
#include<climits>
#include<queue>
#include<map>
#include<set>
#include<sstream>
#include<cassert>
#include<bitset>
#include<list>
#include<unordered_map>
#define lowbit(x) (x&-x)
using namespace std;
typedef long long LL;
typedef unsigned long long ull;
template<typename T>
inline void read(T &x)
T f=1;x=0;
char ch=getchar();
while(0==isdigit(ch))if(ch=='-')f=-1;ch=getchar();
while(0!=isdigit(ch)) x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
x*=f;
template<typename T>
inline void write(T x)
if(x<0)x=~(x-1);putchar('-');
if(x>9)write(x/10);
putchar(x%10+'0');
const int inf=0x3f3f3f3f;
const int N=1e5+100;
const int pcnt=65;
const int pri[]=2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,97,101,103,107,109,113,127,131,137,139,149,151,157,163,167,173,179,181,191,193,197,199,211,223,227,229,233,239,241,251,257,263,269,271,277,281,283,293,307,311,313;
int sum[65][N];
vector<int>node[N];
int main()
#ifndef ONLINE_JUDGE
// freopen("data.in.txt","r",stdin);
// freopen("data.out.txt","w",stdout);
#endif
// ios::sync_with_stdio(false);
int w;
cin>>w;
while(w--)
memset(sum,0,sizeof(sum));
for(int i=0;i<N;i++)
node[i].clear();
int n,m;
read(n),read(m);
for(int i=1;i<=n;i++)
int x;
read(x);
for(int j=0;j<pcnt;j++)
while(x%pri[j]==0)
x/=pri[j];
sum[j][i]++;
if(x>1)
node[x].push_back(i);
for(int j=0;j<pcnt;j++)
for(int i=1;i<=n;i++)
sum[j][i]+=sum[j][i-1];
while(m--)
int l,r,x;
bool flag=true;
read(l),read(r),read(x);
for(int i=0;i<pcnt;i++)
if(x%pri[i]==0)
int cnt=0;
while(x%pri[i]==0)
x/=pri[i];
cnt++;
flag&=(sum[i][r]-sum[i][l-1])>=cnt;
if(x>1)
int pos=*lower_bound(node[x].begin(),node[x].end(),l);
if(pos>r)
flag=false;
puts(flag?"Yes":"No");
return 0;
以上是关于2021年中国大学生程序设计竞赛 女生专场 - 热身赛 Problem C. 口算训练(质因子分解)的主要内容,如果未能解决你的问题,请参考以下文章
2021年中国大学生程序设计竞赛 女生专场 - 热身赛 Problem C. 口算训练(质因子分解)
"巴卡斯杯" 中国大学生程序设计竞赛 - 女生专场(重现)解题思路