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的更多相关文章

  1. hdu 4542 数论 + 约数个数相关 腾讯编程马拉松复赛

    题目:http://acm.hdu.edu.cn/showproblem.php?pid=4542 小明系列故事--未知剩余系 Time Limit: 500/200 MS (Java/Others) ...

  2. hdu 4542 打表+含k个约数最小数

    http://acm.hdu.edu.cn/showproblem.php?pid=4542 给出一个数K和两个操作 如果操作是0,就求出一个最小的正整数X,满足X的约数个数为K. 如果操作是1,就求 ...

  3. hdu 4542 小明系列故事——未知剩余系

    小明系列故事——未知剩余系 题意:操作0表示某数有n个约数,操作1为某数有n个非约数:n <= 47777,若是存在小于2^62的数符合,则输出该数,否则若是不存在输出Illegal,若是大于2 ...

  4. hdu 4542 "小明系列故事——未知剩余系" (反素数+DFS剪枝)

    传送门 参考资料: [1]:https://blog.csdn.net/acdreamers/article/details/25049767 题意: 输入两个数 type , k: ①type = ...

  5. hdu 4542 小明系列故事——未知剩余系 反素数 + 打表

    小明系列故事——未知剩余系 Time Limit: 500/200 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) Prob ...

  6. HDU 4542 小明系列故事——未知剩余系 (数论|反素数)

    分析 kuangbin的blog已经讲的很好了,我做一点补充 1.当做x*y>z的比较时,如果x \(\ast\) y过大,可以写成x>z/y 2.分解质因数时选择用f[][0]保存质数, ...

  7. hdu4542 && ZOJ2562(反素数)

    反素数: 对于任何正整数,其约数个数记为,例如,如果某个正整数满足:对任意的正整 数,都有,那么称为反素数. 有两个特点: 1.一个反素数的质因子必是从2开始的质数 2.如果,那么必有 最常见的问题如 ...

  8. HDOJ 2111. Saving HDU 贪心 结构体排序

    Saving HDU Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total ...

  9. 【HDU 3037】Saving Beans Lucas定理模板

    http://acm.hdu.edu.cn/showproblem.php?pid=3037 Lucas定理模板. 现在才写,noip滚粗前兆QAQ #include<cstdio> #i ...

随机推荐

  1. Windows(7/8/10)搭建kibana 6.x版本(elasticsearch的可视化服务)

    在搭建kibana之前,我们先了解下什么是kibana Kibana 是一款开源的数据分析和可视化平台,它是 Elastic Stack 成员之一,设计用于和 Elasticsearch 协作.您可以 ...

  2. 拼接sql ()

    SELECT ID,FORMNAME,NODENAME,SEQUENCE, NAME, STATE, NOWTIMES,      WORK.FQREALNAME||'('||FQDEPT.FULLN ...

  3. [Luogu 1312] noip11 Mayan游戏

    [Luogu 1312] noip11 Mayan游戏 Problem: Mayan puzzle是最近流行起来的一个游戏.游戏界面是一个 7 行5 列的棋盘,上面堆放着一些方块,方块不能悬空堆放,即 ...

  4. python网络爬虫数据中的三种数据解析方式

    一.正则解析 常用正则表达式回顾: 单字符: . : 除换行以外所有字符 [] :[aoe] [a-w] 匹配集合中任意一个字符 \d :数字 [0-9] \D : 非数字 \w :数字.字母.下划线 ...

  5. Android:EditText属性大全

    一.inputType属性inputType属性在EditText输入值时启动的虚拟键盘的风格有着重要的作用.比如有时需要虚拟键盘只为字符或只为数字. <span style="fon ...

  6. android UI 操作 不要在子线程中操作UI

    不管是android ,还是 ios ,请不要在子线程中操作UI,有时有些崩溃,从报错上看不出什么原因,就有可能是子线程操作了UI:切记,切记! 请放在主线程例: activity.runOnUiTh ...

  7. 注解配置AOP切面编程

    1.导入先关jar包 2.编写applicationContext.xml,配置开启注解扫描和切面注解扫描 <?xml version="1.0" encoding=&quo ...

  8. [转]c# 对密码执行散列和 salt 运算方法

    本文转自:http://www.cnblogs.com/CnBlogFounder/archive/2008/07/04/1235690.html 大家对密码执行散列和Salt运算一定不陌生.两个Vi ...

  9. centos7安装mysql和mysql-connector-c++

    最近为了搭建自己的开发环境,又一次在centos7上面开始安装mysql和c++的访问环境,特此记录一下搭建过程,方便以后查阅 一.安装mysql centos7 默认安装了mariaDB,导致不能安 ...

  10. html5——颜色

    CSS2 1.opacity,可以设置透明度,但是父盒子设置了透明度会影响子盒子 CC3 1.transparent属性,但是不可改变透明值 2.rgba():r--red g--green b--b ...