分解质因数c语言设计思路,深入分析C语言分解质因数的实现方法
首先来看一个最简单的C语言实现质因数分解的列子:
#include
void main( )
{
int data, i = 2;
scanf("%d", &data);
while(data > 1)
{
if(data % i == 0)
{
printf("%d ", i);
data /= i;
}
else i++;
}
}
原理&&方法把一个合数分解为若干个质因数的乘积的形式,即求质因数的过程叫做分解质因数,分解质因数只针对合数
求一个数分解质因数,要从最小的质数除起,一直除到结果为质数为止。分解质因数的算式的叫短除法,和除法的性质差不多,还可以用来求多个个数的公因式:
以24为例:
2 -- 24
2 -- 12
2 -- 6
3 (3是质数,结束)
得出 24 = 2 × 2 × 2 × 3 = 2^3 * 3
代码
可先用素数筛选法,筛选出符合条件的质因数,然后for循环遍历即可,通过一道题目来show一下这部分代码
题目1
题目描述:
求正整数N(N>1)的质因数的个数。
相同的质因数需要重复计算。如120=2*2*2*3*5,共有5个质因数。
输入:
可能有多组测试数据,每组测试数据的输入是一个正整数N,(1
输出:
对于每组数据,输出N的质因数的个数。
样例输入:
120
样例输出:
5
提示:
注意:1不是N的质因数;若N为质数,N是N的质因数。
ac代码
#include
int main()
{
int n, count, i;
while (scanf("%d", &n) != EOF) {
count = 0;
for (i = 2; i * i <= n; i ++) {
if(n % i == 0) {
while (n % i == 0) {
count ++;
n /= i;
}
}
}
if (n > 1) {
count ++;
}
printf("%d\n", count);
}
return 0;
}
深入理解
我所谓的深入理解,就是通过4星的题目来灵活运用分解质因数的方法,题目如下
题目2
题目描述:
给定n,a求最大的k,使n!可以被a^k整除但不能被a^(k+1)整除。
输入:
两个整数n(2<=n<=1000),a(2<=a<=1000)
输出:
一个整数.
样例输入:
6 10
样例输出:
1
思路a^k和n!都可能非常大,甚至超过long long int的表示范围,所以也就不能直接用取余操作判断它们之间是否存在整除关系,因此我们需要换一种思路,从分解质因数入手,假设两个数a和b:
a = p1^e1 * p2^e2 * ... * pn^en, b = p1^d1 * p2^d2 * ... * pn^dn
, 则b除以a可以表示为:
b / a = (p1^d1 * p2^d2 * ... * pn^dn) / (p1^e1 * p2^e2 * ... * pn^en)
若b能被a整除,则 b / a必为整数,且两个素数必护质,则我们可以得出如下规律:
若a存在质因数px,则b必也存在该质因数,且该素因数在b中对应的幂指数必不小于在a中的幂指数
另b = n!, a^k = p1^ke1 * p2^ke2 * ... * pn^ken,因此我们需要确定最大的非负整数k即可。要求得该k,我们只需要依次测试a中每一个素因数,确定b中该素因数是a中该素因数的幂指数的多少倍即可,所有倍数中最小的那个即为我们要求得的k
分析到这里,剩下的工作似乎只是对a和n!分解质因数,但是将n!计算出来再分解质因数,这样n!数值太大。考虑n!中含有素因数p的个数,即确定素因数p对应的幂指数。我们知道n!包含了从1到n区间所有整数的乘积, 这些乘积中每一个p的倍数(包括其本身)都对n!贡献至少一个p因子,且我们知道在1到n中p的倍数共有n/p个。同理,计算p^2,p^3,...即可
代码
#include
#include
#include
#define N 1001
int prime[N], size;
/**
* 素数筛选法进行预处理
*/
void initProcess()
{
int i, j;
for (prime[0] = prime[1] = 0, i = 2; i < N; i ++) {
prime[i] = 1;
}
size = 0;
for (i = 2; i < N; i ++) {
if (prime[i]) {
size ++;
for (j = 2 * i; j < N; j += i) {
prime[j] = 0;
}
}
}
}
int main(void)
{
int i, n, a, k, num, count, base, tmp, *ansbase, *ansnum;
// 预处理
initProcess();
while (scanf("%d %d", &n, &a) != EOF) {
ansbase = (int *)calloc(size, sizeof(int));
ansnum = (int *)calloc(size, sizeof(int));
// 将a分解质因数
for (i = 2, num = 0; i < N && a != 1; i ++) {
if (prime[i] && a % i == 0) {
ansbase[num] = i;
ansnum[num] = 0;
while (a != 1 && a % i == 0) {
ansnum[num] += 1;
a = a / i;
}
num ++;
}
}
// 求最小的k
for (i = 0, k = 0x7fffffff; i < num; i ++) {
base = ansbase[i];
count = 0;
while (base <= n) {
count += n / base;
base *= ansbase[i];
}
tmp = count / ansnum[i];
if (tmp < k) k = tmp;
}
printf("%d\n", k);
}
return 0;
}
/**************************************************************
Problem: 1104
User: wangzhengyi
Language: C
Result: Accepted
Time:0 ms
Memory:916 kb
****************************************************************/
约数个数定理对于一个大于1的正整数n可以分解质因数:
n = p1^a1 * p2^a2 * p3^a3 * ... * pn^an
, 则n的正约数的个数为:
(a1 + 1) * (a2 + 1) * ... *(an + 1)
.其中p1,p2,..pn都是n的质因数,a1, a2...an是p1,p2,..pn的指数
证明
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^ak的约数有(ak+1)个
故根据乘法原理:n的约数的个数就是
(a1+1)*(a2+1)*(a3+1)*…* (ak+1)
题目3
题目描述:
输入n个整数,依次输出每个数的约数的个数
输入:
输入的第一行为N,即数组的个数(N<=1000)
接下来的1行包括N个整数,其中每个数的范围为(1<=Num<=1000000000)
当N=0时输入结束。
输出:
可能有多组输入数据,对于每组输入数据,
输出N行,其中每一行对应上面的一个数的约数的个数。
样例输入:
5
1 3 4 6 12
样例输出:
1
2
3
4
6
代码
#include
#include
#define N 40000
typedef long long int lint;
int prime[N], size;
void init()
{
int i, j;
for (prime[0] = prime[1] = 0, i = 2; i < N; i ++) {
prime[i] = 1;
}
size = 0;
for (i = 2; i < N; i ++) {
if (prime[i]) {
size ++;
for (j = 2 * i; j < N; j += i)
prime[j] = 0;
}
}
}
lint numPrime(int n)
{
int i, num, *ansnum, *ansprime;
lint count;
ansnum = (int *)malloc(sizeof(int) * (size + 1));
ansprime = (int *)malloc(sizeof(int) * (size + 1));
for (i = 2, num = 0; i < N && n != 1; i ++) {
if (prime[i] && n % i == 0) {
ansprime[num] = i;
ansnum[num] = 0;
while (n != 1 && n % i == 0) {
ansnum[num] += 1;
n /= i;
}
num ++;
}
}
if (n != 1) {
ansprime[num] = n;
ansnum[num] = 1;
num ++;
}
for (i = 0, count = 1; i < num; i ++) {
count *= (ansnum[i] + 1);
}
free(ansnum);
free(ansprime);
return count;
}
int main(void)
{
int i, n, *arr;
lint count;
init();
while (scanf("%d", &n) != EOF && n != 0) {
arr = (int *)malloc(sizeof(int) * n);
for (i = 0; i < n; i ++) {
scanf("%d", arr + i);
}
for (i = 0; i < n; i ++) {
count = numPrime(arr[i]);
printf("%lld\n", count);
}
free(arr);
}
return 0;
}
/**************************************************************
Problem: 1087
User: wangzhengyi
Language: C
Result: Accepted
Time:190 ms
Memory:1068 kb
****************************************************************/
分解质因数c语言设计思路,深入分析C语言分解质因数的实现方法相关推荐
- c语言设计思路和有点不足,C语言总结报告
1.当初你是如何做出选择计算机专业的决定的? 经过一个学期,你的看法改变了么,为什么? 你觉得计算机是你喜欢的领域吗,它是你擅长的领域吗? 为什么? 当初报考计算机专业,是看到计算机专业在当今社会有良 ...
- c语言设计匀速直线运动,C语言课程设计指导书2015版(试用版)201506.doc
C语言课程设计指导书2015版(试用版)201506精选 C语言课程设计 指导书 信息技术教学与实验中心 2015-04 目 录 第1章 C语言课程设计基本要求-------------- 1 §1. ...
- 单片机c语言设计数据,单片机C语言设计.doc
单片机C语言设计 第5章 单片机C语言程序设计 C-51的特点 C语言作为一种非常方便的语言而得到广泛的支持,很多硬件开发都用C语言编程,如:各种单片机.DSP.ARM等. C语言程序本身不依赖于机器 ...
- C语言设计的管理系统,C语言设计订单管理系统探索.doc
PAGE \* MERGEFORMAT30 C语言课程设计 随米打印订单管理系统 学 院: 计算机与信息科学学院 学生姓名: 谢润发 指导教师: 王新祥 职称 教授 专 业: 网络工程 班 级: 15 ...
- html语言设计表格,HTML标记语言——表格标记
HTML标记语言--表格标记 互联网 发布时间:2008-10-17 18:55:24 作者:佚名 我要评论 点击这里返回网页教学网 HTML教程 栏目. 上文:标记语言--标题 原文出处 ...
- c语言设计题目代码,C语言课程的设计题目.doc
计算机综合训练题目 指导教师: 尹航 训练班级: 完成时间: 2011-2012学年 第1学期 第4周 目 录 01课程设计题目:教师信息管理系统(由 完成)2 02课程设计题目:销售信息统计(由 完 ...
- 用html语言设计如下网页,HTML语言与网页制作教案(5篇可选)
HTML语言与网页制作教案 城步职业中学:刘红兰 教学内容:HTML语言与网页制作 教学目的: 1.掌握输入和运行HTML文件的方法 2.掌握HTML基本的结构标记 3.掌握html在网页特效中的简单 ...
- html语言设计一个网页,HTML语言与网页设计.ppt
HTML语言与网页设计 陈庆章 qzchen@ 2011年05月30日 HTML入门? HTML语言 一种进行网页设计的语言 一种标签式的程序设计语言 一种纯文本式的语言 一种简单易学的语言 一种可以 ...
- 用c语言设计一个菜单界面_用这7个方法设计菜单,让餐厅的利润暴涨
现在很多餐厅看到的菜单上大多是这样设计的:本店主打招牌菜.最佳人气冠军.主厨推荐.必吃美食等. 表面看菜单很全,但是顾客看着菜单却不知道点什么好,另外也会增加顾客的点餐时间,对翻台率有影响...... ...
最新文章
- 从零开始学电脑办公_从零开始学办公软件之办公软件介绍
- MFC模块的动态链接库DLL以及静态链接库LIB编译后的调用
- Oracle中的in 和 not in
- C++匿名对象调用成员函数
- 2020Alibaba数学预选赛第2轮试题
- 计算机二维全息图原理,三维信息加密如何使用计算全息进行
- OPenCV的中cvErode和cvDilate腐蚀和膨胀函数
- 内核线程和用户线程(SMP)
- CodeForces - 589D
- Android 滑动界面实现---Scroller类别 从源代码和开发文档了解(让你的移动布局)...
- Eclipse下配置Tomcat内存
- WM_NOTIFY的使用
- 全国计算机一级成绩分配,计算机一级ms分值分配
- pytorch使用万能模板
- 苹果系统中英文输入法切换_苹果输入法怎么切换_苹果系统如何切换输入法-win7之家...
- PyG教程(2):图数据
- 网络复现之基于TPS的STN网络
- STM32蜂鸣器实例详解
- 登录管理(权限管理) token
- 王峻涛:垂直电子商务将有很大发展
热门文章
- 算法与程序设计(四):贪心算法
- [echarts] 引入png图片或者图标
- lisp倒入excel数据画图_Excel表格数据在cad画出图形-如何将EXCEL中的数据导入到CAD中,绘制成曲线...
- java小车赛跑_Java模拟赛跑过程
- MySQL入门代码(附案例)
- Redis持久化中RDB和AOF有哪些区别?你知道吗!
- 北京Yeslab安全实验室 现任明教教主
- Python学习之调换顺序
- USB摄像头运行ORB-SLAM2
- oracle 10g安装时准备从如下地址启动然后退出了,RedHat AS4 下安装oracle10g