HDU 2601An easy problem-素数的运用,暴力求解
Time Limit: 3000MS | Memory Limit: 32768KB | 64bit IO Format: %I64d & %I64u |
Description
One day Teddy met a old man in his dream , in that dream the man whose name was“RuLai” gave Teddy a problem :
Given an N , can you calculate how many ways to write N as i * j + i + j (0 < i <= j) ?
Teddy found the answer when N was less than 10…but if N get bigger , he found it was too difficult for him to solve.
Well , you clever ACMers ,could you help little Teddy to solve this problem and let him have a good dream ?
Input
Output
Sample Input
2
1
3
Sample Output
0
1
能够表示为N=i*j+i+j
所以能够化作N+1=(i+1)*(j+1);
如此就有两种思路去做,一暴力枚举,从这里能够看出来i<=sqrt(N+1);
所以一个循环就能够解决,第一份代码就是。可是耗时才一点就过题目提供的3s了。怎么办,是否还有
更好的解决方式呢,有的,N+1=(i+1)*(j+1)能够知道N+1是i+1以及j+1的倍数
如此就能够转换成求解约数的个数(约数是什么,请读者自己百度了解)
N+1=a1^p1*a2^p2*a3^p3....an^pn
当中ai是代表着质数,这个的意思是不论什么大于1的数都能够转换为有限个质数因子的乘积
如此。能够用排列组合来求。第一种有p1+1选择(能够选择0...p1)另外一种有p2+1选择(能够选择0...p2)....第n种有pn+1选择(能够选择0...pn)
所以约数个数ans=(p1+1)*(p2+1)*(p3+1)*(p4+1)*....*(pn+1)(里面还有减去1和n由于他们不属于题目要求范围)
当N+1是全然平方数的话,那么除了1以及N+1本身外,唯有最中间的约数是仅仅计算了一次,其它的数都反复的计算了两次,
当N+1不是全然平方数的话,那么除了1以及N+1本身外,其它的数都反复的计算了两次
所以能够分开推断输出,也能够直接转换输出就像代码中一样,(ans+1)/2-1,当为全然平方数时,我们须要加一除二才干使正确的结果
至于减去一,就是前面的去掉1和n这两个不符合条件的数
当为不全然平方数,ans/2-1就能够了,可是为了合成一个式子,(ans+1)/2-1,是能够取代ans/2-1的
为什么,由于(4+1)/2==4/2,这是不会影响终于结果的。
/*
Author: 2486
Memory: 1416 KB Time: 2823 MS
Language: G++ Result: Accepted
*/
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long LL;
const int maxn=1e5;
int t;
LL n;
int main() {
scanf("%d",&t);
while(t--) {
scanf("%I64d",&n);
if(n==0||n==1) {
printf("0\n");
continue;
}
int cnt=0;
for(int i=1; i<=sqrt(n); i++) {
if((n+1)%(i+1)==0&&(n+1)/(i+1)>=i+1)cnt++;
}
printf("%d\n",cnt);
}
}
/*
Author: 2486
Memory: 1592 KB Time: 46 MS
Language: G++ Result: Accepted
*/
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
const int maxn=100000+5;
LL prime[maxn];
bool vis[maxn];
int T,cnt;
LL N;
void primes() { //初始化素数列表
cnt=0;
for(int i=2; i<maxn; i++) {
if(vis[i])continue;
prime[cnt++]=i;
for(int j=i*2; j<maxn; j+=i) {
vis[j]=true;
}
}
}
void solve(LL n) {
LL ans=1;
for(int i=0; prime[i]*prime[i]<=n; i++) {
if(n%prime[i]==0) {
int s=0;
while(n%prime[i]==0)n/=prime[i],s++;
ans*=(s+1);
}
if(n==1)break;
}
if(n>1)ans*=2;
printf("%I64d\n",(ans+1)/2-1);
}
int main() {
primes();
scanf("%d",&T);
while(T--) {
scanf("%I64d",&N);
N++;
solve(N);
}
return 0;
}
HDU 2601An easy problem-素数的运用,暴力求解的更多相关文章
- HDU 5475An easy problem 离线set/线段树
An easy problem Time Limit: 8000/5000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) ...
- hdu 5475 An easy problem(暴力 || 线段树区间单点更新)
http://acm.hdu.edu.cn/showproblem.php?pid=5475 An easy problem Time Limit: 8000/5000 MS (Java/Others ...
- HDU 5475 An easy problem 线段树
An easy problem Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/showproblem.php?pi ...
- 数据结构(主席树):HDU 4729 An Easy Problem for Elfness
An Easy Problem for Elfness Time Limit: 5000/2500 MS (Java/Others) Memory Limit: 65535/65535 K (J ...
- HDU 2132 An easy problem
http://acm.hdu.edu.cn/showproblem.php?pid=2132 Problem Description We once did a lot of recursional ...
- HDOJ(HDU) 2132 An easy problem
Problem Description We once did a lot of recursional problem . I think some of them is easy for you ...
- hdu 2055 An easy problem (java)
问题: 開始尝试用字符做数组元素.可是并没实用. 在推断语句时把a z排出了. An easy problem Time Limit: 1000/1000 MS (Java/Others) Me ...
- HDU 5475:An easy problem 这题也能用线段树做???
An easy problem Time Limit: 8000/5000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) ...
- HDU 2098 分拆素数和
HDU 2098 分拆素数和 Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768K (Java/Others) [题目描述 ...
随机推荐
- Python自动发送邮件(可带附件)
Python内置对SMTP的支持,可以发送纯文本邮件.HTML邮件以及带附件的邮件. Python对SMTP支持有smtplib和email两个模块,email负责构造邮件,smtplib负责发送邮件 ...
- hlgoj 1766 Cubing
模拟.下图是我做的小模型. #include <iostream> #include <stdio.h> #include <queue> #include < ...
- 【软考5】解释型 or 编译型
导读:在上篇博客中,我们了解到,目前的编码语言经过不断的发展,已经经历了机器语言--汇编语言--高级语言的过程.虽然我们的编码语言在不停的升级,但作为计算机来说,它始终是一个只能理解0和1构成的机器语 ...
- 【转】Eric's并发用户数估算与Little定律的等价性
转自:http://www.cnblogs.com/hundredsofyears/p/3360305.html 在国内性能测试的领域有一篇几乎被奉为大牛之作的经典文章,一个名叫Eric Man Wo ...
- 使用Unity做2.5D游戏教程(一)
最近在研究Unity 3D,看了老外Marin Todorov写的教程很详细,就翻译过来以便自己参考,翻译不好的地方请多包涵. 如果你不了解2.5D游戏是什么,它基本上是个3D游戏而你可以想象是压扁的 ...
- 【bzoj3489】A simple rmq problem 三维KD-tree
题目描述 因为是OJ上的题,就简单点好了.给出一个长度为n的序列,给出M个询问:在[l,r]之间找到一个在这个区间里只出现过一次的数,并且要求找的这个数尽可能大.如果找不到这样的数,则直接输出0.我会 ...
- POJ——3159Candies(差分约束SPFA+前向星+各种优化)
Candies Time Limit: 1500MS Memory Limit: 131072K Total Submissions: 28071 Accepted: 7751 Descrip ...
- LightOJ——1066Gathering Food(BFS)
1066 - Gathering Food PDF (English) Statistics Forum Time Limit: 2 second(s) Memory Limit: 32 MB W ...
- jQuary的相关动画效果
第一种:该方法隐藏所有 <p> 元素: <html> <head> <script type="text/javascript" src= ...
- msp430项目编程52
msp430综合项目---扩展项目二52 1.电路工作原理 2.代码(显示部分) 3.代码(功能实现) 4.项目总结