Topcoder SRM 661 (Div.1) 250 MissingLCM - 数论
【题意】
给你一个数N(1<=N<=10^6),要求最小的M(M>N),使得lcm(n+1,n+2,...m)=lcm(1,2,3,...,m)
【思路】
手速太慢啦,等敲完代码的时候发现比赛已经结束了
一开始我想直接枚举m,并判断lcm(1,..,m)与lcm(n+1,n+2,...,m)是否相等,但发现,当求到lcm(1,...,40)的时候就爆LL了
显然不能这样求
也就是说,要求出具体lcm(1,2,...,m)的值是很困难的
怎么求
可以把它分解质因数,分解成几个质数相乘的形式
判断lcm(1,...,m)和lcm(n+1,n+2,...,m)的质因数是否完全一样。
但是仅仅1~1000000的质数有8万个
枚举m再枚举质数显然吃不消。
然而我注意到有一条性质(不知道算不算)
假设有质数K,可以求出t,使得K的t次方刚好小于m(K^t<=m)
那么lcm(1,2,...,m)分解质因数中一定而且最多有t个质数K连乘,
这样就可以很快地吧lcm(1,2,...,m)分解质因数
那么怎么把lcm(n+1,n+2,...,m)分解质因数呢
仍然假设质数K,可以求出最大的t,以及一个常数c(1<=c<K),使得 n+1<=c*K^t<=m
那么lcm(n+1,n+2,...,m)分解质因数中一定而且最多有t个质数K连乘。
比如说质数3,n=16,m=22,可以求的c=2,t=2,即17<=2*3^2=18<=22,这样lcm(17,18,19,20,21,22)中最多有2个质数3连乘
既然知道怎么求lcm(n+1,n+2,..,m)和lcm(1,2,..,m)了
来探讨一下怎么求最小的m吧
我们想让这两个lcm分解质因数后完全一样,也就是说连乘的质数个数也完全相等。
也就是说,对于每个质数K都可以满足,存在c和最大的t 使得n+1<=c*K^t<=m
对于大于n小于m的质数,我们假设是P,那么一定n+1<=P<=m,一定可以满足条件
所以我们就只看小于等于n的质数就可以了
因为要使每个小于N的质数K,都存在c和最大的t 使得n+1<=c*K^t<=m,
我们枚举每一个质数,并求得c和t,使得刚好c*K^t>=n,
答案就取最大的c*K^t,即 max( c*K^t )
这样lcm(1,2,...,m)和lcm(n+1,n+2,...,m)的分解质因数后均至少有t个质数K。
如果最终m有 k^(t+1)<=m,那么这个K^(t+1)>n一定成立,故仍满足条件
#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<set>
#include<map>
#include<stack>
#include<vector>
#include<queue>
#include<string>
#include<sstream>
#define eps 1e-9
#define ALL(x) x.begin(),x.end()
#define INS(x) inserter(x,x.begin())
#define FOR(i,j,k) for(int i=j;i<=k;i++)
#define MAXN 1005
#define MAXM 40005
#define INF 0x3fffffff
using namespace std;
typedef long long LL;
LL i,j,k,n,m,x,y,T,big,cas,num;
bool flag; LL cur,ans; bool prim[];
LL ver[];
void GetPrim(LL size)
{
LL m=sqrt(size+0.5);
memset(prim,,sizeof(prim));//可以根据情况进行清空操作
num=;//把找到的质数存入ver数组中,num为ver数组的长度 //如果要获得质数数组,i就枚举到size,如果仅仅是prim数组,就枚举到m
for (LL i=;i<=size;i++)
{
if (!prim[i])
{
ver[++num]=i;
if (i<=m) for (LL j=i*i;j<=size;j+=i) prim[j]=;
}
}
} class MissingLCM
{
public:
int getMin(int N)
{
LL n=N;
GetPrim(n);
LL ans=n+;
for (i=num;i>=;i--)
{
LL u=ver[i];
for (j=;j*u<=n;j*=u); ans=max(ans,(n/j+)*j);
}
return ans;
}
};
Topcoder SRM 661 (Div.1) 250 MissingLCM - 数论的更多相关文章
- Topcoder SRM 656 (Div.1) 250 RandomPancakeStack - 概率+记忆化搜索
最近连续三次TC爆零了,,,我的心好痛. 不知怎么想的,这题把题意理解成,第一次选择j,第二次选择i后,只能从1~i-1.i+1~j找,其实还可以从j+1~n中找,只要没有被选中过就行... [题意] ...
- TopCoder SRM 596 DIV 1 250
body { font-family: Monospaced; font-size: 12pt } pre { font-family: Monospaced; font-size: 12pt } P ...
- 【topcoder SRM 702 DIV 2 250】TestTaking
Problem Statement Recently, Alice had to take a test. The test consisted of a sequence of true/false ...
- TopCoder SRM 560 Div 1 - Problem 1000 BoundedOptimization & Codeforces 839 E
传送门:https://284914869.github.io/AEoj/560.html 题目简述: 定义"项"为两个不同变量相乘. 求一个由多个不同"项"相 ...
- TopCoder SRM 667 Div.2题解
概览: T1 枚举 T2 状压DP T3 DP TopCoder SRM 667 Div.2 T1 解题思路 由于数据范围很小,所以直接枚举所有点,判断是否可行.时间复杂度O(δX × δY),空间复 ...
- 竞赛图的得分序列 (SRM 717 div 1 250)
SRM 717 DIV 1 中 出了这样一道题: 竞赛图就是把一个无向完全图的边定向后得到的有向图,得分序列就是每个点的出度构成的序列. 给出一个合法的竞赛图出度序列, 要求构造出原图(原题是求(u, ...
- [topcoder]SRM 646 DIV 2
第一题:K等于1或者2,非常简单.略.K更多的情况,http://www.cnblogs.com/lautsie/p/4242975.html,值得思考. 第二题:http://www.cnblogs ...
- [topcoder]SRM 633 DIV 2
第一题,http://community.topcoder.com/stat?c=problem_statement&pm=13462&rd=16076 模拟就可以了. #includ ...
- Topcoder SRM 648 (div.2)
第一次做TC全部通过,截图纪念一下. 终于蓝了一次,也是TC上第一次变成蓝名,下次就要做Div.1了,希望div1不要挂零..._(:зゝ∠)_ A. KitayutaMart2 万年不变的水题. # ...
随机推荐
- 简单易懂, JUnit 框架问答
本文算是一个关于Junit4相关的知识分享,但是不同于网上大段的源码分析,模式学习文章,我想通过问答的形式,引出代码来简明阐述JUnit4是如何实现需要的功能的. 考虑到任何一个框架,都是为了解决问题 ...
- 写个自动安装JDK的shell脚本
#!/bin/bash ################################################# # # INSTALL JDK AUTOMATICALLY # # auth ...
- [BZOJ 1018] [SHOI2008] 堵塞的交通traffic 【线段树维护联通性】
题目链接:BZOJ - 1018 题目分析 这道题就说明了刷题少,比赛就容易跪..SDOI Round1 Day2 T3 就是与这道题类似的..然而我并没有做过这道题.. 这道题是线段树维护联通性的经 ...
- C语言嵌入式系统编程修炼之六:性能优化
使用宏定义 在C语言中,宏是产生内嵌代码的唯一方法.对于嵌入式系统而言,为了能达到性能要求,宏是一种很好的代替函数的方法. 写一个"标准"宏MIN ,这个宏输入两个参数并返回较小的 ...
- android 自定义命名空间
一.统一的用户界面是可以使得应用程序更友好.要做到用户界面的统一,我们就必须用到风格(style)和主题(theme).自定义一个View的方法步骤如下:1.首先,在values文件夹下定义一个att ...
- oracle的存储过程语法(转)
1.ORA-00942: table or view does not exist 指的你要操作的表尚未存在,需要先create出来先. 2.ORA-00922: missing or invalid ...
- leetcode面试准备:Decode Ways
1 题目 A message containing letters from A-Z is being encoded to numbers using the following mapping: ...
- TCP/IP 目录导航
用了近二十天的时间,把一本800页的书看完,感觉收获还是很大的.对网络,对这些协议有了深刻的认识! 知道了路由器与交换机的区别. 知道了IP地址的特点. 知道了网络的分层,物理层,数据链路层,网络层, ...
- HDwiki文件上传导致远程代码执行漏洞
漏洞版本: HDwiki(2011) 漏洞描述: 互动维客开源系统(HDwiki)作为中国第一家拥有自主知识产权的中文维基(Wiki)系统,由互动在线(北京)科技有限公司于2006 年11月28日正式 ...
- 通过超链接打开App应用
URL schemes, 通过超链接打开App应用 var mobileAppInstall = (function () { var ua = navigator.userAgent, loadIf ...