题意:

解法:

f a ( x ) 的 逆 函 数 是 l o g a b a n s = ∑ a = 2 n ( a ∑ b = a n ⌊ f a − 1 ( b ) ⌋ ⌈ f b − 1 ( a ) ⌉ ) = ∑ a = 2 n ( a ∑ b = a n ⌊ l o g a ( b ) ⌋ ⌈ l o g b ( a ) ⌉ ) 当 b > = a 时 , ⌈ l o g b ( a ) ⌉ = 1 恒 成 立 那 么 a n s = ∑ a = 2 n ( a ∑ b = a n ⌊ l o g a ( b ) ⌋ 当 b < a 2 时 , ⌊ l o g a ( b ) ⌋ = 1 1. 当 a > n 时 , a 2 > n , 由 于 b < = n , 因 此 b < a 2 恒 成 立 , 此 时 ⌊ l o g a ( b ) ⌋ = 1 恒 成 立 , 此 时 a n s = ∑ a = n + 1 n a ∑ b = a n 1 = ∑ a = n + 1 n a ∗ ( n − a + 1 ) = ∑ a = n + 1 n ( n ∗ a − a 2 + a ) = n ∑ a = n + 1 n a − ∑ a = n + 1 n a 2 + ∑ a = n + 1 n a = ( n + 1 ) ∑ a = n + 1 n a − ∑ a = n + 1 n a 2 发 现 两 项 分 别 是 1 次 方 前 n 项 和 与 2 次 方 前 n 项 和 , 可 以 O ( 1 ) 计 算 . 2. 当 a < = n 时 , O ( s q r t n ) = 1 e 6 枚 举 a , 式 子 内 层 O ( l o g ) 计 算 即 可 , 复 杂 度 O ( 1 e 6 ∗ l o g ) . f_a(x)的逆函数是log_{a}b\\ ans=\sum_{a=2}^n(a\sum_{b=a}^n\lfloor f_a^{-1}(b)\rfloor\lceil f_b^{-1}(a)\rceil)\\ =\sum_{a=2}^n(a\sum_{b=a}^n\lfloor log_a(b)\rfloor\lceil log_b(a)\rceil)\\ 当b>=a时,\lceil log_b(a)\rceil=1恒成立\\ 那么ans=\sum_{a=2}^n(a\sum_{b=a}^n\lfloor log_a(b)\rfloor\\ 当b<a^2时,\lfloor log_a(b)\rfloor=1\\ 1.当a>\sqrt n时,a^2>n,由于b<=n,因此b<a^2恒成立,\\ 此时\lfloor log_a(b)\rfloor=1恒成立,\\ 此时ans=\sum_{a=\sqrt n+1}^na\sum_{b=a}^n1=\sum_{a=\sqrt n+1}^na*(n-a+1)\\ =\sum_{a=\sqrt n+1}^n (n*a-a^2+a)\\ =n\sum_{a=\sqrt n+1}^n a-\sum_{a=\sqrt n+1}^n a^2+\sum_{a=\sqrt n+1}^n a\\ =(n+1)\sum_{a=\sqrt n+1}^n a-\sum_{a=\sqrt n+1}^n a^2\\ 发现两项分别是1次方前n项和与2次方前n项和,可以O(1)计算.\\ 2.当a<=\sqrt n时,O(sqrt n)=1e6枚举a,式子内层O(log)计算即可,\\ 复杂度O(1e6*log). fa​(x)的逆函数是loga​bans=a=2∑n​(ab=a∑n​⌊fa−1​(b)⌋⌈fb−1​(a)⌉)=a=2∑n​(ab=a∑n​⌊loga​(b)⌋⌈logb​(a)⌉)当b>=a时,⌈logb​(a)⌉=1恒成立那么ans=a=2∑n​(ab=a∑n​⌊loga​(b)⌋当b<a2时,⌊loga​(b)⌋=11.当a>n ​时,a2>n,由于b<=n,因此b<a2恒成立,此时⌊loga​(b)⌋=1恒成立,此时ans=a=n ​+1∑n​ab=a∑n​1=a=n ​+1∑n​a∗(n−a+1)=a=n ​+1∑n​(n∗a−a2+a)=na=n ​+1∑n​a−a=n ​+1∑n​a2+a=n ​+1∑n​a=(n+1)a=n ​+1∑n​a−a=n ​+1∑n​a2发现两项分别是1次方前n项和与2次方前n项和,可以O(1)计算.2.当a<=n ​时,O(sqrtn)=1e6枚举a,式子内层O(log)计算即可,复杂度O(1e6∗log).

code:

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int maxm=2e6+5;
const int mod=998244353;
int n;
int ppow(int a,int b,int mod){int ans=1%mod;a%=mod;while(b){if(b&1)ans=ans*a%mod;a=a*a%mod;b>>=1;}return ans;
}
int cal(int n){//1次方前n项和return n%mod*((n+1)%mod)%mod*ppow(2,mod-2,mod)%mod;
}
int cal2(int n){//2次方前n项和return n%mod*((n+1)%mod)%mod*((n*2+1)%mod)%mod*ppow(6,mod-2,mod)%mod;
}
void solve(){cin>>n;int ans=0;//a*a<=nfor(int i=2;i*i<=n;i++){for(int l=i,t=1;l<=n;l*=i,t++){int r=min(n,l*i-1);ans+=(r-l+1)*t%mod*i%mod;ans%=mod;}}//a*a>n;int a=1;while(a*a<=n)a++;ans+=((n+1)%mod)*((cal(n)-cal(a-1))%mod)%mod;ans-=(cal2(n)-cal2(a-1))%mod;//ans=(ans%mod+mod)%mod;cout<<ans<<endl;
}
signed main(){ios::sync_with_stdio(0);solve();return 0;
}

2019ICPC银川 F.Function!(数学)相关推荐

  1. F#,印度数学教育与优秀编程

    F#,印度数学教育与优秀编程 gisspace 转帖 这个周末,我终于可以坐下来做我在一年来一直承诺要做的事情:我开始应用F#.这次体验相当有趣.十多年来,我第一次不能只是通过坐着瞧一些样本代码来理解 ...

  2. F. Function!

    题目:F. Function! 总结:对于b>=a,logba往上取,那么一定是1.那么就只需要考虑logab,q =sqrt(n),对于 q+1~n,每一段的答案都是 i*(n-i+1),可以 ...

  3. 2019ICPC(银川) - Function!(数论+数学分块)

    题目链接:点击查看 题目大意:给出公式,其中,则,现在给出n(<=1e12),求出答案对998244353取模后的答案 题目分析:若用暴力实现上述公式,只需要两层for循环,时间复杂度为n*n, ...

  4. [转]F#,印度数学教育与优秀编程

    这个周末,我终于可以坐下来做我在一年来一直承诺要做的事情:我开始应用F#.这次体验相当有趣.十多年来,我第一次不能只是通过坐着瞧一些样本代码来理解一门语言. 这次体验强化了当前通用编程语言的可悲状况. ...

  5. 2019icpc银川打铁站 赛后总结

    经过 10.20上午9点比赛开始,非常紧张,瞄了一眼A.B,果断觉得B更简单,开始读B,是个签到性质的思维题,给每行.每列的数分别加上一个非负整数,遮盖其中一个数,问这个数是多少.由于过于紧张,B大概 ...

  6. CodeForces - 1301C Ayoub's function(数学)

    题目链接:点击查看 题目大意:规定函数 f(s) 为01字符串 s 中至少包含一个 1 的所有子串的数量,比如 f(10101) = 12 ,其中十二个子串分别为(1,2),(1,3),(1,4),( ...

  7. 2019ICPC(银川) - Take Your Seat(概率公式)

    题目链接:点击查看 题目大意:朵哈要坐飞机去旅游,出发的时候包括朵哈在内一共有n个人要上飞机,这班飞机有n个座位,第i个人的座位号是i,朵哈的座位号是1. 现在n个人按照座位号是从1到n的顺序上飞机, ...

  8. [美赛F奖][数学建模][经验贴]2021美赛F奖的那些事

    写在前面 2021美赛都过去半年了,一直也在忙各种各样的事情,刚好上学期有一门项目管理的课程,课程论文写的就是美赛经验,偷个懒直接改下排版复制上来,以作留念 贴个奖状~~ 一.引言 项目是一个组织为实 ...

  9. 2023牛客寒假算法基础集训营2 -- E-Tokitsukaze and Function(数学 二分)

    题目如下: T o k i t s u k a z e Tokitsukaze Tokitsukaze 有一个函数 f ( x ) = ⌊ x n ⌋ + x − 1 f(x)=⌊\frac{x}{n ...

最新文章

  1. 线程池ThreadPoolExecutor
  2. 混合办公时代来了?携程试点每周两天居家办公,76%的员工主动报名!
  3. linux命令:返回上一次目录
  4. 【Spring Boot 分享】开源项目【8个】
  5. 7-7 字符串修改 (15 分)
  6. 设置linearlayout最大高度_技术案例 | 排烟口个数与挡烟垂壁高度的关系探讨
  7. Spring异常重试框架Spring Retry
  8. java初级面试总结
  9. 数字化转型:中国企业数字化转型趋势
  10. webstorm 使用技巧
  11. 看完就学会系列,小小一篇文章教会你利用Python网络爬虫抓取王者荣耀图片(建议收藏)
  12. ckplayer html5 添加广告,ewebeditor下利用ckplayer增加html5 (mp4)全平台的支持
  13. MATLAB小波变换图像处理简单示例
  14. 【WangEditor】使用富文本编辑器 WangEditor 实现用户自定义图片大小(改js源码)
  15. python关机usb端口_使用python关闭windows中的特定USB端口
  16. Python基础嘚吧嘚
  17. 《工程学导论》读书笔记第二章工程与科学
  18. 科技新品 | 卡西欧高性能金属腕表;佳能大幅面打印机;三星2亿像素传感器
  19. 用安卓手机远程管理linux,支持SSH密钥登录
  20. 软件测试pdca循环的好处,认识PDCA循环重要性

热门文章

  1. 学习笔记69—金蝶财务软件安装教程(KIS12.3,win10)
  2. Ionic+Angular 开发移动应用一:环境搭建
  3. python-图片批量处理成正方形
  4. 2021-10-27
  5. 那些40岁的程序员都去哪了
  6. ODOO10 富文本编辑器SummerNote输入联想短句,输入快捷短语
  7. 第二讲-单元格格式设置
  8. Ubuntu开启,关闭防火墙
  9. cast计算机英语,机械类英语词汇汇总
  10. 解决Deepin系统QQ图片不显示的问题