1173: 这是物理学的奇迹!!

题目描述

goagain在做物理电学实验时需要一个2Ω的电阻,但是他发现他的实验台上只剩下了3Ω,4Ω,5Ω,6Ω的电阻若干,于是goagain把两个4Ω的电阻并联起来,发现效果也不错,于是他惊呼,这是物理学的奇迹!!

实验结束之后goagain思考了一下,其实用1个3Ω的电阻和一个6Ω的电阻并联一下也是不错的,但是似乎没有其他的方法可以并联两个电阻得到等效的2Ω电阻了.

goagain就想啊,如果我需要一个k Ω的电阻,但是我的手头上只有k+1 Ω,k+2 Ω,.....,一直到正无穷的整数Ω电阻,那么我有多少种办法把他们其中的两个并联起来,得到一个k Ω的电阻呢?

goagain仔细想了想,终于发现自己的智力不足以解决这样的问题,于是他不得不感叹:这真是物理学的奇迹!

输入要求

有多组测试数据,每组数据包含一个正整数k(k<2^16),处理到文件结束

输出要求

对于每组测试数据,输出一个正整数,代表方法的总数

假如输入

123

应当输出

122

提示

1/2=1/3+1/6=1/4+1/4

1/3=1/4+1/12=1/6+1/6


由公式 1/k=1/a+1/b ,转换成 1/k=1/(c+k)+1/(d+k)。进行化简。得到这样的一条式子:  k*k = c*d;  另N=k*k。  则N=c*d。很显然的,题目转换成了求能整除N的数的个数。

由因式分解可知。

设pi为N的因子。有 N=p1^a * p2^b * p3^c ..... 等。  可以知道。p1可以取0到a个,p2可以取0到b的个,p3可以取0到c个。因此可以得出这样的式子。  能整除N的数的个数=(a+1)*(b+1)*(c+1)*....

因此。我们只需快速的对N进行因式分解,即可快速的算出答案,不过还得留意N是n*n的情况,需要多加一次。

最后记得将答案除以2.因为ab重复了。  当n=1时。是个特殊情况

筛法求素数的模板getprm()快速获得素数

  

求 k*k分解时有个小小的优化。。就是可以求出每个k的a,b,c最后算k*k时,是2a,2b,2c。这个优化还是比较明显的。从1100MS掉到了220MS

当然因为k < 2^16次。我们可以直接设计个方法用数组记录答案。理论上存在的。不过设计起来比较麻烦。

/*
* @user ipqhjjybj
* @Time
* @data 20130701
*/
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <ctime> #include <iostream>
#include <cmath>
#include <algorithm>
#include <numeric>
#include <utility> #include <cstring>
#include <vector>
#include <stack>
#include <queue>
#include <map>
#include <string>
using namespace std; #define inf 0x3f3f3f3f
#define MAXN 800000 #define clr(x,k) memset((x),(k),sizeof(x))
#define clrn(x,k) memset((x),(k),(n+1)*sizeof(int))
#define cpy(x,k) memcpy((x),(k),sizeof(x))
#define Base 10000 typedef vector<int> vi; #define foreach(it,c) for(vi::iterator it = (c).begin();it != (c).end();++it) #define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
bool is[MAXN];int prm[MAXN];
int k;
int getprm(int n){
int i,j,k;
for(i=0;i<n;i++) is[i]=true;
is[0]=is[1]=false;
k=0;
for(i=2;i<n;i++){
if(is[i]){
prm[k++]=i;
for(j=i+i;j<n;j+=i)
is[j]=false;
}
}
return k;
} int fenjie(int n){
int i,zn,j,ans=1;
for(i = 0,zn=prm[i];prm[i]<=n;++i,zn = prm[i]){
if(n%zn==0){
j=1;
n/=zn;
while(!(n%zn)){
n/=zn;j++;
}
ans*=(j+j+1); //我们分解的可是n*n哦。。所以是2*j
} }
return ans+1; //当n*n情况时!
}
int main(){
long long n;
k=getprm(MAXN); while(scanf("%lld",&n)!=EOF){
if(n==1)printf("1\n");
else printf("%d\n",fenjie(n)>>1);
}
return 0;
}

NEU 1173: 这是物理学的奇迹!! 分解质数的更多相关文章

  1. Acwing-197-阶乘分解(质数)

    链接: https://www.acwing.com/problem/content/199/ 题意: 给定整数 N ,试把阶乘 N! 分解质因数,按照算术基本定理的形式输出分解结果中的 pi 和 c ...

  2. javascrip 求最大公因数(分解质数法)发生的问题

    //这是求一个数的质因数,例如:12=2*2*3 其中2,3都是质数.function primeArray(n, array) { array = new Array(); for (var i = ...

  3. 快速质因数分解及素性测试&ABC142D

    首先,这个整数的标准分解非常的显然易见对吧: 一般我们要把一个数分解成这个样子我们可以这样写: #include<cstdio> ],w[],k; void factorize(int n ...

  4. 蒟蒻修养之cf橙名计划

    因为太弱,蒟蒻我从来没有上过div1(这就是今年的最后愿望啊啊啊啊啊)已达成................打cf几乎每次都是fst...........所以我的cf成绩图出现了惊人了正弦函数图像.. ...

  5. NOIP2017 国庆郑州集训知识梳理汇总

    第一天 基础算法&&数学 day1难度测试 如果要用一个词来形容上午的测试,那真是体无完肤.  成绩: 题目 成绩 评价 T1 50 一般 T2 10 大失所望 T3 0 差 基础算法 ...

  6. CSUST 四月选拔赛个人题解

    这场比赛演的逼真,感谢队友不杀之恩 总结:卡题了赶紧换,手上捏着的题尽快上机解决 http://csustacm.com:4803/ 1113~1122 1113:六学家 题意:找出满足ai+aj=a ...

  7. 洛谷P2480 [SDOI2010]古代猪文(卢卡斯定理+中国剩余定理)

    传送门 好吧我数学差的好像不是一点半点…… 题目求的是$G^{\sum_{d|n}C^d_n}mod\ 999911659$ 我们可以利用费马小定理$a^{k}\equiv a^{k\ mod\ (p ...

  8. 博客索引and题目列表

    目录 笔记整理 计划 要学的东西 缺省源 要做的题 搜索 高斯消元 矩阵 排列组合 2019.7.9 2019.7.10 kmp ac自动机 2019.7.11 2019.7.15 笔记整理 1.同余 ...

  9. F - Minimum Sum LCM

    LCM (Least Common Multiple) of a set of integers is defined as the minimum number, which is a multip ...

随机推荐

  1. Qt核心剖析: 寻找 QObject 的源代码

    本来打算把<Qt学习之路>作为一个类似教程的东西,所以就不打算把一些关系到源代码的内容放在那个系列之中啦.因此今天就先来看一个新的开始吧!这个系列估计不会进展很快,因为最近公司里面要做 f ...

  2. Qt核心剖析:信息隐藏(三篇)

    http://devbean.blog.51cto.com/448512/335550 http://devbean.blog.51cto.com/448512/325581 http://devbe ...

  3. 娓娓道来c指针 (0)c语言的梦魇:c指针

    (0)c语言的梦魇:c指针 序 c语言中有一个重点:c指针.它也是一个难点.当然,这是一句废话:重点往往也是难点.在c标准中,对指针的定义是这种: 指针的类型是derived from其他类型,也就是 ...

  4. [置顶] cocos2d-x 3.0游戏开发xcode5帅印博客教学 003.[HoldTail]游戏世界以及背景画面

    cocos2d-x 3.0游戏开发xcode5帅印博客教学 003.[HoldTail]游戏世界以及背景画面 写给大家的前言,在学习cocos2d-x的时候自己走了很多的弯路,也遇到了很多很多问题,不 ...

  5. delphi高手突破之异常及错误处理

    什么是异常?为什么要用它? 所谓“异常”是指一个异常类的对象.Delphi的VCL中,所有异常类都派生于Exception类.该类声明了异常的一般行为.性质.最重要的是,它有一个Message属性可以 ...

  6. POJ 2404 Jogging Trails

    Jogging Trails Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 2122   Accepted: 849 Des ...

  7. jquery mobile自己定义webapp开发实例(一)——前言篇

    用jquery mobile做了一段时间的webapp开发,准备用自己的一个小demo做一个模块化的分享 点击demo演示 手机演示二维码: 此demo已经是比較老的版本号,用户体验流畅度确实还存在非 ...

  8. element.style覆盖了我的样式!!

    原文:element.style覆盖了我的样式!! 有时候在写css时,显示效果会出现非常诡异的效果 不知道有没有遇到这种 css: #logo{ border: solid 1px blue; } ...

  9. javascript 中对this关键字的一些理解

    var self = this 可能会觉得奇怪为什么我会定义一个_self的变量, 因为在js里,this不用对于其他的对象语言,他的解析过程与运行过程中this会改变的.这里简单说说js里this的 ...

  10. HDU 4070 Phage War

    贪心,t 大的放到前面...因为感染所有cell需要的phage的总数是一定的,所以产生phage需要的时间是一定的,只需要考虑用来感染的时间,这样考虑的话,把 t 小的放后面的话,可以发现总时间的最 ...