ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

HDU 5608 function(杜教筛)

2021-06-13 14:04:29  阅读:261  来源: 互联网

标签:function dn HDU frac cdot sum rfloor 杜教 5608


目录

链接

描述

  • 已 知 N 2 − 3 N + 2 = ∑ d ∣ N f ( d ) , 求 ∑ i = 1 n f ( i ) . 已知 N^2-3N+2 = \sum_{d|N}f(d), 求\sum_{i=1}^{n}f(i). 已知N2−3N+2=∑d∣N​f(d),求∑i=1n​f(i).

分析

杜教筛

  • 迪立克利卷积 g ∗ f g*f g∗f

( g ∗ f ) ( n ) = ∑ d ∣ n g ( d ) ⋅ f ( n d ) = ∑ d ∣ n g ( n d ) ⋅ f ( d ) (g*f)(n) = \sum_{d|n}g(d)\cdot f(\frac{n}{d}) = \sum_{d|n}g(\frac{n}{d})\cdot f(d) (g∗f)(n)=d∣n∑​g(d)⋅f(dn​)=d∣n∑​g(dn​)⋅f(d)

  • 杜教筛

∑ i = 1 n ( g ∗ f ) ( i ) = ∑ i = 1 n ∑ d ∣ i g ( d ) ⋅ f ( i d ) = ∑ d = 1 n g ( d ) ⋅ S ( ⌊ n d ⌋ ) \sum_{i=1}^{n}(g*f)(i)= \sum_{i=1}^{n}\sum_{d|i}g(d)\cdot f(\frac{i}{d})=\sum_{d=1}^{n}g(d)\cdot S(\lfloor\frac{n}{d}\rfloor) i=1∑n​(g∗f)(i)=i=1∑n​d∣i∑​g(d)⋅f(di​)=d=1∑n​g(d)⋅S(⌊dn​⌋)
= g ( 1 ) ⋅ S ( n ) + ∑ d = 2 n g ( d ) ⋅ S ( ⌊ n d ⌋ ) ⇒ =g(1)\cdot S(n)+\sum_{d=2}^{n}g(d)\cdot S(\lfloor\frac{n}{d} \rfloor) \Rightarrow =g(1)⋅S(n)+d=2∑n​g(d)⋅S(⌊dn​⌋)⇒
g ( 1 ) ⋅ S ( n ) = ∑ i = 1 n ( g ∗ f ) ( i ) − ∑ d = 2 n g ( d ) ⋅ S ( ⌊ n d ⌋ ) g(1)\cdot S(n)=\sum_{i=1}^{n}(g*f)(i)-\sum_{d=2}^{n}g(d)\cdot S(\lfloor\frac{n}{d} \rfloor) g(1)⋅S(n)=i=1∑n​(g∗f)(i)−d=2∑n​g(d)⋅S(⌊dn​⌋)

建模

  • 条件

∑ d ∣ N f ( d ) ⋅ 1 = ( f ∗ I ) ( n ) = ( I ∗ f ) ( n ) = ∑ d ∣ N 1 ⋅ f ( n d ) , 式 中 I ( n ) = 1 \sum_{d|N}f(d)\cdot 1 = (f*I)(n)=(I*f)(n)=\sum_{d|N}1\cdot f(\frac{n}{d}) ,式中I(n)=1 d∣N∑​f(d)⋅1=(f∗I)(n)=(I∗f)(n)=d∣N∑​1⋅f(dn​),式中I(n)=1
即 : N 2 − 3 N + 2 = ∑ d ∣ N 1 ⋅ f ( n d ) 即:N^2-3N+2=\sum_{d|N}1\cdot f(\frac{n}{d}) 即:N2−3N+2=d∣N∑​1⋅f(dn​)

  • 结论

S ( N ) = ∑ i = 1 N ( i 2 − 3 i + 2 ) − ∑ d = 2 N 1 ⋅ S ( ⌊ N d ⌋ ) = ∑ i = 1 n ( i 2 − 3 i + 2 ) − ∑ d = 2 N S ( ⌊ N d ⌋ ) S(N) =\sum_{i =1}^{N}(i^2-3i+2)-\sum_{d=2}^{N}1\cdot S(\lfloor\frac{N}{d} \rfloor)=\sum_{i =1}^{n}(i^2-3i+2)-\sum_{d=2}^{N}S(\lfloor\frac{N}{d} \rfloor) S(N)=i=1∑N​(i2−3i+2)−d=2∑N​1⋅S(⌊dN​⌋)=i=1∑n​(i2−3i+2)−d=2∑N​S(⌊dN​⌋)

标签:function,dn,HDU,frac,cdot,sum,rfloor,杜教,5608
来源: https://blog.csdn.net/jpphy0/article/details/117874271

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有