2021牛客寒假算法基础集训营1 J. 一群小青蛙呱蹦呱蹦呱(lcm思维)
发布日期:2021-06-30 10:28:23
浏览次数:2
分类:技术文章
本文共 1310 字,大约阅读时间需要 4 分钟。
直接写不好写
考虑 l c m lcm lcm是取每个数分解质因子后,每种质因子的最高次幂
那我们看,青蛙会把所有只包含一种质因子的数筛掉
那么对于质因子 p p p,包含最多幂 p p p的数显然是 2 p x < = n 2p^x<=n 2px<=n
解出这个 x x x就是对 l c m lcm lcm的贡献
特殊的, p = 2 p=2 p=2时,接触 3 p x < = n 3p^x<=n 3px<=n的最大 x x x
#includeusing namespace std;typedef long long ll;const int mod = 1e9+7;const int maxn = 160000000;bool vis[maxn+10],ok[maxn+10];int prime[maxn],top,n;void make_prime(){ for(int i=2;i<=maxn;i++) { if( !vis[i] ) prime[++top] = i; for(int j = 1;i*prime[j]<=maxn;j++) { vis[i*prime[j]] = 1; if( i%prime[j]==0 ) break; } }}int gcd(int a,int b){ return b==0?a:gcd(b,a%b); }int quick(int x,int n){ int ans = 1; for( ; n ; n>>=1,x=1ll*x*x%mod ) if( n&1 ) ans = 1ll*ans*x%mod; return ans;}int lcm(int a,int b){ return 1ll*a*b%mod*quick(gcd(a,b),mod-2)%mod; }int main(){ cin >> n; make_prime(); int ans = 1; for(int j=2;j<=top;j++) { ll now = 2*prime[j]; int ci = 0; while( now<=n ) { ci++; now = now*prime[j]; ans = 1ll*ans*prime[j]%mod; } } //7的话是6 ll now = 3*2; while( now<=n ) { now = now*2; ans = 1ll*ans*2%mod; } if( ans==1 ) cout << "empty"; else cout << ans;}//如果刚好是素数的次幂,会被筛选掉//
转载地址:https://issue-is-vegetable.blog.csdn.net/article/details/113577102 如侵犯您的版权,请留言回复原文章的地址,我们会给您删除此文章,给您带来不便请您谅解!
发表评论
最新留言
能坚持,总会有不一样的收获!
[***.219.124.196]2024年04月06日 22时24分30秒
关于作者
喝酒易醉,品茶养心,人生如梦,品茶悟道,何以解忧?唯有杜康!
-- 愿君每日到此一游!
推荐文章
159.用ASP实现无组件上传/下载文件
2019-04-30
160.可以生成标准的.sql文件,在mysql下导入
2019-04-30
161.SQL Server到Oracle连接服务器的实现
2019-04-30
162.SQL Server到SQLBASE连接服务器的实现
2019-04-30
163.SQL Server到SYBASE连接服务器的实现
2019-04-30
164.二进制文件存取-案例
2019-04-30
165.作业使用案例
2019-04-30
166.常用作业定义的T-SQL模板
2019-04-30
167.使用作业定时启动婷数据库的案例
2019-04-30
168.使用作业异步调用存储过程的案例
2019-04-30
169.实现秒级作业的案例
2019-04-30
170.解决计算机名修改或作业移植导致的服务名称问题
2019-04-30
171.操作SQL SERVERAGENT服务的扩展存储过程
2019-04-30
rtf格式的一些说明,转载的,我找到的rtf资料中比较实用的一片文章了
2019-04-30
RTF文件格式编码说明
2019-04-30
RTF 语法1
2019-04-30
RTF文件格式研究报告
2019-04-30
RichEdit的用法总结
2019-04-30
BCB 中测量Richedit 的文本总行高
2019-04-30
C++ Builder组件属性详解
2019-04-30