传送门

文章目录

  • 题意:
  • 思路:

题意:


思路:

首先nnn很小的话可以暴力连边,让后染个色求一个颜色最多的即可。但是这个题显然不行,由于是三次方,所以考虑质因子入手。
首先很容易就能想到将所有的数的质因子的幂次模上333,之后他对应的乘起来为三次方数的数是唯一的。因为对于同一个质数来说,aaa的幂次为0,1,20,1,20,1,2对应的bbb的幂次为0,2,10,2,10,2,1,所以现在问题就变成了如何快速将所有数的幂次模333,并且求出来他对应的数,然后答案就是二者取一个maxmaxmax即可,因为二者一定是一个取一个不取,即染上不同的色。下面考虑如何将所有数幂次模333。
我们可以筛出来200020002000以内的质数的三次幂,让后每次遍历跑一遍即可。这样筛出来的数xxx的质因子幂次都<3<3<3,但是我们要求他对应的数怎么办呢?显然不能直接分解质因子,这样的复杂度还是(2e9)\sqrt{(2e9)}(2e9)​的,这里有一个巧妙的做法,就是对x∗xx*xx∗x再进行一次上面的分解,得出来的数即为xxx对应的数,因为x∗xx*xx∗x可以将原来111变成222,222变成4mod3=14\bmod 3=14mod3=1,正符合上面的规律。
复杂度约为O(n∗200)O(n*200)O(n∗200)。

// Problem: 牛牛的最大兴趣组
// Contest: NowCoder
// URL: https://ac.nowcoder.com/acm/contest/7604/C
// Memory Limit: 524288 MB
// Time Limit: 2000 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];
LL prime[N],cnt;
map<LL,int>mp,has;
bool st[N];void get_prime(int n) {for(int i=2;i<=n;i++) {if(!st[i]) {prime[++cnt]=1ll*i*i*i;for(int j=i+i;j<=n;j+=i)st[j]=true;}}
}LL divide(LL x) {for(int i=1;i<=cnt;i++) {if(x%prime[i]==0) {while(x%prime[i]==0) x/=prime[i];}if(prime[i]>x) return x;}return x;
}int main()
{//  ios::sync_with_stdio(false);
//  cin.tie(0);get_prime(1300);cout<<cnt<<endl;scanf("%d",&n);for(int i=1;i<=n;i++) {scanf("%d",&a[i]);a[i]=divide(a[i]);mp[a[i]]++;has[a[i]]=divide(1ll*a[i]*a[i]);}int ans=0;for(auto &x:mp) {if(x.X==1) {ans++;continue;}LL fa=has[x.X];ans+=max(mp[fa],x.Y);x.Y=0; mp[fa]=0;}cout<<ans<<endl;return 0;
}
/**/

nowcoder 牛牛的最大兴趣组 质因子 + 思维相关推荐

  1. 牛客 - 牛牛的最大兴趣组(思维+数论)

    题目链接:点击查看 题目大意:给出 nnn 个数,要求选出最多的数,使得任意两个数的乘积不能是三次平方数,三次平方数,诸如23=8,33=272^3=8,3^3=2723=8,33=27 题目分析:这 ...

  2. 质因子个数————质因子分解

    文章目录 1 题目 2 解析 3 参考代码 1 题目 质因子个数 时间限制 1000 ms 内存限制 32768 KB 代码长度限制 100 KB 判断程序 Standard (来自 小小) 题目描述 ...

  3. Codeforces Round #601 (Div. 2) E2. Send Boxes to Alice (Hard Version) 思维 + 质因子

    传送门 文章目录 题意: 思路: 题意: 大体题意跟easyeasyeasy版本差不多,就是hardhardhard版本的aaa范围更大.见这里Codeforces Round #601 (Div. ...

  4. 等式(分解质因子求因子个数)

    链接: https://www.nowcoder.com/acm/contest/90/F 来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 32768K,其他语言6553 ...

  5. 求n!中含有质因子p的个数

    定理:  中含有质因子p的个数为  ,其中  int cal(int n, int p) {int ans = 0;while (n != 0) {ans += n / p;n /= p; //相当与 ...

  6. 求一个数的因子个数/因子和/质因子 C/C++实现

    求一个数的因子个数时间复杂度O√n ll get_number(ll x){ll num=0;for(ll i=1;i*i<=x;i++){if(x%i==0) num+=2; if(i*i== ...

  7. 阶乘分解质因数[经典题组合数学枚举质因子]

    引入问题: 给定整数NNN,试把阶乘 N!N!N! 分解质因数,按照算术基本定理的形式输出分解结果中的 pip_ipi​ 和 cic_ici​即可. N!N!N!分解质因数后的结果,共若干行,每行一对 ...

  8. UVA-10212 The Last Non-zero Digit. 分解质因子+容斥定理

    这个是参考了别人之后的代码,POJ上0MS过了.Orz......对于一个序列在提取了2,5之后,例如1,2,3,4,5,6,7,8,9,10,我们可以将其中的奇数和偶数分开来对待,对于偶数序列2,4 ...

  9. PAT甲级1059 Prime Factors :[C++题解]分解质因子

    文章目录 题目分析 题目链接 题目分析 暴力求质因数 下面i就是质因子,s是质因子i的阶数. 暴力的时间复杂度O(n),会超时 void divide(int n){for(int i=2;i< ...

最新文章

  1. 专访王田苗:机器人是“刚需”,市场正处于逆周期增长
  2. 亚马逊贝索斯伸出橄榄枝后,巴菲特也力挺特朗普
  3. PX Deq: Execution Msg等待事件
  4. C#实现快速重命名JPEG文件
  5. 【Java 并发编程】线程池机制 ( ThreadPoolExecutor 线程池构造参数分析 | 核心线程数 | 最大线程数 | 非核心线程存活时间 | 任务阻塞队列 )
  6. 【2018/10/11完成】 ubuntu切换国内镜像源,加速apt-get下载速度
  7. mybatis plus条件拼接
  8. 在html中取消或加上超链接的下划线
  9. 使用码云git的webhook实现生产环境代码的自动pull
  10. python实现邮件客户端_利用python实现简单的邮件发送客户端示例
  11. pythontuple数据类型_Python数据类型之Tuple元组实例详解
  12. 解决办法:dpkg: 错误: 无法打开软件包的 info 文件 /var/lib/dpkg/available 以便读取: 没有那个文件或目录
  13. 深入浅出通信原理笔记(1)
  14. 大家好,我是谢公子,来自深信服—深蓝攻防实验室
  15. 微软笔试题 回忆(回文方面)
  16. Python制作安卓游戏外挂
  17. Jetson Nano 关闭开启图形界面减少内存占用
  18. TI达芬奇系列TMS320DM8148浮点DSP C674x + ARM Cortex-A8高性能视频处理器;
  19. 钉钉开放文档——JSAPI鉴权失败
  20. python生成关键词

热门文章

  1. ora22813操作数值超出系统的限制_最新:华为“鸿蒙”操作系统终于面世!一旦遭到限制,将随时启用...
  2. 王道408数据结构——第三章 栈和队列
  3. 编写程序在窗口中写出自己名字的拼音缩写_各类英文缩写:全称居然这么朴素?网友:最后一个我笑了...
  4. 上班骚扰男同事被抓到......
  5. 暴击!被初中生碾压智商!这份被国家数学集训队采用的初中奥数资料究竟有多厉害?...
  6. 法国一家夜总会上演机器人钢管舞
  7. 强烈推荐!孩子的科普从这套全球畅销250万册的最酷科学书起步
  8. 老板啥都懂,还天天套路我?!
  9. 12星座程序员写代码
  10. python双除_python2和python3 中的除法的区别