内容 :

Hanks 博士是 BT (Bio-Tech,生物技术)  领域的知名专家,他的儿子名叫 Hankson。现在,刚刚放学回家的 Hankson 正在思考一个有趣的问题。 
     今天在课堂上,老师讲解了如何求两个正整数 c1和 c2 的最大公约数和最小公倍数。现在 Hankson 认为自己已经熟练地掌握了这些知识,他开始思考一个“求公约数”和“求公倍数”之类问题的“逆问题” ,这个问题是这样的:已知正整数 a0,a1,b0,b1,设某未知正整数 x 满足: 
1. x 和 a0 的最大公约数是 a1; 
2. x 和b0 的最小公倍数是 b1。 
     Hankson 的“逆问题”就是求出满足条件的正整数 x。但稍加思索之后,他发现这样的x 并不唯一,甚至可能不存在。因此他转而开始考虑如何求解满足条件的 x 的个数。请你帮助他编程求解这个问题。 

输入说明 :

第一行为一个正整数 n,表示有 n 组输入数据。接下来的 n 行每行一组输入数据,为四个正整数 a0,a1,b0,b1,每两个整数之间用一个空格隔开。输入数据保证 a0能被 a1 整除,b1 能被 b0整除。

输出说明 :

共n 行。每组输入数据的输出结果占一行,为一个整数。 
对于每组数据:若不存在这样的 x,请输出 0; 若存在这样的 x,请输出满足条件的 x 的个数。

范例输入 :

若題目沒有特別說明,則應該以多測資的方式讀取,若不知如何讀取請參考 a001 的範例程式。
241 1 96 28895 1 37 1776

范例输出 :

62

提示 :

第一组输入数据,x 可以是 9、18、36、72、144、288,共有6个。 
第二组输入数据,x 可以是 48、1776,共有 2 个。

对于 50%的数据,保证有 1≤a0,a1,b0,b1≤10000且 n≤100。 
对于 100%的数据,保证有 1≤a0,a1,b0,b1≤2,000,000,000且 n≤2000。

出处 :

NOIP2009提高组第二题 (管理員:liouzhou_101)
此题有几点要注意:
1.分解因数要注意效率,将条件设为:
while i<=sqrt(n)+1 do

而不是 while n<>1 do
然后在后面加上:
if n<>1 thenbegininc(count);a[count].data:=n;a[count].exp:=1;end;

效率大大增加。
2.数字相乘时要留心越界
求最小公倍数时,一定要写成
if (b0 div gcd(n,b0)*n)=b1 then f:=true;

而不是
if (b0*n div gcd(n,b0))=b1 then f:=true;

乍看好像没有差别,但要注意到后者的b0*n已超出longint的范围了!
程序:
program d778;
var a0,a1,b0,b1,i,count,x,m,k,res:longint;a:array[1..1000]of recorddata,exp:longint;end;b:array[0..1000000]of longint;
procedure deal(n:longint);
var i:longint;
begincount:=0;i:=2;while i<=sqrt(n)+1 dobeginif n mod i=0 thenbegininc(count);a[count].data:=i;a[count].exp:=0;while n mod i=0 dobegininc(a[count].exp);n:=n div i;end;end;inc(i);end;if n<>1 thenbegininc(count);a[count].data:=n;a[count].exp:=1;end;
end;
procedure dfs(n:longint);
var t,i:longint;
beginif n>count thenbegininc(b[0]);b[b[0]]:=x;endelsebeginif n=1 then x:=1;for i:=0 to a[n].exp dobegint:=x;dfs(n+1);x:=t*a[n].data;end;end;
end;
function gcd(a,b:longint):longint;
beginif a mod b=0 then gcd:=belse gcd:=gcd(b,a mod b);
end;
function f(n:longint):boolean;
var t:longint;
beginf:=false;t:=gcd(n,a0);if t=a1 thenif (b0 div gcd(n,b0)*n)=b1 then f:=true;
end;
beginreadln(m);for k:=1 to m dobeginfillchar(a,sizeof(a),0);readln(a0,a1,b0,b1);if b1=1 then begin count:=1;a[1].data:=1;a[1].exp:=0;end else deal(b1);b[0]:=0;dfs(1);res:=0;for i:=1 to b[0] do if f(b[i]) then inc(res);writeln(res);end;
end.

转载于:https://www.cnblogs.com/lixihan/archive/2010/09/17/1829612.html

d778: NOIP2009 2.Hankson的趣味题相关推荐

  1. Hankson的趣味题

    Hankson的趣味题 时间限制: 1 Sec  内存限制: 128 MB [提交][状态][讨论版] 题目描述 已知正整数a0,a1,b0,b1.设某未知正整数x 满足: 1. x 和a0 的最大公 ...

  2. 约数——Hankson的趣味题

    Hankson的趣味题 Hanks博士是BT(Bio-Tech,生物技术)领域的知名专家,他的儿子名叫Hankson. 现在,刚刚放学回家的Hankson正在思考一个有趣的问题. 今天在课堂上,老师讲 ...

  3. java蓝桥杯练习 Hankson的趣味题

    java蓝桥杯练习 Hankson的趣味题 资源限制 时间限制:1.0s 内存限制:64.0MB 问题描述 Hanks 博士是BT (Bio-Tech,生物技术) 领域的知名专家,他的儿子名叫Hank ...

  4. 洛谷1072 Hankson 的趣味题

    https://www.luogu.org/problem/show?pid=1072 Hankson 的趣味题 正解:素数筛+因数分解+乘法原理.简单地来说就是我不会! 但是考虑: 1.从b1入手, ...

  5. P1072 [NOIP2009 提高组] Hankson 的趣味题

    题目描述 Hanks 博士是 BT(Bio-Tech,生物技术) 领域的知名专家,他的儿子名叫 Hankson.现在,刚刚放学回家的 Hankson 正在思考一个有趣的问题. 今天在课堂上,老师讲解了 ...

  6. [NOIp 2009]Hankson的趣味题

    Description Hanks 博士是 BT (Bio-Tech,生物技术) 领域的知名专家,他的儿子名叫 Hankson.现在,刚刚放学回家的 Hankson 正在思考一个有趣的问题. 今天在课 ...

  7. 【noip】HankSon的趣味题

    描述 Hanks 博士是 BT (Bio-Tech,生物技术) 领域的知名专家,他的儿子名叫 Hankson.现 在,刚刚放学回家的 Hankson 正在思考一个有趣的问题. 今天在课堂上,老师讲解了 ...

  8. CODEVS 1172 Hankson 的趣味题

    题目描述 Description Hanks 博士是BT (Bio-Tech,生物技术) 领域的知名专家,他的儿子名叫Hankson.现 在,刚刚放学回家的Hankson 正在思考一个有趣的问题. 今 ...

  9. 2017.4.22 hankson的趣味题 思考记录

    这个题非常的不容易 ,第一遍交爆零了 首先根据最大公约数和最小公倍数可知: 这个数x一定包含a1的所有质因数 这个数x一定不包含a0/a1的所有质因数 这个数一定包含b1/b0的所有质因数 这个数可能 ...

最新文章

  1. Redis 笔记(10)— 发布订阅模式(发布订阅单个信道、订阅信道后的返回值分类、发布订阅多个信道)
  2. C#最小化到托盘+双击托盘恢复+禁止运行多个该程序
  3. mybatis count返回null_Mybatis属性示例-Properties的三种配置方式
  4. Centos7安装Oracle JDK
  5. ASP.NET之Ajax系列(三)
  6. QT int QString char转换
  7. Tina-TI 电路仿真 多参数交流传输特性分析 频率与相位增益的关系 选择控制对象
  8. 徐明星斥资超4亿香港买壳:史玉柱女儿现身股东行列
  9. java三角函数计算器_java 计算器代码能实现三角函数和阶乘功能
  10. 特斯拉为什么要“干掉”保险丝和继电器?
  11. win7 flash 提示加载失败
  12. WCDMA 基本信令流程
  13. Python入门教程(从小白到高级)(新版)
  14. awk&sed替换字符串
  15. Linux 误删文件恢复命令及方法!
  16. 车辆检测(视觉分类器训练)
  17. win10软件拒绝访问删不掉_win10结束进程时拒绝访问的处理办法
  18. ssti练习之[CSCCTF 2019 Qual]FlaskLight 1
  19. 无线通信网络优化的自动路测系统设计(Matlab代码实现)
  20. 傅里叶变换与拉普拉斯变换的概念理解

热门文章

  1. python numpy 数组元素周围_Python numpy数组元素属性及访问
  2. dbnetlib sqlserver不存在或拒绝访问_部署IIS+PHP+SQL server环境
  3. java字符后移_java把字符串参数往后移3位后输出
  4. idea 启动tomcat 工程_如何在IDEA中创建web项目并且部署到Tomcat中
  5. 全国计算机等级查询系统贵州省,贵州计算机等级考试成绩查询入口
  6. python智能工厂_智能工厂关键技术应用 第八讲 智能工厂的Python编程应用
  7. 北邮 鲁鹏老师 视频笔记
  8. 千万级别数据查询优化_从千万级数据查询来聊一聊索引结构和数据库原理
  9. mysql load 占位符,mysql:用不存在的数据的占位符初始化摘要表
  10. android dialog activity 窗口全透明,Android Dialog形式的Activity