[POI2010]Divine Divisor

题目大意:

给你\(m(m\le600)\)个数\(a_i(a_i\le10^{18})\)。\(n=\prod a_i\)。现在要你找到一个最大的\(k\)使得\(\exists d\ne1,d^k|n\),并求出有多少\(d\)满足这样的条件。

思路:

首先线性筛预处理出\(10^6\)以内的所有质数,用这些质数除\(a_i\),剩下的\(a_i\)分为以下\(4\)种情况:

  1. \(a_i=1\),表示\(a_i\)的所有素数均被找出。
  2. \(a_i=p^2\),可以判断\(\lfloor\sqrt{a_i}\rfloor\)是否等于\(\lceil\sqrt{a_i}\rceil\),是的话就说明这是两个\(>10^6\)的质数平方。
  3. \(a_i=p\),可以使用Miller-Rabin算法判断是否为质数。
  4. \(a_i=pq\),对于这样的数,可以与其它所有数求一遍\(\gcd\)。若\(\gcd\ne1\)就说明我们成功分解了它的质因数。否则虽然我们不能知道它的质因数到底是什么,但是我们可以知道它与其它数没有共同的质因数,因此我们只需要统计出现的次数,而不需要关心其具体数值。

对于每个质数,我们统计其出现次数\(cnt[i]\)。第一个答案就是\(\max\{cnt[i]\}\)。若有\(k\)个质数的出现次数为\(\max\{cnt[i]\}\),则第二个答案就是\(2^k-1\)。

由\(k\)可能会很大,需要写高精度。

但是我们可以注意到,若不考虑\(-1\),答案就是\(2\)的幂。用浮点数来储存不会丢失精度,且\(-1\)后不会发生退位。因此可以先用浮点数计算\(2^k\),转成字符串,再在最后一位\(-1\)。

源代码:

#include<map>
#include<cmath>
#include<ctime>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<cstdlib>
#include<algorithm>
typedef long long int64;
typedef __int128 int128;
inline int64 getint() {
register char ch;
while(!isdigit(ch=getchar()));
register int64 x=ch^'0';
while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
return x;
}
const int M=601,LIM=1e6+1,P=78499;
bool vis[LIM];
int p[P],b[M];
int64 a[M];
std::map<int64,int> cnt,cnt2;
inline void sieve() {
vis[1]=true;
for(register int i=2;i<LIM;i++) {
if(!vis[i]) p[++p[0]]=i;
for(register int j=1;j<=p[0]&&i*p[j]<LIM;j++) {
vis[i*p[j]]=true;
if(i%p[j]==0) break;
}
}
}
inline int64 montgomery(int64 a,int64 k,const int64 &mod) {
int64 ret=1;
for(;k;k>>=1) {
if(k&1) ret=(int128)ret*a%mod;
a=(int128)a*a%mod;
}
return ret;
}
inline bool miller_rabin(const int64 &x) {
for(register int i=0;i<5;i++) {
const int64 a=(int64)rand()*rand()%(x-2)+2;
if(montgomery(a,x-1,x)!=1) return false;
}
return true;
}
char ans[1000];
int main() {
sieve();
srand(time(NULL));
const int m=getint();
for(register int i=1;i<=m;i++) {
a[i]=getint();
for(register int j=1;j<=p[0]&&a[i]!=1;j++) {
while(a[i]%p[j]==0) {
a[i]/=p[j];
cnt[p[j]]++;
}
}
if(a[i]==1) continue;
if(floor(sqrt(a[i]))==ceil(sqrt(a[i]))) {
cnt[sqrt(a[i])]+=2;
b[i]=1;
continue;
}
if(miller_rabin(a[i])) {
cnt[a[i]]++;
b[i]=2;
continue;
}
}
for(register int i=1;i<=m;i++) {
if(a[i]==1||b[i]) continue;
for(register int j=1;j<=m;j++) {
if(a[i]==a[j]||a[j]==1) continue;
const int64 d=std::__gcd(a[i],a[j]);
if(d==1) continue;
cnt[d]++;
cnt[a[i]/d]++;
goto Next;
}
cnt2[a[i]]++;
Next:;
}
int ans1=0,ans2=0;
for(register std::map<int64,int>::iterator i=cnt.begin();i!=cnt.end();i++) {
ans1=std::max(ans1,i->second);
}
for(register std::map<int64,int>::iterator i=cnt2.begin();i!=cnt2.end();i++) {
ans1=std::max(ans1,i->second);
}
for(register std::map<int64,int>::iterator i=cnt.begin();i!=cnt.end();i++) {
if(i->second==ans1) ans2++;
}
for(register std::map<int64,int>::iterator i=cnt2.begin();i!=cnt2.end();i++) {
if(i->second==ans1) ans2+=2;
}
printf("%d\n",ans1);
sprintf(ans,"%.Lf",ldexpl(1,ans2));
ans[strlen(ans)-1]--;
puts(ans);
return 0;
}

[POI2010]Divine Divisor的更多相关文章

  1. BZOJ2082 : [Poi2010]Divine divisor

    将所有数分解质因数,那么第一问就是求指数的最大值,第二问就是$2^{指数最大的质数个数}-1$. 首先将$10^6$以内的质因数全部找到,那么剩下部分的因子$>10^6$,且只有3种情况: 1. ...

  2. 【BZOJ2082】【POI2010】Divine divisor 假的pollard-rho

    题目大意:给你$m$个数$a_i$,定义$n=\Pi_{i=1}^{m}a_i$.将$n$分解质因数为$\Pi p_i^{k_i} $,$p_i$是质数.请输出$2^{max(k_i)}-1$,以及存 ...

  3. POI2010题解

    POI2010题解 我也不知道我为什么就开始刷POI了 有些题目咕掉了所以不完整(我都不知道POI到底有多少题) [BZOJ2079][Poi2010]Guilds (貌似bz跟洛谷上的不是一个题?) ...

  4. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  5. [UCSD白板题] Greatest Common Divisor

    Problem Introduction The greatest common divisor \(GCD(a, b)\) of two non-negative integers \(a\) an ...

  6. BZOJ 2080: [Poi2010]Railway 双栈排序

    2080: [Poi2010]Railway Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 140  Solved: 35[Submit][Statu ...

  7. 【BZOJ】【2084】【POI2010】Antisymmetry

    Manacher算法 啊……Manacher修改一下就好啦~蛮水的…… Manacher原本是找首尾相同的子串,即回文串,我们这里是要找对应位置不同的“反回文串”(反对称?233) 长度为奇数的肯定不 ...

  8. BZOJ2086: [Poi2010]Blocks

    题解: 想了想发现只需要求出最长的一段平均值>k即可. 平均值的问题给每个数减去k,判断是否连续的一段>0即可. 然后我们发现如果i<j 且 s[i]<s[j],那么 j 对于 ...

  9. BZOJ2083: [Poi2010]Intelligence test

    2083: [Poi2010]Intelligence test Time Limit: 10 Sec  Memory Limit: 259 MBSubmit: 241  Solved: 96[Sub ...

随机推荐

  1. javascript中=、==与===的区别

    1.等号 =赋值运算符,给变量赋值 var a="1"; 2.相等和不相等操作符 相等操作符由==表示,若两个操作数相等,则返回true:不相等操作符由!=表示,若两个操作数不相等 ...

  2. Object的公用方法们

    如图所示,Object一共有10种方法: 下面详细描述: 1.public Object() 方法,默认构造函数方法,当新建一个Object对象的时候,调用这个方法向堆区申请一片内存: 2.priva ...

  3. MinnowBoard

    MinnowBoard https://github.com/RafaelRMachado/MinnowBoard https://github.com/RafaelRMachado https:// ...

  4. C基础 常用设计模式粗解

    引言 面向对象, 设计模式是现代软件开发基石. C的面向过程已经很简洁, 但不代表C就没有面向对象.(libuv框架中C面向对象用的很多) 因为思想是互通的.全当熟悉一下那些常用的设计模式.先假定有一 ...

  5. 2.C 基础

    C 基础 原文地址:http://rypress.com/tutorials/objective-c/c-basics OC 可以说是C语言的一个超集,这样你可以无缝的和C语言结合编程也就是你可以这两 ...

  6. [ python ] 类中的一些特殊方法

    item系列 __getitem__(self, item) 对象通过 object[key] 触发 __setitem__(self, key, value) 对象通过 object[key] = ...

  7. Jmeter接口测试示例

    如果是Web,需要使用badboy进行录制,今天讲的是接口,因此可以不用录制. (1)新建测试计划 (2)添加http请求默认值 (3)添加http信息头管理器 (4)添加token的正则表达式:&q ...

  8. P2327

    code[class*="language-"] { padding: .1em; border-radius: .3em; white-space: normal; backgr ...

  9. bzoj 2434 AC自动机 + fail指针建树 + 树状数组

    思路:我们先跟着它给定的字符串走把字典树建出来,求出fail指针,我们考虑两个字符串 A和B, 如果想要求B中有多少A的子串,转换一下就是有多少个B的前缀的后缀包含A,这个在AC自动机 的状态图中很容 ...

  10. 35.Spark系统运行内幕机制循环流程

    一:TaskScheduler原理解密 1,  DAGScheduler在提交TaskSet给底层调度器的时候是面向接口TaskScheduler的,这符合面向对象中依赖抽象而不依赖的原则,带来底层资 ...