首页 > 其他分享 >CSP模拟-19

CSP模拟-19

时间:2023-08-12 15:12:38浏览次数:25  
标签:lfloor ... lcm 19 ll 整除 include CSP 模拟

前言

emm.....考场其实想到T2正解的思路了,但是不会优化,导致有些拉胯少了50分。还有就是说数学题我是真不行,向上次 \(fengwu\) 的筛我不会,这会最简单的容斥想这么老半天学这么老半天都不会,着实是有些废物了。

T1 十年之约

一道很简单的数学题QAQ,但我就是不会,我真服了。

首先对于 \(f(i)=k\),则我们的 \(i\) 一定能被 \(1,2,...,k-1\) 的所有数整除并且一定不能被 \(k\) 整除,于是就有了我们下面的式子:\(lcm(1,2,3,...,k-1) \mid i\) 且 \(k \nmid i\)。

对于一个集合里面任意的数 \(i\) 在满足 \(lcm(1,2,3,...,k-1) \mid i\) 的情况下它只会有两种情况,能被 \(k\) 整除和不能被 \(k\) 整除。所以用一个小小的容斥就能求出来题目要求的 \(k\) 的个数啦 \(k\) :\(\left\lfloor \frac{i}{lcm(1,2,3,...,k-1)} \right\rfloor - \left\lfloor \frac{i}{lcm(1,2,3,...,k)} \right\rfloor\)。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>

using namespace std;
typedef long long ll;
const ll mod = 1000000007;

ll t, n;

inline void solve() {
    ll lcm = 1, last, ans = 0;
    for (ll i = 2; lcm <= n; ++ i) {
        last = lcm;
        lcm = lcm * i / __gcd(lcm, i);
        ans = (ans + i * (n / last - n / lcm) % mod + mod) % mod;
    }
    printf("%lld\n", ans);
}

int main() {
    scanf("%lld", &t);
    while (t --) {
        scanf("%lld", &n);
        solve();
    }
    return 0;
}

标签:lfloor,...,lcm,19,ll,整除,include,CSP,模拟
From: https://www.cnblogs.com/jueqingfeng/p/17624823.html

相关文章

  • Linux磁盘故障,模拟故障及解决思路方法
    每个分区起始位置都有一个inod表索引节点表(类似于目录表)每一个文件都对应一个编号称为索引节点,如果这个空间文件数太多了,记满了,就说明索引节点表耗尽。故障1 该分区不能正常读写或者说只能读不能写了但是又没有满,就代表文件系统有问题,文件系统有问题需要进行修复命令:故障2:索引......
  • SpringBoot复习:(19)Condition接口和@Conditional注解
    Condition接口代码如下:publicinterfaceCondition{ booleanmatches(ConditionContextcontext,AnnotatedTypeMetadatametadata);}它是一个函数式接口,只有一个方法matches用来表示条件是否满足。matches方法中的ConditionContext类对象context可以通过getEnvironment方法获......
  • 2023.8.11 模拟赛
    A询问\(L\lei,j\leR\),其中\(\gcd(i,j)\not=1,i,j\)的对数。莫反先求出\(gcd(i,j)\not=1\)的对数,然后再直接调和级数暴力删去\(i,j\)是倍数的对数即可。BP4334[COI2007]Policija考虑建立圆方树。圆方树是怎么样的呢?圆方树是对于每个点双,都建立一个方点,然后......
  • 19.17RU安装问题汇总
    问题概述19.17RU安装问题汇总一、lib库被其他用户使用二、CRS-1159:Theclustercannotbesettorollingpatchmode三、NoreadorwritepermissiontoORACLE_HOME/.patch_storage四、Datapatch:couldn'topenencmapgbk.enc五、CRS-6706:OracleClusterwareReleasepatch......
  • IIS8.5 Error Code 0x8007007e HTTP 错误 500.19的解决方法
    windowserver2012R2IIS8.5引用:https://www.52jbj.com/yunying/340443.htmlHTTP错误500.19-InternalServerError  无法访问请求的页面,因为该页的相关配置数据无效。    详细错误信息    模块DynamicCompressionModule    通知SendResponse    处......
  • 谷歌2023年4月19日最新更新规则及算法
    多年来,谷歌的的核心排名系统一直致力于奖励能够提供良好网页体验的内容,正如谷歌在 2011年提供的2019年更新,并于去年加入创建实用、可靠、以用户为中心的内容帮助页面的指南中所述。该帮助页面是谷歌的搜索要素的关键资源。会定期推荐希望使用Google搜索取得理想成效的用户,仔......
  • VS2019 C++ 调用python函数/类对象的方法
    1.环境配置VS工程配置要和python一致,安装的python如果是64位的,工程配置也要选成64位的在工程配置中添加包含目录和库目录,添加python环境目录里的include和libs文件夹路径。想要运行的keras-yolo3是在Anaconda中配置的环境,所以相应的文件夹路径可以在Anaconda的环境文件中......
  • 基于matlab模拟RADAR预警雷达
    ✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信。......
  • 2023下半年产品经理NPDP认证8月19日开班
    产品经理国际资格认证NPDP是新产品开发方面的认证,集理论、方法与实践为一体的全方位的知识体系,为公司组织层级进行规划、决策、执行提供良好的方法体系支撑。  【认证机构】 产品开发与管理协会(PDMA)成立于1979年,是全球范围内产品开发与管理专业人士最杰出的倡导者,协助个人、企业......
  • 2023年CSPM-3国标项目管理中级认证报名到这里错不了
    CSPM-3中级项目管理专业人员评价,是中国标准化协会(全国项目管理标准化技术委员会秘书处),面向社会开展项目管理专业人员能力的等级证书。旨在构建多层次从业人员培养培训体系,建立健全人才职业能力评价和激励机制的要求,培养我国项目管理领域复合型人才。  【证书含金量】 ·竞聘优先......