uva 10061 How many zero's and how many digits ?(不同进制阶乘末尾几个0)+poj 1401

2024-09-09 16:48

本文主要是介绍uva 10061 How many zero's and how many digits ?(不同进制阶乘末尾几个0)+poj 1401,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题意是求在base进制下的 n!的结果有几位数,末尾有几个0。

想起刚开始的时候做的一道10进制下的n阶乘末尾有几个零,以及之前有做过的一道n阶乘的位数。

当时都是在10进制下的。


10进制下的做法是:

1. n阶位数:直接 lg(n!)就是得数的位数。

2. n阶末尾0的个数:由于2 * 5 将会在得数中以0的形式存在,所以计算2或者计算5,由于因子中出现5必然出现2,所以直接一趟遍历,找出能整除5的数就行了,5的个数就是末尾0的个数。


推广到base进制下:

1. n阶的位数:log base (n!),换底公式,可以换成 lg(n!) /  lg(base)。

2.n阶末尾0的个数:在n!中找可以凑成base的因子,比如16,可以找1, 16或者2,8或者4,4。

具体做法是n!里面的每一个数都分解成因子相乘的形式,并记录下来,然后到base中去凑,若能凑成,阶乘末尾就多一个0。


代码:

#include <stdio.h>
#include <math.h>
#include <string.h>int n, base;
int a[1000];//记录n中小于等于base的因子数int ZeroNum()
{memset(a, 0, sizeof(a));//将n!中的所有数拆分出因子,并记录下来该因子(j)对应有多少个for (int i = 2; i <= n; i++){int tmp = i;for (int j = 2; j <= tmp && j <= base; j++){while (tmp % j == 0){a[j]++;tmp /= j;}}}//将拆分出的因子枚举,若其能拼凑成一组base,则出现一个0。int res = 0;while (1){int tmp =  base;for (int i = 2; i <= tmp; i++){while (tmp % i == 0 && a[i] > 0){a[i]--;tmp /= i;}}if (tmp == 1)res++;elsebreak;}return res;
}int DigitNum()
{double sum = 0;int res;for (int i = 1; i <= n; i++){sum += log(i);}res = sum / log(base) + 1;return res;
}int main()
{#ifdef LOCALfreopen("in.txt", "r", stdin);#endif // LOCALwhile (scanf("%d%d", &n, &base) == 2){printf("%d %d\n", ZeroNum(), DigitNum());}return 0;
}

大神带你飞。

不理解这个求末尾0的做法:

add 2014.12.30:

正在学习数论,看完书再翻回来看这个代码就懂了。

这是素数的一个基本定理,即:n!的素因子分解种的素数p的幂为:(n / p) + (n / p ^ 2) + (n / p ^ 3) + …

在十进制下,能凑成末尾0的因子是2,5,而对于n!,在因式分解中,2的因子个数要大于5的因子个数(++先到2后到5,so有5必有2)。

所以如果存在一个因子5,必然对应着n!末尾的一个0,所以就变成了求n!中5的因子数,直接由性质得出,见下题poj1401.


转换成base进制下,自然成了找到乘积为base中最大的那个素因子,然后算有几个,最后转换下进制就行了。


#include <stdio.h>
#include <math.h>int main()
{#ifdef LOCALfreopen("in.txt", "r", stdin);#endif // LOCALint n, base;while (scanf("%d%d", &n, &base) == 2){int num;double sum = 0;//计算位数for (int i = 1; i <= n; i++){sum += log(i);}num = sum / log(base) + 1;//计算末位0的个数int cnt2 = 0;int cnt1;int maxprime;
//找最大素因子
for (int i = 2; i <= base; i++){cnt1 = 0;while (base % i == 0){maxprime = i;cnt1++;base /= i;}}while (n){n /= maxprime;cnt2 += n;}printf("%d %d\n", cnt2 / cnt1, num);}return 0;
}

poj 1401:

#include <stdio.h>
#include <math.h>int findZero(int n, int base)
{int cnt2 = 0;int cnt1;int maxprime;for (int i = 2; i <= base; i++){cnt1 = 0;while (base % i == 0){maxprime = i;cnt1++;base /= i;}}while (n){n /= maxprime;cnt2 += n;}return cnt2 / cnt1;
}int main()
{
#ifdef LOCAL//freopen("in.txt", "r", stdin);
#endif // LOCALint ncase;scanf("%d", &ncase);while (ncase--){int n;scanf("%d", &n);printf("%d\n", findZero(n, 10));}return 0;
}


这篇关于uva 10061 How many zero's and how many digits ?(不同进制阶乘末尾几个0)+poj 1401的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1151728

相关文章

每天认识几个maven依赖(ActiveMQ+activemq-jaxb+activesoap+activespace+adarwin)

八、ActiveMQ 1、是什么? ActiveMQ 是一个开源的消息中间件(Message Broker),由 Apache 软件基金会开发和维护。它实现了 Java 消息服务(Java Message Service, JMS)规范,并支持多种消息传递协议,包括 AMQP、MQTT 和 OpenWire 等。 2、有什么用? 可靠性:ActiveMQ 提供了消息持久性和事务支持,确保消

2. c#从不同cs的文件调用函数

1.文件目录如下: 2. Program.cs文件的主函数如下 using System;using System.Collections.Generic;using System.Linq;using System.Threading.Tasks;using System.Windows.Forms;namespace datasAnalysis{internal static

【Prometheus】PromQL向量匹配实现不同标签的向量数据进行运算

✨✨ 欢迎大家来到景天科技苑✨✨ 🎈🎈 养成好习惯,先赞后看哦~🎈🎈 🏆 作者简介:景天科技苑 🏆《头衔》:大厂架构师,华为云开发者社区专家博主,阿里云开发者社区专家博主,CSDN全栈领域优质创作者,掘金优秀博主,51CTO博客专家等。 🏆《博客》:Python全栈,前后端开发,小程序开发,人工智能,js逆向,App逆向,网络系统安全,数据分析,Django,fastapi

usaco 1.2 Palindromic Squares(进制转化)

考察进制转化 注意一些细节就可以了 直接上代码: /*ID: who jayLANG: C++TASK: palsquare*/#include<stdio.h>int x[20],xlen,y[20],ylen,B;void change(int n){int m;m=n;xlen=0;while(m){x[++xlen]=m%B;m/=B;}m=n*n;ylen=0;whi

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

uva 10055 uva 10071 uva 10300(水题两三道)

情歌两三首,水题两三道。 好久没敲代码了为暑假大作战热热身。 uva 10055 Hashmat the Brave Warrior 求俩数相减。 两个debug的地方,一个是longlong,一个是输入顺序。 代码: #include<stdio.h>int main(){long long a, b;//debugwhile(scanf("%lld%lld", &

hdu 2602 and poj 3624(01背包)

01背包的模板题。 hdu2602代码: #include<stdio.h>#include<string.h>const int MaxN = 1001;int max(int a, int b){return a > b ? a : b;}int w[MaxN];int v[MaxN];int dp[MaxN];int main(){int T;int N, V;s

poj 1511 Invitation Cards(spfa最短路)

题意是给你点与点之间的距离,求来回到点1的最短路中的边权和。 因为边很大,不能用原来的dijkstra什么的,所以用spfa来做。并且注意要用long long int 来存储。 稍微改了一下学长的模板。 stack stl 实现代码: #include<stdio.h>#include<stack>using namespace std;const int M

poj 3259 uva 558 Wormholes(bellman最短路负权回路判断)

poj 3259: 题意:John的农场里n块地,m条路连接两块地,w个虫洞,虫洞是一条单向路,不但会把你传送到目的地,而且时间会倒退Ts。 任务是求你会不会在从某块地出发后又回来,看到了离开之前的自己。 判断树中是否存在负权回路就ok了。 bellman代码: #include<stdio.h>const int MaxN = 501;//农场数const int

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n