Strange Definition CodeForces - 1471D

题意:

定义数字 x 和 y 是“相邻”的当且仅当 lcm(x,y)/gcd(x,y) 是一个平方数。
给定一个长度为 n 的数组 a。
每过一秒,数组 a 会发生变化:ai 会变成数组 a 中与其“相邻”的所有数字的乘积。
定义 di 为数组 a 中与 ai “相邻” 的数字个数。
定义数组 a 的美丽值为 max1≤i≤ndi
给出 q 个询问,每次询问给出当前时间 w,问当前数组 a 的美丽值。

题解:

参考题解
不会做。。emmm。。。看完题解恍然大悟
题目中w的范围0≤w≤1018,那说明不可能真的一轮一轮进行,因为w太大了,可能会存在a数组发生优先次变化后,答案就固定了
题目所定义的“相邻”经过简化后可得
在这里插入图片描述
分母为(gcd(x,y))2,分子为x * y,为了让整体变成平方数,那说明x * y就是平方数,根据唯一分解定理,任何一个数都可以分解成质因数相乘的形式,x * y 为平方数,说明x和y对应的质因子的幂次奇偶性是相同的。(就比如x可以分解出一个3,那么y肯定也要分解出一个3,这样才能凑出平方)
现在我们对x和y进行处理,将x和y中偶数次幂的质因子删除掉,只留下奇数次幂的质因子。(因为偶数次的本身就是平方数了,对题目没有影响,有影响的是奇数)
处理后的到X和Y,如果X = = Y,x * y 就是平方数,否则就不是
接下来我们来处理美丽值这个问题
在a数组中没出现次数最多的数字的出现次数就是数组a的美丽值
a数组每轮都是会发生改变的,我们前面说了只有X== Y时,X与Y才是相邻的,那么与a[i]相邻的有d[i]个数,a[i]就变成了a[i]d[i]
我们接着对处理后的A数组进行质因数删除,只保留奇数次幂的质因子
如果d[i]是偶数,那么a[i]就变成了1,否则还是a[i]
我们发现变成1的数,之后再也不会发生改变了,相当于固定了,而没变成1的数永远都不是1,(是它本身)
因为没变成1的数说明是奇数个,那么与A[i]相邻的数也是奇数个,相当于A[i]变化后是A[i]d[i],d[i]为奇数,而A[i]d[i]经过质因数删除操作后又得到A[i],所以永远不可能是1
说明:变化1次和变化多次的答案是一样的
也就是当w>=1时,我们只需要比较不变化的1的个数和变化1次后1的个数,输出最大值就行
变化后1的数量 = 还未变化时出现次数为偶数的非1的数字 + 原本1的数量
题目非常巧妙,好好揣摩

代码:

情况数组时最好不要用memset了,特别是大数据。。。。
不然
在这里插入图片描述

#include<bits/stdc++.h>
typedef long long ll;
using namespace std;
inline int read(){int s=0,w=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();//s=(s<<3)+(s<<1)+(ch^48);return s*w;
}
const int maxn=1e6+9;
int prime[maxn],tot,tag[maxn];
void Prime(int N){tag[0]=tag[1]=1;for(int i=2;i<=N;i++){if(tag[i]==0)prime[tot++]=i;for(int j=0;j<tot&&i*prime[j]<N;j++){tag[i*prime[j]]=1;if(i%prime[j]==0)break;}}
}
int a[maxn],sum[maxn];
vector<int>vec;
int main()
{int T;T=read();Prime(1000000);while(T--){vec.clear();int n;n=read();int maxx=0;//maxx表示出现此处最多的数字的出现次数 int one=0;//记录没变化时1的数量 //memset(a,0,sizeof(a));for(int i=1;i<=n;i++)a[i]=read();for(int i=1;i<=n;i++){int tmp=1;//tmp为奇数个质因子的乘积 int j=0;while(a[i]>=prime[j]){int num=0;while(a[i]%prime[j]==0){num++;a[i]/=prime[j];}if(num%2==1)tmp=tmp*prime[j];//如果该质因数分解除奇数j++;if (!tag[a[i]]) {tmp *= a[i];break;}}//sum用来处理后,有多少数是相同的 vec.push_back(tmp); maxx=max(maxx,++sum[tmp]);}int even=0;for(int x:vec){if(x!=1){if(sum[x]%2==0)///未变化时出现次数为偶数且非 1 的数字数量even++; }else{one++;}}for (int v : vec) {sum[v] = 0;}int q;q=read();while(q--){ll w;scanf("%lld",&w);if(w==0)printf("%d\n",maxx);else printf("%d\n",max(maxx,one+even));}}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318329.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

学习三分 (概念 + 模板 + 例题:曲线)

这好像是我第一次尝试写一个新知识入门 而不是习题解 文章目录三分概念模板例题&#xff1a;曲线题目题解代码实现三分概念 我们都知道&#xff0c;二分是在一个单调函数&#xff08;即一次函数&#xff09;上通过每次查找折半的方式&#xff0c;对答案进行搜索查找。那么&am…

迎元旦,庆surging 1.0发布

一位摄影程序员的独白每个人都有爱好&#xff0c;都有释放压力的活动&#xff0c;而我也不例外&#xff0c;我除了每天上班&#xff0c;周末就会约一群好友去拍妹子&#xff0c;成家后&#xff0c;就改为拍虫子&#xff0c;一拍就到了30岁&#xff0c;到了30岁就感觉到了中年的…

YBTOJ:求好元素(哈希表)

文章目录题目描述解析代码题目描述 解析 如果枚举m,n,p的话是n3的 会超时 但我们注意到右边查询只有O(n) 这就很不平衡 所以考虑把p移到右边&#xff0c;预处理枚举m、n存到哈希表中 查询枚举i、p 这样就把复杂度均摊降到了n2 但是本题数据较强 所以我们得保证哈希表的单词查询…

Strange Shuffle CodeForces - 1471E(交互题)

交互题 这类型不同于普通的题。 可以理解为有个问题需要你解决&#xff0c;你通过输入某些东西表示你要问系统的问题&#xff0c;这时系统会回答你的问题。在代码中的回答方式就是会输入某个东西就是系统给你的答案&#xff0c;通过这些信息你可以得到问题的解你是不可以自己测…

学习KMP (概念 + 模板 + 例题: 子串查找)

我又回来了&#xff0c;感jio这几天有点勤啊&#xff01;&#xff01; 这一次我带着KMP来了&#xff0c; 文章目录KMP介绍模板例题&#xff1a; 子串查找题目暴力题解KMP题解代码实现KMP介绍 KMP&#xff0c;即 Knuth-Morris-Pratt 字符串查找算法&#xff0c;由Donald Knuth…

直播预告 - 博时基金DevOps体系建设和自动化测试分享

最近几年&#xff0c;基金行业发展比较快&#xff0c;业务范围从传统公募到大资管&#xff0c;业务地域从中国大陆到全球化&#xff0c;在互联网金融浪潮中扮演了重要角色&#xff0c;金融科技又带来新的挑战和机遇。据毕马威2014年研究报告显示&#xff0c;由于新技术、人口变…

YBTOJ:公共子串(KMP)

文章目录题目描述解析代码题目描述 解析 我的做法是二分答案&#xff0c;暴力枚举子串&#xff0c;对其他的字符串进行KMP。 复杂度logL*L2 没有一次AC&#xff0c;因为忘记了字符串清零&#xff01;&#xff0c;否则它会和上次二分的字符串没被覆盖到的部分连在一起。 之前遇…

Asp.Net Core + Docker 搭建

一、Docker安装&#xff08;1&#xff09;机器准备 本次实验&#xff0c;我是在Win10上,使用Vmware Workstation虚拟机安装的Ubuntu 18.04.l LTS版本进行搭建。&#xff08;2&#xff09;安装Docker# yum install docker注意:如果提示未安装yum 使用命令:apt-get install yum&a…

数数字

今天中午1小时&#xff0c;定时两道&#xff0c;我就知道自己的太阳降落了。。。orz 文章目录题目题解代码实现题目 PB 带来了若干只蒟蒻。 众所周知&#xff0c;NTF 是数论学会的会长&#xff0c;于是 PB 准备用数字击败 NTF&#xff0c;以证明 PB 比 NTF 更强。 于是 PB 准…

YBTOJ:前缀数组(KMP)

文章目录题目描述解析代码题目描述 解析 题面脸上写着5个大字&#xff1a;我是KMP 但是本题没有自己做出来。。。我一开始的思路其实很接近题解了&#xff0c;只是被我舍弃了qwq。 后来卡在暴力nL2的瓶颈上&#xff0c;用了个倍增的诡异操作搞到了nLlogL&#xff0c;但是n5&am…

西安活动 | 2019年1月13号 拥抱开源, 又见.NET 线下交流活动报名进行中

随着.NET Core的发布和开源&#xff0c;.NET又重新回到人们的视野。.NET Core的下个3.0即将release&#xff0c;加入非常多的新功能&#xff0c;越来越拥抱变化&#xff0c;DevOps和Microservice的最佳实践已经在.NET Core落地&#xff0c;比如 Ocelot网关、GrpcConsul 服务注册…

[ZOJ 3203] 灯泡

今天这个**又带着题走来了&#xff01;&#xff01; 文章目录[TOC] 题目题解代码实现题目 相比 wildleopard 的家&#xff0c;他的弟弟 mildleopard 比较穷。他的房子是狭窄的而且在他的房间里面仅有一个灯泡。每天晚上&#xff0c;他徘徊在自己狭小的房子里&#xff0c;思考如…

YBTOJ:字符串题(KMP)

文章目录题目描述解析1.p[i]>0:2.p[i]0&#xff1a;代码题目描述 解析 我个人做起来很费劲的一道题&#xff0c;用vector&#xff0c;并查集等等东西搞来搞去过掉了&#xff08;竟然只WA了一次 &#xff09; 看题解思路就一下子清晰了&#xff0c;还是对KMP的理解不到位。 …

Asp.Net Core 轻松学-正确使用分布式缓存

前言本来昨天应该更新的&#xff0c;但是由于各种原因&#xff0c;抱歉&#xff0c;让追这个系列的朋友久等了。上一篇文章 在.Net Core 使用缓存和配置依赖策略 讲的是如何使用本地缓存&#xff0c;那么本篇文章就来了解一下如何使用分布式缓存&#xff0c;通过本章&#xff0…

计算几何初步

两点之间距离 欧氏距离 即欧几里得距离。 平面内两点的距离为 \[\sqrt{(x_1-x_2)^2(y_1-y_2)^2} \]立体空间内两点的距离为 \[\sqrt{(x_1-x_2)^2(y_1-y_2)^2(z_1-z_2)^2} \]\(\dots\) \(n\) 维空间内两点的距离为 \[\sqrt{\sum_{i1}^{n}{(x_1-x_2)^2}} \]曼哈顿距离 二维空间内…

树链剖分概念及模板 + 例题 [POJ3237 tree + 软件包管理器]

文章目录概念模板例题1&#xff1a;软件包管理器题目题解代码实现例题2&#xff1a;POJ3237 tree题目题解代码实现概念 树链剖分主要是用于解决以下这两个问题。 1、更改树上点x到点y的最短路径上的所有结点的值 2、查询树上点x到点y的最短路径上的所有结点的和。 在讲树链剖分…

YBTOJ:字符串匹配(KMP)

文章目录题目描述解析代码题目描述 解析 看了题解。。。 这题的关键在于可以变换匹配的一个充要条件&#xff1a; 每个字符与前一个相同字符的距离相同 这个搞出来之后就可以以它为关键字进行KMP了 注意&#xff01; 当与前一个字符的距离超过匹配长度时&#xff0c;是没有意义…

TransactionScope事务处理方法介绍及.NET Core中的注意事项

作者&#xff1a;依乐祝原文链接&#xff1a;https://www.cnblogs.com/yilezhu/p/10170712.html今天在写CzarCms的UnitOfWork的使用使用到了这个TransactionScope事务&#xff0c;因此对它进行了相关资料的查阅并记录如下&#xff0c;希望对大伙在.NET Core中使用有所帮助。写在…

YBTOJ:字符匹配(KMP)

文章目录题目描述解析代码题目描述 解析 显然应该要尝试套kmp的板子 关键是如何套 也就是那个判断匹配的条件是什么的问题 本题的关键是当kmp匹配时&#xff0c;匹配位之前的所有位大小关系的顺序都是匹配的&#xff0c;所以我们只需要看当前位即可 考虑对b预处理出3个数组&a…

历史上的今天(history)+ 勇者斗恶龙(dragon)

朋友们我来了&#xff0c;好久都没有更新了&#xff0c;手实在痒痒&#xff0c;不撸两道&#xff0c;内心过意不去 文章目录A&#xff1a;历史上的今天&#xff08;history&#xff09;题目题解代码实现B&#xff1a; 勇者斗恶龙&#xff08;dragon&#xff09;题目描述题解代…