描述

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 的个数;

样例

样例输入

2
41 1 96 288
95 1 37 1776

样例输出

6
2

限制

每个测试点1s

来源

NOIP 2009

这个是2009年的题,我只想说,这道题虽然不是很难,但是真心恶心。
写法比较多,我这里写得比较长,网上有些写得特别短的,我蒟蒻比不上。
我这里就简单的说一下我的做法:
首先就是把200000000−−−−−−−−√\sqrt{200000000}内的质数筛出来,然后建立3个数组分别表示可以取的质因子,必须取的质因子,一定不能取的质因子都找出来(什么都给你了,随便乱搞就都出来了)
然后继续乱搞,比如某一个质因子同时存在于必须去和不能去,就爆0,然后如果可取和不能去同时在就不能去,就像这样把可以取的更新到最后,然后假设这个数组为p1,p2,p3,p4,pnp_1,p_2,p_3,p_4,p_n那么答案就是(p1+1)∗(p2+1)∗(p3+1)∗……∗(pn+1)(p_1+1)*(p_2+1)*(p_3+1)*……*(p_n+1)
就这样乱搞,就出来了
然后下面是我的代码,乱搞的太长了,仅供参考

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<set>
#include<queue>
#include<algorithm>
#include<vector>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<stack>
#define INF 2100000000
#define ll long long
#define clr(x)  memset(x,0,sizeof(x));
#define M 5000using namespace std;int pri[M],Eshai[45000],t;void get_pri()
{for(int i=2;i<=45000;i++)if(!Eshai[i])for(int j=i+i;j<=45000;j+=i)Eshai[j]=1;t=0;for(int i=2;i<=45000;i++)if(!Eshai[i])pri[++t]=i;
}void get_yin(int k,int num[])
{int i=1;while(pri[i]<=k){while(k%pri[i]==0){num[i]++;k/=pri[i];}i++;if(i>t){num[0]=k;return ;}}
}void get_may_must(int a,int b,int may[],int must[])
{int A[M],B[M];clr(A);clr(B);get_yin(a,A);get_yin(b,B);if(A[0]==B[0]&&A[0]!=0)may[0]=1;for(int i=1;i<=t;i++){if(B[i]>A[i])must[i]=B[i];if(A[i]==B[i])may[i]=A[i];}}void get_cant(int a,int b,int cant[])
{int A[M],B[M];clr(A);clr(B);get_yin(a,A);get_yin(b,B);for(int i=1;i<=t;i++)if(A[i]>B[i])cant[i]=B[i]+1;
}int n;int main()
{freopen("son.in","r",stdin);freopen("son.out","w",stdout);get_pri();cin>>n;for(int i=1;i<=n;i++){int a1,a2,b1,b2,bo=0;scanf("%d%d%d%d",&a1,&a2,&b1,&b2);int may[M],cant[M],gcd[M],must[M];clr(may);clr(cant);clr(gcd);clr(must)get_may_must(b1,b2,may,must);get_cant(a1,a2,cant);get_yin(a2,gcd);for(int i=1;i<=t;i++){if(must[i]>=cant[i]&&cant[i]){cout<<0<<'\n';bo=1;break;}if(cant[i])may[i]=0;if(gcd[i]<=may[i])may[i]-=gcd[i];}if(bo)continue;ll ans=1;for(int i=0;i<=t;i++)ans*=may[i]+1;printf("%d\n",ans);}return 0;
}

大概就是这个样子,如果有什么问题,或错误,请在评论区提出,谢谢。

【noip】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. [NOIp 2009]Hankson的趣味题

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

  6. d778: NOIP2009 2.Hankson的趣味题

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

  7. P1072 [NOIP2009 提高组] 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. Elasticsearch内存分配设置详解
  2. Android SDK Manager国内下载缓慢的问题
  3. 数据库-ADONET-在数据集DataSet中使用关系对象DataRelation处理关系
  4. Python 使用控制台运行带有相对路径的指令时,是以运行文件为基准,还是以控制台当前路径为基准
  5. Linux 字符设备驱动开发基础(一)—— 编写简单 LED 设备驱动
  6. Leetcode 82. Remove Duplicates from Sorted List II
  7. Python 3.8 稳定版正式发布,新特性全面解读
  8. 牛客网编程题04--字符串处理
  9. Linux 命令源码 —— cat
  10. php mysql sql model_PHP数据库模型
  11. 六石管理学:半弹性工作时间
  12. boxplot函数--Matplotlib
  13. KHV0031-himall3.0商城异常类(一)
  14. 忍者安全渗透系统(NINJITSU OS V3)的安装详细过程,亲测新旧vm版本都可安装,附带下载来源
  15. 泰拉服务器生存模式物品修改,荒岛求生修改存档生存物品方法
  16. 成都python数据分析师职业技能_合格大数据分析师应该具备的技能
  17. STM32——打地鼠
  18. 【codevs2980】买帽子 动态规划
  19. 爱普生Epson Expression Photo XP-850 一体机驱动
  20. WPF XAML X名称空间详解

热门文章

  1. java jsp中的日历表,jsp日历表格怎么做
  2. 关于谷歌搜索的大体_彭世瑜_新浪博客
  3. 所有Java程序员必须掌握的十款开源工具!
  4. PS学习之合成特效:被风沙侵蚀的动物们
  5. Flyway自学之路-04(springboot结合Flyway)
  6. XML与JSON(超级详解)
  7. 图文并茂:JVM内存布局彻底懂了!
  8. 百分号用计算机语言怎么表示,c++语言中 百分号如何输出 如要求输出50%
  9. 中国公司为 Wine ARM 实现申请专利并获批准
  10. OCPC系列 - OCPC介绍扫盲贴来啦