唯一分解定理


Uva10791

题意:

输入整数n,要求至少两个正整数,使得他们的最小公倍数为n,且这些整数的和最小

解法:

首先假设我们知道了一系列数字a1,a2,a3……an,他们的LCM是n,那么什么时候他们是最优解呢,当他们两两互质的时候

为了方便我们以两个数来说明问题。

a和b的LCM是n,GCD是m,那么n=a/m*b , 它们的和就是sum=a+b;

如果m不为1(即a和b不互质),那么我们为什么不优化一下,将a变为a=a/m呢?,改变后a和b的LCM依然是n,但是他们的和显然减少了

所以我们得到最重要的一个性质,要想a1,a2,a3……an的和最小,要保证他们两两互质,只要存在不互质的两个数,就一定可以近一步优化

那我们怎么保证两两互质呢?方法其实很简单,直接分解质因子

例如24=2*2*2*3 , 只能分解为8和3,因为这里有3个2,这3个2必须在一起,如果分开了这3个2,这出现有两个数会有一个公共的质因子2,并且会使这两个数的LCM不是24

再例如72=2*2*2*3*3,只能分为8和9,因为3个2和2个3都不能分开,他们必须在一次

所以,我们将一个数n分解为质因子后,顺便做一个处理,在除干净一个质因子的同时,将他们乘起来作为一个因子,处理完后会得到多个因子,他们之间同样满足两两互质的性质

然后是进一步的分析

例如264600=8*27*25*49  , 只是由3个2,3个3,2个5,2个7,处理后得到的因子,那么8,27,25,49的LCM是264600,并且两两互质,他们还要不要处理呢?不需要了,直接将他们加起来就是我们要的答案!为什么呢?可以将8,27,25,49这些数字乘起来,无论怎样乘都好,最后得到的数字它们的LCM依然是n,但是乘起来再相加显然比直接相加要大得多!

所以我们已经得到了这个问题的解法

1.将一个数分解成质因子,将相同的因子乘起来作为一个处理后的因子

2.将处理后得到的多个因子直接相加就是答案

3.因为题目说只要需要两个数字,所以对于1和素数我们需要小心。对于素数,我们只能分解出一个因子就它自己,对于1一个因子都分解不出来(我们不把1当做因子),他们的答案都是n+1,因为只有1和n的LCM是n

 /*
唯一分解定理的应用,work_quality_factor就是分解质因数的板子
将一个数分解质因数,将他们所有相同的因子乘起来作为一个新的因子,最后的和就是这些因子和
*/
#include<cstdio>
#include<cmath>
using namespace std;
typedef long long ll;
const int maxn = ;
ll fac[maxn], frq[maxn]; ll poww(ll a, ll b) {
ll ans = , base = a;
while (b != ) {
if (b & != )
ans *= base;
base *= base;
b >>= ;
}
return ans;
} ll work_quality_factor(ll n, ll quality_fac[], ll frequency[])
{//n是待分解的数,quality_fac[]会存放它包含的质因子,而frequency[]存放对应次数
//如q_f[k]=7,fre[k]=2就表示质因数分解后里面包含有7,且次数是2
//函数返回有几种质因子,比如分解了25就返回1,分解28返回2
ll res, temp, i;
res = ;
temp = n;
for (i = ; i*i <= temp; i++)
if (temp%i == )
{
quality_fac[res] = i;
frequency[res] = ;
while (temp%i == )
{
temp = temp / i;
frequency[res]++;
}
res++;
}
if (temp > )
{
quality_fac[res] = temp;
frequency[res++] = ;
}
return res;
} int main() {
ll n; int kase = ;
while (scanf("%lld", &n) && n) {
ll num = work_quality_factor(n, fac, frq);
ll ans = ;
if (num == || num == ) ans = n + ;
else {
for (int i = ; i < num; i++) {
ans += poww(fac[i], frq[i]);
}
}
printf("Case %d: %lld\n", kase++, ans);
}
return ;
}

Uva10791 唯一分解定理模板的更多相关文章

  1. B - Common Divisors (codeforces)数论算法基本定理,唯一分解定理模板

    You are given an array aa consisting of nn integers. Your task is to say the number of such positive ...

  2. UVA10791-Minimum Sum LCM(唯一分解定理基本应用)

    原题:https://vjudge.net/problem/UVA-10791 基本思路:1.借助唯一分解定理分解数据.2.求和输出 知识点:1.筛法得素数 2.唯一分解定理模板代码 3.数论分析-唯 ...

  3. FZU 1075 分解素因子【数论/唯一分解定理/分解素因子裸模板】

    [唯一分解定理]:https://www.cnblogs.com/mjtcn/p/6743624.html 假设x是一个正整数,它的值不超过65535(即1<x<=65535),请编写一个 ...

  4. NOIP2009Hankson 的趣味题[唯一分解定理|暴力]

    题目描述 Hanks 博士是 BT (Bio-Tech,生物技术) 领域的知名专家,他的儿子名叫 Hankson.现 在,刚刚放学回家的 Hankson 正在思考一个有趣的问题. 今天在课堂上,老师讲 ...

  5. UVA - 10375 Choose and divide[唯一分解定理]

    UVA - 10375 Choose and divide Choose and divide Time Limit: 1000MS   Memory Limit: 65536K Total Subm ...

  6. uva10375 Choose and Divide(唯一分解定理)

    uva10375 Choose and Divide(唯一分解定理) 题意: 已知C(m,n)=m! / (n!*(m-n!)),输入整数p,q,r,s(p>=q,r>=s,p,q,r,s ...

  7. 1341 - Aladdin and the Flying Carpet ---light oj (唯一分解定理+素数筛选)

    http://lightoj.com/volume_showproblem.php?problem=1341 题目大意: 给你矩形的面积(矩形的边长都是正整数),让你求最小的边大于等于b的矩形的个数. ...

  8. UVA 10375 Choose and divide【唯一分解定理】

    题意:求C(p,q)/C(r,s),4个数均小于10000,答案不大于10^8 思路:根据答案的范围猜测,不需要使用高精度.根据唯一分解定理,每一个数都可以分解成若干素数相乘.先求出10000以内的所 ...

  9. 唯一分解定理 poj 1365

    一行代表一个数 x 给你底数和指数 求x-1的唯一分解定理的底数和指数 从大到小输出 #include<stdio.h> #include<string.h> #include ...

随机推荐

  1. python strip()方法使用

    描述 python strip() ,用于去除述字符串头尾指定字符(默认为空格或换行符)或字符序列. 注意:此方法只能去除头尾的空格或是换行符,不能去除中间的. 语法: str.strip([char ...

  2. python中常见的报错信息

    python中常见的报错信息 在运行程序时常会遇到报错提示,报错的信息会提示是哪个方向错的,从而帮助你定位问题: 搜集了一些python最重要的内建异常类名: AttributeError:属性错误, ...

  3. codeforces 922 B. Magic Forest(枚举、位运算(异或))

    题目链接:点击打开链接 Imp is in a magic forest, where xorangles grow (wut?) A xorangle of order n is such a no ...

  4. 《Python学习手册 第五版》 -第1章 问答环节

    第一章的主要内容是解疑答惑的,这个部分也是很适合初学者的,回答了大部分初学者所关注的问题 1.为什么使用Python,或者说Python的优点, 作者是分为两个部分来谈的,人们使用中的经验总结以及Py ...

  5. js + php服务器推送see(自定义推送时间)

    //javascript code var source = new EventSource("./main.php"); source.onmessage=function(ev ...

  6. 杭电1007-----C语言实现

    这道题花了好久的时间才做出来,刚开始没有思路,最后看了网上的解答,好难得样子,每次都没有看完,但是掌握了大概思想,今天试着做了一下,已ac 主要思想:先将点对按照x排序,再在x排好序的基础上按照y来排 ...

  7. 国产安全自主可控IT智能运维管理解决方案

    新一轮科技革命和产业变革席卷全球,大数据.云计算.物联网.人工智能.区块链等新技术不断涌现,数字经济正深刻地改变着人类的生产和生活方式,作为经济增长新动能的作用日益凸显.伴随增长的,还有网络中不断涌现 ...

  8. 实训第八天 有关python orm 的学习记录 常用方法01

    沿用第七天的数据库,数据库现在是这样的: 配置好主路由include子路由 子路由引入views 在views页面定义test测试请求如下: def test(request): # 1.all()方 ...

  9. Go1.14发布了,快来围观新的特性啦

    如期而至,Go1.14发布了,和往常一样,该版本保留了Go 1兼容性的承若,这个版本的大部分更新在工具链 .运行时库的性能提升方面,总的来说,还是在已有的基础上不断优化提成,大家期待的泛型还没有到来, ...

  10. StarUML之一、UML的相关基本概念

    为什么用UML 项目需要,在项目开发实现前期进行框架技术设计(条条大路通罗马通罗马,画图或者写代码都可以,适合就可以!). 用户的交互我们通常用Axure(原型设计)体现, 框架和功能结构设计则用UM ...