传送门

文章目录

  • 题意:
  • 思路:

题意:

给你nnn个数,每个数的因子个数不超过777个,选出最少的数使其乘积为平方数。
n≤1e5n\le 1e5n≤1e5

思路:

由于因子不超过777个,所以由约数个数(1+p1)∗(1+p2)∗...∗(1+pn)(1+p_1)*(1+p_2)*...*(1+p_n)(1+p1​)∗(1+p2​)∗...∗(1+pn​)得其质因子个数不超过222个,启发我们用质因子来做这个题。
一个显然的性质就是平方数的质因子幂次mod2\bmod 2mod2之后这个平方数为111。所以我们将每个数的质因子幂次mod2\bmod 2mod2之后留下为幂次111的质因子。现在问题就变成了选出尽可能少的点使其乘积的质因子幂次为偶数,考虑将其转换成图上的问题。
如果一个数本身就是平方数那么直接输出111就好啦,下面讨论质因子个数不为000的情况。
先考虑每个数的质因子都有两个的情况,考虑以质数为每个点,将每个数的两个质因子连边,那么画出来图之后可以发现最小环即为答案。
如果质因子只有一个怎么办呢?考虑建立一个虚拟节点111,将其与111连边,之后跑最小环就好啦。
那么问题来了,最小环?FloydFloydFloyd?显然是不可以的,我们需要发现我们构建的图的一些特殊性质。
首先边权都为111且为无向图,那么我们可以枚举起点,遍历每条边,记一个depthdepthdepth,每当碰到的depthdepthdepth被更新过了,那么ans=min(ans,depth[u]+depth[v]+1)ans=min(ans,depth[u]+depth[v]+1)ans=min(ans,depth[u]+depth[v]+1)即可,但是复杂度是n2n^2n2的仍然不可接受。
继续考虑优化起点,由于>1000>1000>1000的质数不会相互连边,所以只需枚举≤1000\le1000≤1000的质数即可,复杂度nan\sqrt ana​,可以接受。

// Problem: E. Ehab's REAL Number Theory Problem
// Contest: Codeforces - Codeforces Round #628 (Div. 2)
// URL: https://codeforces.com/contest/1325/problem/E
// Memory Limit: 256 MB
// Time Limit: 3000 ms
//
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=1000010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n;
int a[N];
int e[N],ne[N],h[N],idx;
int prime[N],cnt,ans=INF;
int depth[N];
bool st[N];
vector<int>v;void add(int a,int b) {e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}void divide(int x) {int a,b; a=b=-1;for(int i=2;i<=x/i;i++) {if(x%i==0) {int cnt=0;while(x%i==0) x/=i,cnt++;if(cnt%2==0) continue;if(a==-1) a=i;else if(b==-1) b=i;}}if(x>1) {if(a==-1) a=x;else b=x;}if(a==-1) {puts("1");exit(0);}if(b==-1) add(1,a),add(a,1);else add(a,b),add(b,a);
}void get_prime(int n) {for(int i=2;i<=n;i++) if(!st[i]) {for(int j=i+i;j<=n;j+=i) st[j]=1;prime[++cnt]=i;}
}void bfs(int st) {queue<int>q1,q2;memset(depth,-1,sizeof(depth));q1.push(st); q2.push(-1);depth[st]=0;while(q1.size()) {int a=q1.front(); q1.pop();int b=q2.front(); q2.pop();for(int i=h[a];~i;i=ne[i]) {int j=e[i];if((i^1)==b) continue;if(depth[j]==-1) depth[j]=depth[a]+1,q1.push(j),q2.push(i);else ans=min(ans,depth[a]+depth[j]+1);}}
}int main()
{//  ios::sync_with_stdio(false);
//  cin.tie(0);get_prime(2000);memset(h,-1,sizeof(h));scanf("%d",&n);for(int i=1;i<=n;i++) {scanf("%d",&a[i]);divide(a[i]);}for(int i=1;i<=cnt;i++) bfs(i);printf("%d\n",ans==INF? -1:ans);return 0;
}
/**/

Codeforces Round #628 (Div. 2) E. Ehab‘s REAL Number Theory Problem 巧妙的质因子建图相关推荐

  1. Codeforces Round #628 (Div. 2) F. Ehab‘s Last Theorem dfs树

    传送门 文章目录 题意: 思路: 题意: 给你个nnn个点mmm条边的图,可以选择完成以下两个任务中的一个: (1)(1)(1)找出大小恰好为n\sqrt nn​的一个独立集. (2)(2)(2)找出 ...

  2. 【Codeforces Round #525(Div. 2)】Ehab and another another xor problem(思维+异或)

    题目链接 D. Ehab and another another xor problem time limit per test 1 second memory limit per test 256 ...

  3. Codeforces Round #628 (Div. 2)C - Ehab and Path-etic MEXs

    C - Ehab and Path-etic MEXs 题目链接 简要描述: 从一个点到另一个点出发,所经过的边组成一个集合, 这些集合中所不包含的最小自然数为mex(u,v),要求使所有的mex(u ...

  4. Codeforces Round #628 (Div. 2) C. Ehab and Path-etic MEXs

    C. Ehab and Path-etic MEXs 题目链接-C. Ehab and Path-etic MEXs Note The tree from the second sample: 题目大 ...

  5. Codeforces Round #628 (Div. 2) D. Ehab the Xorcist

    D. Ehab the Xorcist 题目链接-D. Ehab the Xorcist 题目大意 给定两个数字 u 和 v ,构造一个最短的序列a满足a1⊕a2⊕-⊕an=u 并且 a1+a2+-+ ...

  6. Codeforces Round #628 (Div. 2) C. Ehab and Path-etic MEXs(贪心+思维)

    题目链接 思路:这个题做法貌似很多,我用的就是从叶子节点开始赋值,不过听大佬说的从度大于2的节点赋值也可以.比赛的时候wa test19了,特判一下n=2就行了. #include<bits/s ...

  7. Codeforces Round #628 (Div.2) C.Ehab and Path-etic MEXs(树,思维)

    传送门 题意: 给一颗n个结点的数,然后n-1条边,我们要做的就是把0-n-2,这n-1个数赋给n-1条边,然后使得所有MEX(u,v)最大值最小,输出每条边赋的值 MEX(u,v)是u到v这条路径上 ...

  8. Codeforces Round #525 (Div. 2) D Ehab and another another xor problem

    题目:http://codeforces.com/contest/1088/problem/D 题目大意:有0 <= a , b <= 230的a,b两个整数,定义一次操作(c , d)的 ...

  9. Codeforces Round #649 (Div. 2)C. Ehab and Prefix MEXs[排列的构造]

    C. Ehab and Prefix MEXs 题目大意: 解题思路:题目说保证a数组是非递减的,那么如果某位置a[i]!=a[i−1]a[i]!=a[i-1]a[i]!=a[i−1]那么这个位置ii ...

最新文章

  1. 自动驾驶软件工程课程之SLAM(2)
  2. 我的 Windows Vista 上安装的工具
  3. bluez 设置绑定pin码_国家工信部紧急提醒:一定要设置这个密码!
  4. linux目录空间内存,Linux 目录结构:内存文件夹
  5. 【转】java string类的方法及说明
  6. 文献学习(part76-A)--To cluster, or not to cluster: An analysis of clusterability methods
  7. C# - 为引用类型重定义相等性 - 继承相关
  8. 【汇编语言】8086汇编语言的debug中,t和p命令的区别
  9. 微小宝公众号排行榜_榜单 广东省技工院校微信公众号排行榜(第51期)
  10. 研究生从事学术研究的20条法则:导师错的时候不多!
  11. 前端开发 容器标签 0226
  12. 9行代码提高少样本学习泛化能力! ICLR2021 Oral,代码已开源
  13. CentOS 6.9 下安装DB2
  14. mysql 能用compute_SQL compute by 的使用
  15. 数据分析之数据透视表(Office2019)
  16. ArcGIS Pro与ArcGis区别
  17. 部署ROS2 Bouncy版本时遇到的一些问题
  18. 阿里云服务器学生应该怎么购买教程以及将远程主机配置成服务器教程
  19. 100个软硬件都要懂的示波器基础知识
  20. PHPBB网站如何搬迁新服务器,phpbb 安装教程

热门文章

  1. 中国古代数学有多牛?
  2. 男朋友和女朋友的区别。。。
  3. 全网爆红!全新洗脑神曲:程序员Disco
  4. 湘乡江南计算机学校,湘乡职业中等专业学校2021年招生录取分数线
  5. python tkinter进度条_在python3.7中更新tkinter进度条
  6. mysql isam 食物_MySQL常用存储引擎之MyISAM
  7. html li 做瀑布流,js实现瀑布流效果(自动生成新的内容)
  8. mysql group 条件,mysql - mysql group by,两个条件,限制1 - SO中文参考 - www.soinside.com...
  9. java gzip 解压文件_Java实现文件压缩与解压[zip格式,gzip格式]
  10. 7-6 0-1背包 (20 分)(思路加详解+网格做法+动态规划)Come Baby !!!!!!!!!!!!!!