HDU 4542
T_T终于让我过了,坑啊,竟然时限是200ms。。。
我是预处理出不整除了个数的,因为这个较容易一点。利用算术基本定理,f=p1^a1*p2^a2......
所以,整除它的个数就是(a1+1)*(a2+1)......
开始预处理时,利用线性筛来先求素数,再计算,。。。呃,果断TLE了。后来发度娘告诉我一种令人发指的方法,那就是直接两个循环来求,比较令当前因数为i,循环求能被i整除的。。靠O(n^2)的啊,竟然比我的还要快。
另个求整除的,就直接DFS了。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define LL __int64
using namespace std;
const LL INF=(1ll<<62)+1;
//bool isprime[50005];
//int mprime[50005],cp=0;
int prime[20]={2,3,5,7,11,13,17,19,23,29,31,37,41,43,53,59};
//int cnt[50005],fac[100];
int findx[50005];
/*
void Init(){
memset(isprime,false,sizeof(isprime));
for(int i=2;i<=50000;i++){
if(!isprime[i]) mprime[cp++]=i;
for(int j=0;j<cp;j++){
if(mprime[j]>50000/i) break;
isprime[i*mprime[j]]=true;
if(i%mprime[j]==0) break;
}
}
}
/*
void WorkBefore(){
int facnt=0;
memset(cnt,0,sizeof(cnt));
memset(findx,0,sizeof(findx));
for(int i=2;i<=50000;i++){
facnt=0;
int ti=i;
for(int k=0;ti/mprime[k]>=mprime[k];k++){
if(ti%mprime[k]==0){
fac[facnt]=0;
while(ti%mprime[k]==0){
ti/=mprime[k];
fac[facnt]++;
}
facnt++;
}
}
if(ti!=1)
fac[facnt]=1,facnt++;
int e=1;
for(int k=0;k<facnt;k++){
e*=(fac[k]+1);
}
cnt[i]=i-e;
if(findx[cnt[i]]==0) findx[cnt[i]]=i;
}
}
*/ void dfs(int pos,int k,LL tans,LL &ans){
if(k==1&&tans<ans) ans=tans;
if(k<=1||tans>=ans) return;
LL t=1;
for(int i=1;i<=62;i++){
t*=prime[pos];
if(k<i+1) break;
if(k%(i+1)==0&&ans/t>=tans){
dfs(pos+1,k/(i+1),tans*t,ans);
}
}
} void Init(){
for(int i=1;i<=50000;i++) findx[i]=i;
for(int i=1;i<=50000;i++){
for(int j=i;j<=50000;j+=i) findx[j]--;
if(!findx[findx[i]]) findx[findx[i]]=i;
findx[i]=0;
}
} int main(){
Init();
// WorkBefore();
// cout<<"YES"<<endl;
int T,t=0,type,k;
scanf("%d",&T);
while(T--){
scanf("%d%d",&type,&k);
printf("Case %d: ",++t);
if(type==1){
int ans=-1;
if(findx[k]==0)
printf("Illegal\n");
else printf("%d\n",findx[k]);
}
else{
LL ans=INF;
dfs(0,k,1ll,ans);
if(ans==INF)
printf("INF\n");
else
printf("%I64d\n",ans);
}
}
return 0;
}
HDU 4542的更多相关文章
- hdu 4542 数论 + 约数个数相关 腾讯编程马拉松复赛
题目:http://acm.hdu.edu.cn/showproblem.php?pid=4542 小明系列故事--未知剩余系 Time Limit: 500/200 MS (Java/Others) ...
- hdu 4542 打表+含k个约数最小数
http://acm.hdu.edu.cn/showproblem.php?pid=4542 给出一个数K和两个操作 如果操作是0,就求出一个最小的正整数X,满足X的约数个数为K. 如果操作是1,就求 ...
- hdu 4542 小明系列故事——未知剩余系
小明系列故事——未知剩余系 题意:操作0表示某数有n个约数,操作1为某数有n个非约数:n <= 47777,若是存在小于2^62的数符合,则输出该数,否则若是不存在输出Illegal,若是大于2 ...
- hdu 4542 "小明系列故事——未知剩余系" (反素数+DFS剪枝)
传送门 参考资料: [1]:https://blog.csdn.net/acdreamers/article/details/25049767 题意: 输入两个数 type , k: ①type = ...
- hdu 4542 小明系列故事——未知剩余系 反素数 + 打表
小明系列故事——未知剩余系 Time Limit: 500/200 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others) Prob ...
- HDU 4542 小明系列故事——未知剩余系 (数论|反素数)
分析 kuangbin的blog已经讲的很好了,我做一点补充 1.当做x*y>z的比较时,如果x \(\ast\) y过大,可以写成x>z/y 2.分解质因数时选择用f[][0]保存质数, ...
- hdu4542 && ZOJ2562(反素数)
反素数: 对于任何正整数,其约数个数记为,例如,如果某个正整数满足:对任意的正整 数,都有,那么称为反素数. 有两个特点: 1.一个反素数的质因子必是从2开始的质数 2.如果,那么必有 最常见的问题如 ...
- HDOJ 2111. Saving HDU 贪心 结构体排序
Saving HDU Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total ...
- 【HDU 3037】Saving Beans Lucas定理模板
http://acm.hdu.edu.cn/showproblem.php?pid=3037 Lucas定理模板. 现在才写,noip滚粗前兆QAQ #include<cstdio> #i ...
随机推荐
- NetCore Netty 框架 BT.Netty.RPC 系列随讲 二 WHO AM I 之 NETTY/NETTY 与 网络通讯 IO 模型之关系?
一:NETTY 是什么? Netty 是什么? 这个问题其实百度上一搜一堆. 这是官方话的描述:Netty 是一个基于NIO的客户.服务器端编程框架,使用Netty 可以确保你快速和简单的开发出一个 ...
- 【洛谷3321_BZOJ3992】[SDOI2015]序列统计(原根_多项式)
题目: 洛谷3321 分析: 一个转化思路比较神(典型?)的题-- 一个比较显然的\(O(n^3)\)暴力是用\(f[i][j]\)表示选了\(i\)个数,当前积在模\(m\)意义下为\(j\)的方案 ...
- Patch 21352635 - Database Patch Set Update 11.2.0.4.8
一.CPU和PSU 近日,将数据库从9.2.0.6升级到11.2.0.4后,发现11.2.0.4通过DBLINK访问其他的9i库时发生ORA-02072错误,通过Google找到解决方案,即升级到PS ...
- 关于linux下的.a文件与 .so 文件
连续几天终于将一个又一个问题解决了,这里说其中一个问题 描述问题:使用多线程pthread的时候,(我用的IDE,CODEBOLCKS)编译后发现直接弹出窗口,程序还没有被Build..巴拉巴拉,然后 ...
- 检查阿里云ssl证书到期情况
#!/usr/bin/env python # -*- coding: utf-8 -*- # @Time : 2019-06-10 16:00 # @Author : Anthony.long # ...
- <stddef.h>
Common definitions 定义类型: ptrdiff_t 两指针相减的结果,signed integer size_t sizeof操作符的结果,unsigned integer max_ ...
- Android 检查手机上是否安装了指定的软件(根据包名检测)
Android检查手机上是否安装了指定的软件(根据包名检测) /** * 检查手机上是否安装了指定的软件 * @param context * @param packageName * @return ...
- spring 将配置文件中的值注入 属性
1.编写配置文件 #债权转让 #默认周期 必须大于0 credit.defaultDuration=1 #最小转让金额(元) credit.minBidAmount=1.00 #最小转让时间 到期时间 ...
- vue编辑回显问题
真是疯了,vue怪毛病真多 就下面这玩意儿,多选组合框,新增的时候好用的不行不行的,到了编辑的时候,要回显数据,怪毛病一堆一堆的 首先,回显的时候要传一个数组,但是这个数组里的元素得是字符串类型的,如 ...
- 使用form标签时注意事项
今天写程序的时候,使用了一个form:select标签,然后系统一直报错 原因找了好久也没找到 后来咨询得知, 在使用form:标签的时候 前后的form标签要写成<form:form> ...