方法一:朴素算法:O(n)。

 #include<bits/stdc++.h>
using namespace std;
int get_num(int n){
int num=;
for(int i=;i<=n;++i)
if(n%i==)num++;
return num;
}
int get_sum(int n){
int tot=;
for(int i=;i<=n;++i)
if(n%i==)tot+=i;
return tot;
}
int main(){
int n;
while(cin>>n){
cout<<get_num(n)<<endl;//求n的因子个数
cout<<get_sum(n)<<endl;//求n的因子数之和
}
return ;
}

方法二:约数个数定理:O(√n)。

 #include<bits/stdc++.h>
using namespace std;
int get_num(int n){
int tot=;
for(int i=;i*i<=n;++i){
if(n%i==){
int x=;
while(n%i==){
n/=i;
x++;
}
tot*=(x+);
}
}
if(n>)tot*=;
return tot;
}
int get_sum(int n){
int tot=;
for(int i=;i*i<=n;i++){
if(n%i==){
int mul=;
while(n%i==){
n/=i;
mul*=i;
}
tot*=(mul*i-)/(i-);
}
}
if(n>)tot*=(n+);
return tot;
}
int main(){
int n;
while(cin>>n){
cout<<get_num(n)<<endl;//求n的因子个数
cout<<get_sum(n)<<endl;//求n的因子数之和
}
return ;
}

实战例题:

题解报告:hdu 2521 反素数

Problem Description

反素数就是满足对于任意i(0<i<x),都有g(i)<g(x),(g(x)是x的因子个数),则x为一个反素数。现在给你一个整数区间[a,b],请你求出该区间的x使g(x)最大。

Input

第一行输入n,接下来n行测试数据
输入包括a,b, 1<=a<=b<=5000,表示闭区间[a,b].

Output

输出为一个整数,为该区间因子最多的数.如果满足条件有多个,则输出其中最小的数.

Sample Input

3
2 3
1 10
47 359

Sample Output

2 6 240

Hint

2的因子为:1 2

10的因子为:1 2 5 10
解题思路:求区间[a,b]中最小的元素x,使得所含的因子个数最大,由于数据比较小,暴力枚举区间[a,b]并套一下公式即可。时间复杂度为O($n\sqrt{n}$)。
AC代码(15ms):
 #include<bits/stdc++.h>
using namespace std;
int n,a,b,maxnum,maxindex;
int get_num(int n){//求n的因子个数
int tot=;
for(int i=;i*i<=n;++i){
if(n%i==){
int x=;
while(n%i==){
n/=i;
x++;
}
tot*=(x+);
}
}
if(n>)tot*=;
return tot;
}
int main(){
while(cin>>n){
while(n--){
cin>>a>>b;maxnum=;maxindex=a;
for(int i=a;i<=b;++i){
int tmp=get_num(i);
if(tmp>maxnum){maxindex=i;maxnum=tmp;}
}
cout<<maxindex<<endl;
}
}
return ;
}

题解报告:hdu 1215 七夕节

Problem Description

七夕节那天,月老来到数字王国,他在城门上贴了一张告示,并且和数字王国的人们说:"你们想知道你们的另一半是谁吗?那就按照告示上的方法去找吧!"
人们纷纷来到告示前,都想知道谁才是自己的另一半.告示如下:
数字N的因子就是所有比N小又能被N整除的所有正整数,如12的因子有1,2,3,4,6.
你想知道你的另一半吗?

Input

输入数据的第一行是一个数字T(1<=T<=500000),它表明测试数据的组数.然后是T组测试数据,每组测试数据只有一个数字N(1<=N<=500000).

Output

对于每组测试数据,请输出一个代表输入数据N的另一半的编号.

Sample Input

3
2
10
20

Sample Output

1
8
22
解题思路:求n的因子数之和,简单套一下公式,时间复杂度为O($t\sqrt{n}$),最坏大概为3e8。
AC代码一(93ms):
 #include<bits/stdc++.h>
using namespace std;
int t,n;
int get_sum(int n){
int tot=;
for(int i=;i*i<=n;++i){
if(n%i==){
int mul=;
while(n%i==){
n/=i;
mul*=i;
}
tot*=(mul*i-)/(i-);
}
}
if(n>)tot*=(n+);
return tot;
}
int main(){
while(~scanf("%d",&t)){
while(t--){
scanf("%d",&n);
printf("%d\n",get_sum(n)-n);
}
}
return ;
}

AC代码二(109ms):简单打表,时间复杂度为nlogn。

 #include<bits/stdc++.h>
const int maxn=;
using namespace std;
int t,n,sum[maxn];
int main(){
for(int i=;i<maxn/;++i)//因子i
for(int j=i*;j<maxn;j+=i)//j是i的倍数,即j的因子是i
sum[j]+=i;//sum[j]表示其所有因子之和
while(~scanf("%d",&t)){
while(t--){
scanf("%d",&n);
printf("%d\n",sum[n]);
}
}
return ;
}

题解报告:hdu 1999 不可摸数

Problem Description

s(n)是正整数n的真因子之和,即小于n且整除n的因子和.例如s(12)=1+2+3+4+6=16.如果任何数m,s(m)都不等于n,则称n为不可摸数.

Input

包含多组数据,首先输入T,表示有T组数据.每组数据1行给出n(2<=n<=1000)是整数。

Output

如果n是不可摸数,输出yes,否则输出no

Sample Input

3
2
5
8

Sample Output

yes
yes
no
解题思路:题意就是给一个n(n表示为s(m)的真因子数之和),如果任意数m的真因子之和在n∈[2,1000],那么s(m)就不是不可摸数,否则为不可摸数。打表可以发现数m只需枚举到10^6便可将真因数之和为1000以内的全部找出来。于是这道题也可以直接打印出89个不可摸数,简单标记一下即可。
AC代码一(62ms):
 #include<bits/stdc++.h>
const int maxn=;
using namespace std;
int t,n,sum[maxn];bool f[maxn];
int main(){
memset(f,false,sizeof(f));
for(int i=;i<maxn/;++i)//因子i
for(int j=i*;j<maxn;j+=i)//j是i的倍数,即j的因子是i
sum[j]+=i;//sum[j]表示其所有因子之和
for(int i=;i<maxn;++i)
if(sum[i]<=)f[sum[i]]=true;//如果真因数之和在1000以内,那么sum[i]不是不可摸数
while(~scanf("%d",&t)){
while(t--){
scanf("%d",&n);
if(f[n])printf("no\n");
else printf("yes\n");//n为不可摸数
}
}
return ;
}

AC代码二(0ms):

 #include<bits/stdc++.h>
using namespace std;
int t,n;bool f[];
int main(){
int obj[]={,,,,,,,,,,
,,,,,,,,,,,,
,,,,,,,,,,,,
,,,,,,,,,,,,
,,,,,,,,,,,,
,,,,,,,,,,,,
,,,,,,,,,,,,
,,,,,,};
memset(f,false,sizeof(f));
for(int i=;i<;++i)f[obj[i]]=true;//1000以内共有89个数为不可摸数
while(~scanf("%d",&t)){
while(t--){
scanf("%d",&n);
if(f[n])printf("yes\n");
else printf("no\n");
}
}
return ;
}

题解报告:hdu 1299 Diophantus of Alexandria

Problem Description

Diophantus of Alexandria was an egypt mathematician living in Alexandria. He was one of the first mathematicians to study equations where variables were restricted to integral values. In honor of him, these equations are commonly called diophantine equations. One of the most famous diophantine equation is x^n + y^n = z^n. Fermat suggested that for n > 2, there are no solutions with positive integral values for x, y and z. A proof of this theorem (called Fermat's last theorem) was found only recently by Andrew Wiles.
Consider the following diophantine equation: 
1 / x + 1 / y = 1 / n where x, y, n ∈ N+ (1)Diophantus is interested in the following question: for a given n, how many distinct solutions (i. e., solutions satisfying x ≤ y) does equation (1) have? For example, for n = 4, there are exactly three distinct solutions: 
1 / 5 + 1 / 20 = 1 / 4
1 / 6 + 1 / 12 = 1 / 4
1 / 8 + 1 / 8 = 1 / 4Clearly, enumerating these solutions can become tedious for bigger values of n. Can you help Diophantus compute the number of distinct solutions for big values of n quickly?

Input

The first line contains the number of scenarios. Each scenario consists of one line containing a single number n (1 ≤ n ≤ 10^9). 

Output

The output for every scenario begins with a line containing "Scenario #i:", where i is the number of the scenario starting at 1. Next, print a single line with the number of distinct solutions of equation (1) for the given value of n. Terminate each scenario with a blank line. 

Sample Input

2
4
1260

Sample Output

Scenario #1:
3
Scenario #2:
113
解题思路:首先,根据分式等式的性质可得x、y均大于n,假设y=n+k(k≥1),将其代入1/x+1/y=1/n化简得x=n^2/k+n,因为x∈N+,所以k|n2(k≥1),即求n2的因子个数sum,则最终的答案为:(sum+1)/2,因为必定有一半的x、y这一对是重复了的,只是值互换了而已,举个例子:当n=4时,
x=     8  6   5
k=  1   2  4  8  16(n^2=16的所有因子)
y=     6  8    
显然此时有3组解,也就是满足n2/k == k,(其中k|n2)的解有(sum+1)/2组。
AC代码:
 #include<bits/stdc++.h>
using namespace std;
typedef long long LL;int t;LL n;
LL get_num(LL x){
LL ans=,tp;
for(LL i=;i*i<=x;++i){
if(x%i==){
tp=;
while(x%i==)tp++,x/=i;
ans*=(+*tp);
}
}
if(x>)ans*=;
return ans;
}
int main(){
while(cin>>t){
for(int i=;i<=t;++i){
cin>>n;
cout<<"Scenario #"<<i<<":\n"<<(get_num(n)+)/<<"\n"<<endl;
}
}
return ;
}

题解报告:NYOJ #66 分数拆分

描述

现在输入一个正整数k,找到所有的正整数x>=y,使得1/k=1/x+1/y.

输入

第一行输入一个整数n,代表有n组测试数据。
接下来n行每行输入一个正整数k

输出

按顺序输出对应每行的k找到所有满足条件1/k=1/x+1/y的组合。

样例输入

2
2
12

样例输出

1/2=1/6+1/3
1/2=1/4+1/4
1/12=1/156+1/13
1/12=1/84+1/14
1/12=1/60+1/15
1/12=1/48+1/16
1/12=1/36+1/18
1/12=1/30+1/20
1/12=1/28+1/21
1/12=1/24+1/24

解题思路:由x>=y且x、y均大于k可知1/x<=1/y,1/k-1/y<=1/y,即k<y<=2k,所以只需控制y的枚举范围即可。

AC代码:

 #include<bits/stdc++.h>
using namespace std;
int t,k;
int main(){
while(cin>>t){
while(t--){
cin>>k;
for(int i=k+;i<=*k;++i)//枚举k+1~2k
if(i*k%(i-k)==)printf("1/%d=1/%d+1/%d\n",k,i*k/(i-k),i);//通分,其中(i-k)|(i*k)。
}
}
return ;
}

财大情侣

Time Limit: 2000/1000ms (Java/Others)

Problem Description:

高考完终于来到梦寐以求的财经大学,男三女七有木有!食堂吃饭前后左右三排都是女生有木有!该学校总共有n个人,每个人都有对应自己的魅力值,从1到n。现规定每个人的情商为每个人魅力值的约数(不包括本身)之和。在广财有个很神奇的事情,当一个人的情商刚好等于另一个的魅力值的时候,这两个人就meant to be情侣(不是一男一女怎么办?随便吧)。
如:220:1+2+4+5+10+11+20+22+44+55+110=284
 284:1+2+4+71+142=220
这样魅力值为220和284的就成为一对啦。
给定两个数,a,b,找出所有a和b之间(inclusive)成对的。

Input:

输入包含多组测试数据,每组数据输入两个数a,b (0<=a,b<=100000).

Output:

对于每组测试,输出a和b之间的所有情侣。每对情侣占一行,魅力值小的在前面。对于多对情侣,按情侣中魅力值小的排序。参考sample output。

Sample Input:

200 1300
1 200

Sample Output:

220 284
1184 1210
NO
解题思路:求每个数的所有真因子之和,简单打个表再判断一下即可。
AC代码:
 #include<bits/stdc++.h>
const int maxn=;//范围大一点,否则会出现越界的情况
using namespace std;
int a,b,x,sum[maxn];bool flag;
int main(){
for(int i=;i<maxn/;++i)//求因子数之和
for(int j=i*;j<maxn;j+=i)//j初始为i的2倍,以后以i的步长增长,累加j的真因子
sum[j]+=i;
while(cin>>a>>b){
if(a>b)swap(a,b);flag=false;
for(int i=a;i<=b;++i)//先判断是不是一对情侣,并且前面的因子数之和x要在区间范围内,且前一个数i要小于后一个数x
if(i==sum[x=sum[i]]&&x<=b&&i<x){flag=true;cout<<i<<' '<<x<<endl;}
if(!flag)cout<<"NO"<<endl;
}
return ;
}

求n的因子个数与其因子数之和的更多相关文章

  1. Soldier and Number Game---cf546D(打表求n的素因子个数)

    题目链接:http://codeforces.com/problemset/problem/546/D 题意: 给出一个n,n开始是a!/b!,每次用一个x去整除n得到新的n,最后当n变成1的时候经过 ...

  2. Almost All Divisors(求因子个数及思维)

    ---恢复内容开始--- We guessed some integer number xx. You are given a list of almost all its divisors. Alm ...

  3. LightOj1028 - Trailing Zeroes (I)---求因子个数

    题目链接:http://lightoj.com/volume_showproblem.php?problem=1028 题意:给你一个数 n (1<=n<=10^12), 然后我们可以把它 ...

  4. POJ 2992 Divisors (求因子个数)

    题意:给n和k,求组合C(n,k)的因子个数. 这道题,若一开始先预处理出C[i][j]的大小,再按普通方法枚举2~sqrt(C[i][j])来求解对应的因子个数,会TLE.所以得用别的方法. 在说方 ...

  5. Acdream1084 寒假安排 求n!中v因子个数

    题目链接:pid=1084">点击打开链接 寒假安排 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 128000/64000 ...

  6. Trailing Zeroes (I) LightOJ - 1028(求因子个数)

    题意: 给出一个N 求N有多少个别的进制的数有后导零 解析: 对于一个别的进制的数要转化为10进制 (我们暂且只分析二进制就好啦) An * 2^(n-1) + An-1 * 2^(n-2) + `` ...

  7. 第13届景驰-埃森哲杯广东工业大学ACM程序设计大赛-等式(求$N^2$的因子个数)

    一.题目链接 https://www.nowcoder.com/acm/contest/90/F 二.题面 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 32768K,其他语言655 ...

  8. 求n!中因子k的个数

    思路: 求n的阶乘某个因子k的个数,如果n比较小,可以直接算出来,但是如果n很大,此时n!超出了数据的表示范围,这种直接求的方法肯定行不通.其实n!可以表示成统一的方式. n!=(km)*(m!)*a ...

  9. BZOJ3994:约数个数和(莫比乌斯反演:求[1,N]*[1,M]的矩阵的因子个数)

    Description  设d(x)为x的约数个数,给定N.M,求   Input 输入文件包含多组测试数据. 第一行,一个整数T,表示测试数据的组数. 接下来的T行,每行两个整数N.M. Outpu ...

随机推荐

  1. thinkphp Class 'PDO' not found 错误

    thinkphp Class 'PDO' not found 错误,原因mysql5.7.26缺少pdo驱动,需要安装php的pdo和pdo_mysql扩展 本文以centOS为例 1.进入PHP源码 ...

  2. 【POJ 1159】Palindrome

    [POJ 1159]Palindrome 近期各种题各种奇葩思路已经司空见惯了...又新出个滚动数组= = 该题另一点须要知道 最少须要补充的字母数 = 原序列S的长度 - S和S'的最长公共子串长度 ...

  3. IE8 下背景图标不显示

    如图所示 : 微博微信前方各应有个图标,但是IE8下图标没有显示 css如下 .weibo_icon{background: url(../ieImages/weibo_icon.png)no-rep ...

  4. iOS 如何改变表视图分割线在iOS7中的默认偏移

    - (void)viewDidLoad { [super viewDidLoad]; self.automaticallyAdjustsScrollViewInsets = NO; if ([self ...

  5. ios很好的开源库

    Tim9Liu9/TimLiu-iOS 自己总结的iOS.mac开源项目及库,持续更新.. 目录 UI 下拉刷新 模糊效果 AutoLayout 富文本 图表 表相关与Tabbar 隐藏与显示 HUD ...

  6. POJ2752 Seek the Name, Seek the Fame —— KMP next数组

    题目链接:https://vjudge.net/problem/POJ-2752 Seek the Name, Seek the Fame Time Limit: 2000MS   Memory Li ...

  7. css中链接的4个状态 link、visited 、hover、active 顺序与设置

    link       未点时链接色visited  已点过的链接色hover   光标放入链接上的色active   点击时的样式 顺序一定要按link.visited.hover.active来设置 ...

  8. MYSQL初级学习笔记九:MySQL索引的使用!(视频序号:初级_51)

    知识点十一:索引的使用(51) 什么是索引: 索引的定义: 在关系型数据库中,索引是一种与表有关的数据库结构,它可以使对应于表的SQL语句执行的更快.索引的作用相当于图书的目录,可以 根据目录中的页码 ...

  9. lucene 5的测试程序——API变动太大

    package hello; import java.io.BufferedReader; import java.io.File; import java.io.FileReader; import ...

  10. jetty源码剖析

    最近使用jetty自己写了一个web server,现在闲了花了一天的时间看了一jetty的源代码,主要以server的启动为主线,进行了剖析,经过阅读对jetty的源码大赞,写的简洁.清晰,架构也不 ...