hiho一下 第144周(机会渺茫)解题报告及拓展

最新推荐文章于 2018-07-17 09:18:57 发布
最新推荐文章于 2018-07-17 09:18:57 发布
阅读量158 收藏
点赞数
原文链接：http://www.cnblogs.com/cmyg/p/6683179.html
版权

题目1 : 机会渺茫

时间限制:5000ms
单点时限:1000ms
内存限制:256MB

描述

小Hi最近在追求一名学数学的女生小Z。小Z其实是想拒绝他的，但是找不到好的说辞，于是提出了这样的要求：对于给定的两个正整数N和M，小Hi随机选取一个N的约数N'，小Z随机选取一个M的约数M'，如果N'和M'相等，她就答应小Hi。

小Z让小Hi去编写这个随机程序，到时候她review过没有问题了就可以抽签了。但是小Hi写着写着，却越来越觉得机会渺茫。那么问题来了，小Hi能够追到小Z的几率是多少呢？

输入

每个输入文件仅包含单组测试数据。

每组测试数据的第一行为两个正整数N和M，意义如前文所述。

对于40%的数据，满足1<=N,M<=106

对于100%的数据，满足1<=N,M<=1012

输出

对于每组测试数据，输出两个互质的正整数A和B（以A分之B表示小Hi能够追到小Z的几率）。

样例输入 
3 2
样例输出 
4 1

 

Solution:

 

Way1:

  1 #include <stdio.h>
  2 #include <stdlib.h>
  3 #include <stdbool.h>
  4 
  5 long ans=0,maxn,zhi[100000];
  6 
  7 void Get_Zhi()
  8 {
  9     long i,j;
 10     bool vis[maxn+1];
 11     for (i=2;i<=maxn;i++)
 12         vis[i]=true;
 13     for (i=2;i<=maxn;i++)
 14     {
 15         if (vis[i])
 16         {
 17             ans++;
 18             zhi[ans]=i;
 19         }
 20         for (j=1;j<=ans;j++)
 21         {
 22             if (i*zhi[j]>maxn)
 23                 break;
 24             vis[i*zhi[j]]=false;
 25             if (i%zhi[j]==0)
 26                 break;
 27         }
 28     }
 29 }
 30 
 31 int main()
 32 {
 33     long long n,m,a,b,temp,x,y;
 34     long i,g;
 35     scanf("%lld%lld",&n,&m);
 36     if (n>m)
 37         maxn=(long)sqrt(n);
 38     else
 39         maxn=(long)sqrt(m);
 40     Get_Zhi();
 41     if (n>m)
 42     {
 43         a=n;
 44         b=m;
 45     }
 46     else
 47     {
 48         a=m;
 49         b=n;
 50     }
 51     while (b)
 52     {
 53         temp=b;
 54         b=a%b;
 55         a=temp;
 56     }
 57     x=1;
 58     y=1;
 59     //a
 60     if (a!=1)
 61     {
 62         for (i=1;i<=ans;i++)
 63             if (a%zhi[i]==0)
 64             {
 65                 g=1;
 66                 while (a%zhi[i]==0)
 67                 {
 68                     a/=zhi[i];
 69                     g++;
 70                 }
 71                 x*=g;
 72                 if (a==1)
 73                     break;
 74             }
 75         //a,n,m都有该质数，且系数都为1(若系数大于1，则相乘大于最大值)
 76         if (a!=1)
 77         {
 78             //x*2,y*4->x,y*2
 79             y*=2;
 80             n/=a;
 81             m/=a;
 82         }
 83     }
 84 
 85     //n
 86     if (n!=1)
 87     {
 88         for (i=1;i<=ans;i++)
 89             if (n%zhi[i]==0)
 90             {
 91                 g=1;
 92                 while (n%zhi[i]==0)
 93                 {
 94                     n/=zhi[i];
 95                     g++;
 96                 }
 97                 y*=g;
 98                 if (n==1)
 99                     break;
100             }
101         if (n!=1)
102             y*=2;
103     }
104 
105     //m
106     if (m!=1)
107     {
108         for (i=1;i<=ans;i++)
109             if (m%zhi[i]==0)
110             {
111                 g=1;
112                 while (m%zhi[i]==0)
113                 {
114                     m/=zhi[i];
115                     g++;
116                 }
117                 y*=g;
118                 if (m==1)
119                     break;
120             }
121         if (m!=1)
122             y*=2;
123     }
124 
125     a=y;
126     b=x;
127     while (b)
128     {
129         temp=b;
130         b=a%b;
131         a=temp;
132     }
133     x/=a;
134     y/=a;
135     printf("%lld %lld\n",y,x);
136     return 0;
137 }

 

Way2:

  1 #include <stdio.h>
  2 #include <stdlib.h>
  3 #include <stdbool.h>
  4 
  5 long ans=0,maxn,zhi[100000];
  6 
  7 //打质数表！
  8 
  9 void Get_Zhi()
 10 {
 11     long i,j;
 12     bool vis[maxn+1];
 13     for (i=2;i<=maxn;i++)
 14         vis[i]=true;
 15     for (i=2;i<=maxn;i++)
 16     {
 17         if (vis[i])
 18         {
 19             ans++;
 20             zhi[ans]=i;
 21         }
 22         for (j=1;j<=ans;j++)
 23         {
 24             if (i*zhi[j]>maxn)
 25                 break;
 26             vis[i*zhi[j]]=false;
 27             if (i%zhi[j]==0)
 28                 break;
 29         }
 30     }
 31 }
 32 
 33 int main()
 34 {
 35     long long n,m,y=1;
 36     long i,j,g1=0,g2=0,p[100],q[100],u[100],v[100];
 37     scanf("%lld%lld",&n,&m);
 38     if (n>m)
 39         maxn=(long)sqrt(n);
 40     else
 41         maxn=(long)sqrt(m);
 42     Get_Zhi();
 43 
 44     //n
 45     if (n!=1)
 46     {
 47         for (i=1;i<=ans;i++)
 48             if (n%zhi[i]==0)
 49             {
 50                 g1++;
 51                 p[g1]=i;
 52                 q[g1]=1;
 53                 while (n%zhi[i]==0)
 54                 {
 55                     n=n/zhi[i];
 56                     q[g1]++;
 57                 }
 58                 if (n==1)
 59                     break;
 60             }
 61     }
 62     if (n!=1)
 63         y*=2;
 64 
 65     //m
 66     if (m!=1)
 67     {
 68         for (i=1;i<=ans;i++)
 69             if (m%zhi[i]==0)
 70             {
 71                 g2++;
 72                 u[g2]=i;
 73                 v[g2]=1;
 74                 while (m%zhi[i]==0)
 75                 {
 76                     m/=zhi[i];
 77                     v[g2]++;
 78                 }
 79                 if (m==1)
 80                     break;
 81             }
 82     }
 83     if (m!=1 && m!=n)
 84         y*=2;
 85     i=1;
 86     j=1;
 87     while (i<=g1 && j<=g2)
 88     {
 89         if (p[i]<u[j])
 90         {
 91             y*=q[i];
 92             i++;
 93         }
 94         else if (u[j]<p[i])
 95         {
 96             y*=v[j];
 97             j++;
 98         }
 99         else
100         {
101             if (q[i]>v[j])
102                 y*=q[i];
103             else
104                 y*=v[j];
105             i++;
106             j++;
107         }
108     }
109     while (i<=g1)
110     {
111         y*=q[i];
112         i++;
113     }
114     while (j<=g2)
115     {
116         y*=v[j];
117         j++;
118     }
119     printf("%lld 1\n",y);
120     return 0;
121 }

 

Way3:

  1 #include <stdio.h>
  2 #include <stdlib.h>
  3 #include <stdbool.h>
  4 
  5 long ans=0,maxn,zhi[100000];
  6 
  7 void Get_Zhi()
  8 {
  9     long i,j;
 10     bool vis[maxn+1];
 11     for (i=1;i<=maxn;i++)
 12         vis[i]=true;
 13     for (i=2;i<=maxn;i++)
 14     {
 15         if (vis[i])
 16         {
 17             ans++;
 18             zhi[ans]=i;
 19         }
 20         for (j=1;j<=ans;j++)
 21         {
 22             if (i*zhi[j]>maxn)
 23                 break;
 24             vis[i*zhi[j]]=false;
 25             if (i%zhi[j]==0)
 26                 break;
 27         }
 28     }
 29 }
 30 
 31 int main()
 32 {
 33     long long n,m,a,b,temp,x,y;
 34     long g,i,r=0,t,p[100],q[100];
 35     scanf("%lld%lld",&n,&m);
 36     if (n>m)
 37         maxn=(long)sqrt(n);
 38     else
 39         maxn=(long)sqrt(m);
 40     Get_Zhi();
 41     if (n>m)
 42     {
 43         a=n;
 44         b=m;
 45     }
 46     else
 47     {
 48         a=m;
 49         b=n;
 50     }
 51     while (b)
 52     {
 53         temp=b;
 54         b=a%b;
 55         a=temp;
 56     }
 57     x=1;
 58     y=1;
 59 
 60     n/=a;
 61     m/=a;
 62 
 63     //a
 64     r=0;
 65     if (a!=1)
 66     {
 67     for (i=1;i<=ans;i++)
 68         if (a%zhi[i]==0)
 69         {
 70             r++;
 71             p[r]=i;
 72             q[r]=1;
 73             while (a%zhi[i]==0)
 74             {
 75                 a/=zhi[i];
 76                 q[r]++;
 77             }
 78             x*=q[r];
 79             if (a==1)
 80                 break;
 81         }
 82         //a,n,m都有该质数，且系数都为1(若系数大于1，则相乘大于最大值)
 83         //x*2,y*4->x,y*2
 84         if (a!=1)
 85             y*=2;
 86     }
 87     p[r+1]=1000000;
 88 
 89     //n
 90     if (n!=1)
 91     {
 92         t=1;
 93         for (i=1;i<=ans;i++)
 94             if (n%zhi[i]==0)
 95             {
 96                 g=1;
 97                 while (n%zhi[i]==0)
 98                 {
 99                     n/=zhi[i];
100                     g++;
101                 }
102                 //p[r+1]=1000000;
103                 while (p[t]<i)
104                 {
105                     y*=q[t];
106                     t++;
107                 }
108                 if (p[t]==i)
109                 {
110                     y*=(g+q[t]-1);
111                     t++;
112                 }
113                 else
114                     y*=g;
115                 if (n==1)
116                     break;
117             }
118         if (n!=1)
119             y*=2;
120         while (t<=r)
121         {
122             y*=q[t];
123             t++;
124         }
125     }
126     else
127         y*=x;
128 
129     //m
130     if (m!=1)
131     {
132         t=1;
133         for (i=1;i<=ans;i++)
134             if (m%zhi[i]==0)
135             {
136                 g=1;
137                 while (m%zhi[i]==0)
138                 {
139                     m/=zhi[i];
140                     g++;
141                 }
142                 //p[r+1]=1000000;
143                 while (p[t]<i)
144                 {
145                     y*=q[t];
146                     t++;
147                 }
148                 if (p[t]==i)
149                 {
150                     y*=(g+q[t]-1);
151                     t++;
152                 }
153                 else
154                     y*=g;
155                 if (m==1)
156                     break;
157             }
158         if (m!=1)
159             y*=2;
160         while (t<=r)
161         {
162             y*=q[t];
163             t++;
164         }
165     }
166     else
167         y*=x;
168 
169     a=y;
170     b=x;
171     while (b)
172     {
173         temp=b;
174         b=a%b;
175         a=temp;
176     }
177     x/=a;
178     y/=a;
179     printf("%lld %lld\n",y,x);
180     return 0;
181 }
182 /*
183 input:
184 16 8
185 output:
186 5 1
187 
188 input:
189 60 48
190 output:
191 20 1
192 */

Advance:

 1 #include <stdio.h>
 2 #include <stdlib.h>
 3 #include <stdbool.h>
 4 #define maxn 100000
 5 #define maxs 100000
 6 
 7 long zhi[maxn],z[maxs+1][12],g[maxs+1][12],ansz[maxs+1];
 8 
 9 void Get_Zhi()
10 {
11     long i,j,s,ans=0;
12     long long k;
13     bool vis[maxn+1];
14     for (i=2;i<=maxn;i++)
15         vis[i]=true;
16     for (i=2;i<=maxn;i++)
17     {
18         if (vis[i])
19         {
20             ans++;
21             zhi[ans]=i;
22             for (j=i;j<=maxs;j+=i)
23             {
24                 ansz[j]++;
25                 z[j][ansz[j]]=i;
26                 //为了求约数总数，比正确个数多1
27                 g[j][ansz[j]]=2;
28             }
29             k=(long long)i*i;
30             s=2;
31             while (k<=maxs)
32             {
33                 for (j=k;j<=maxs;j+=k)
34                     g[j][ansz[j]]++;
35                 k=(long long)k*i;
36                 s++;
37             }
38         }
39         for (j=1;j<=ans;j++)
40         {
41             if (i*zhi[j]>maxn)
42                 break;
43             vis[i*zhi[j]]=false;
44             if (i%zhi[j]==0)
45                 break;
46         }
47     }
48 }
49 
50 int main()
51 {
52     long i,j,m,n;
53     double result=0,re;
54     for (i=1;i<=maxs;i++)
55         ansz[i]=0;
56     Get_Zhi();
57     for (n=1;n<=maxs;n++)
58         for (m=1;m<=maxs;m++)
59         {
60             re=1.0;
61             i=1;
62             j=1;
63             while (i<=ansz[n] && j<=ansz[m])
64             {
65                 if (z[n][i]<z[m][j])
66                 {
67                     re/=g[n][i];
68                     i++;
69                 }
70                 else if (z[n][i]>z[m][j])
71                 {
72                     re/=g[m][j];
73                     j++;
74                 }
75                 else
76                 {
77                     if (g[n][i]>g[m][j])
78                         re/=g[n][i];
79                     else
80                         re/=g[m][j];
81                     i++;
82                     j++;
83                 }
84             }
85             while (i<=ansz[n])
86             {
87                 re/=g[n][i];
88                 i++;
89             }
90             while (j<=ansz[m])
91             {
92                 re/=g[m][j];
93                 j++;
94             }
95             result+=re;
96         }
97     printf("%lf\n",result*1.0/maxs/maxs);
98     return 0;
99 }

 

 

转载于:https://www.cnblogs.com/cmyg/p/6683179.html

weixin_30565327
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
hiho-pytorch-base
02-17
【标题】"hiho-pytorch-base"是一个与PyTorch相关的项目，它可能是为了提供一个基础框架或者学习资源，帮助用户快速上手PyTorch的使用。在PyTorch中，用户可以构建和训练深度学习模型，尤其适用于研究和实验环境中。...
hihocoder和leetcode-hiho-and-leet:hiho-and-leet
07-06
【标题】"hihocoder和leetcode-hiho-and-leet:hiho-and-leet" 提供的项目，主要关注的是两个知名的在线编程挑战平台——hihocoder和LeetCode。这两个平台都致力于帮助程序员提升编程技能，尤其是对于算法和数据结构...
HihoCode 1284 : 机会渺茫【约数个数】
DTL66的博客
11-01 445
#1284 : 机会渺茫 Time Limit:5000ms Case Time Limit:1000ms Memory Limit:256MB Description 小Hi最近在追求一名学数学的女生小Z。小Z其实是想拒绝他的，但是找不到好的说辞，于是提出了这样的要求：对于给定的两个正整数N和M，小Hi随机选取一个N的约数N'，小Z随机选取一
hihocoder-1284 机会渺茫(水题)
weixin_30237719的博客
03-27 115
机会渺茫 时间限制:5000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi最近在追求一名学数学的女生小Z。小Z其实是想拒绝他的，但是找不到好的说辞，于是提出了这样的要求：对于给定的两个正整数N和M，小Hi随机选取一个N的约数N'，小Z随机选取一个M的约数M'，如果N'和M'相等，她就答应小Hi。 小Z让小Hi去编写这个随机程序，到时候她review过没有问...
一生的追求
lvye1221的专栏
12-01 819
1 开场白 我来自湖南的一个小村庄，认真读书考大学是我最好的出路。中学的时候，我爸常和我说，“考不上，就在家种田玩泥巴哈；考上了，家里就是砸锅卖铁都会供你读”。 我也很努力，语数外、理化生、政史地等，每一科都非常努力地学习，顺利地考上大学。而我爸也兑现了他的承诺。 来到大学，学习计算机科学与技术专业，深知爸爸的不容易，以及对未来工作的担忧，我大学四年几乎没逃过课、没玩过游戏，没谈
【hihocoder】机会渺茫（map，离散化，概率计算，gcd）
better_space的博客
11-01 715
#1284 : 机会渺茫 时间限制:5000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi最近在追求一名学数学的女生小Z。小Z其实是想拒绝他的，但是找不到好的说辞，于是提出了这样的要求：对于给定的两个正整数N和M，小Hi随机选取一个N的约数N‘，小Z随机选取一个M的约数M‘，如果N‘和M‘相等，她就答应小Hi。 小Z让小Hi去编写这个随机程序，
leetcode答案-CodingProblemReview:这是一个个人leetcode，hiho，kikstart问题解决和文书工作库，
06-30
leetcode,hiho, kikstart problem solving and paper work repository After work like a year, I am thinking my objective is cloud computing, focusing on storage, buildout, deployment. It is a little bit ...
matlab代码做游戏-PAR-11-vs-HiHo-9-Scoring-for-Squash:蒙特卡洛模拟，观察统计数据，以将点对点计分系统
05-21
vs.HiHo-9得分 蒙特卡洛模拟，观察统计数据，以将点对点计分系统与MA Squash在2009年提出的上手交接式计分系统进行对比。 以下代码在计算赢得一场比赛的概率时假设一个简单的参数分布。 该剧本在两个玩家之间进行...
Paladin:准备工作面试。 包含文档，代码，博客，问题等
04-30
此项目是由zhuoyuan, Clint Wu和sosohu三人维护,主要是分享和交流各自的关于找工作相关的内容，包括面试题，解题报告，代码分享，理论知识等等. #主要文件夹 文件夹 内容 doc 主要存放一些文档类的东西，像一些经典...
[offer收割]编程练习赛3-题目2：机会渺茫
andahuzhuang的博客
03-28 502
题目2 : 机会渺茫 时间限制:5000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi最近在追求一名学数学的女生小Z。小Z其实是想拒绝他的，但是找不到好的说辞，于是提出了这样的要求：对于给定的两个正整数N和M，小Hi随机选取一个N的约数N'，小Z随机选取一个M的约数M'，如果N'和M'相等，她就答应小Hi。 小Z让小Hi去编写这
hihocoder offer收割赛。。#1284
Sun.StriKE()的博客
03-27 671
好久没刷题，水一水，反正排不上名次。。这道题记下 我想着蛋疼的做了质因数分解，然后再算的因子个数。。慢的一比，结果导致超时，还不如直接一个for循环搞定。。也是醉了 最后其实就是算两个数的gcd，然后gcd的因子个数除两个数的因子个数乘积就是最后结果，再约分一下就行。 #1284 : 机会渺茫 时间限制:5000ms 单点时限:1000ms 内存限制:256
[数论] 约数个数定理与约数和定理
Uniontake
04-05 5618
约数个数定理对于一个大于1正整数n可以分解质因数： 则n的正约数的个数就是  。其中a1、a2、a3…ak是p1、p2、p3，…pk的指数。约数定理证明首先同上，n可以分解质因数：n=p1^a1×p2^a2×p3^a3*…*pk^ak,由约数定义可知p1^a1的约数有:p1^0, p1^1, p1^2......p1^a1 ，共（a1+1）个;同理p2^a2的约数有（a2+1）个......pk^...
hihoCoder-week144:机会渺茫
trash
04-02 263
描述 小Hi最近在追求一名学数学的女生小Z。小Z其实是想拒绝他的，但是找不到好的说辞，于是提出了这样的要求：对于给定的两个正整数N和M，小Hi随机选取一个N的约数N'，小Z随机选取一个M的约数M'，如果N'和M'相等，她就答应小Hi。 小Z让小Hi去编写这个随机程序，到时候她review过没有问题了就可以抽签了。但是小Hi写着写着，却越来越觉得机会渺茫。那么问题来了，小Hi能够追到小Z的几率是
从零开始的算法题生活(三)
yiranblade的博客
07-17 209
引言 练习的有点频繁了，无奈近来心情太乱，刷刷算法题有益身体健康^_^,本次的题是hihoCoder144周练习题 描述 小Hi最近在追求一名学数学的女生小Z。小Z其实是想拒绝他的，但是找不到好的说辞，于是提出了这样的要求：对于给定的两个正整数N和M，小Hi随机选取一个N的约数N’，小Z随机选取一个M的约数M’，如果N’和M’相等，她就答应小Hi。 小Z让小Hi去编写这个随机程序，到时候她...
hihoCoder | hiho一下144周 机会渺茫
雄性皮卡丘
04-02 478
hiho一下144周题，题意是找出随机取出两个数刚好是n,m公因数的概率。思路是分别找出两个数所有因数，相乘得到取数情况，易得公因数个数。再调用求最大公因数函数，输出化简后的数字就行。map存放n的因数，m每次找到因数时到map去找是否已包含，以此确定公因数，将map理解为一个hashset即可。#include <iostream> #include <map> using namespace s
Hihocoder 147周 小HI的烦恼 bitset(技巧)
Masamiiiiii
04-25 381
点击打开链接 题意:n 暴力O(5*n^2) TLE 设 id[i][j] 第j门课排名为i的学生编号,s[j][i] 第j门课 排名为1..i的是那些人 第i人的答案为 第k门排名高于他的人的集合相交(k=1~5),则ans=s[k][a[i][k]-1]相交后集合中1的个数 #include using namespace std; const int N=3e4+10; in
hiho一下 第144周 机会渺茫
硕鼠的博客
04-01 539
题目1 : 机会渺茫 传送门 时间限制:5000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi最近在追求一名学数学的女生小Z。小Z其实是想拒绝他的，但是找不到好的说辞，于是提出了这样的要求：对于给定的两个正整数N和M，小Hi随机选取一个N的约数N'，小Z随机选取一个M的约数M'，如果N'和M'相等，她就答应小Hi。 小Z让小Hi
2024年公关服务行业分析报告.pptx
最新发布
07-09
行业报告

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值