%I#94 2023年11月14日16:51:09

%S 0,0,2,3,9,9,20,21,32,37,54,50,77,81,96105135170168199217，

%电话：252240294309338350405393464655135415825756681724，

%电话：730819807902906957100910801052116811821254128013771365

%N反西格玛（N）：不除以N的小于N的数之和。

%C a（n）是n的适当非除数之和，三角形A173541中的行和_Omar E.Pol_，2010年5月25日

%如果n在A076617中，则C a（n）可被A000203（n）整除_伯纳德·肖特，2022年4月12日

%H N.J.A.Sloane，N表，N=1..10000的A（N）</a>

%F a（n）=n*（n+1）/2-西格玛（n）=A000217（n）-A000203（n）。

%F a（n）=A024916（n-1）-A153485（n），n>1_Omar E.Pol，2014年6月24日

%F From _Wesley Ivan Hurt_，2014年7月16日，2015年12月28日：（开始）

%F a（n）=总和{i=1..n}i*（天花板（n/i）-地板（n/i））。

%F a（n）=和{k=1..n}（n模k）+（-n模k）。（结束）

%F G.F.：x/（1-x）^3-和{k>=1}k*x^k/（1-x^k）.-_伊利亚·古特科夫斯基，2017年9月18日

%F From _Omar E.Pol_，2021年3月21日：（开始）

%F a（n）=A244048（n）+A004125（n）。

%F a（n）=A153485（n-1）+A004125（n），n>=2。（结束）

%对于p素数，F a（p）=（p-2）*（p+1）/2_伯纳德·肖特，2022年4月12日

%e a（12）=50等于5+7+8+9+10+11=50（1、2、3、4、6不包括在内，因为它们除以12）。

%p A024816:=程序（n）

%pn*（n+1）/2-数字理论[σ]（n）；

%p end程序：#_R.J.Mathar_，2013年8月3日

%t表[n（n+1）/2-除数Sigma[1，n]，{n，55}]（*_Robert G.Wilson v_*）

%t表[Total[Complement[Range[n]，Divisors[n]]]，{n，60}]（*_Harvey P.Dale_，2012年9月23日*）

%t使用[{nn=60}，#[[1]]-#[2]]&/@Thread[{Accumulate[Range[nn]]，DivisorSigma[1，Range[nn]}]]（*哈维·P·戴尔，2014年11月22日*）

%o（PARI）a（n）=n*（n+1）/2-西格玛（n）\\查尔斯·格里特豪斯IV_，2012年3月19日

%o（哈斯克尔）

%o a024816=总和。a173541_row---Reinhard Zumkeller，2014年2月19日

%o（Magma）[n*（n+1）div 2-SumOfDivisors（n）：n in[1..60]]；//_Vincenzo Librandi_，2015年12月29日

%o（Python）

%o从sympy导入divisor_sigma

%o定义A024816（n）：返回（n*（n+1）>>1）-divitor_sigma（n）#_Chai Wah Wu_，2023年4月28日

%o（SageMath）

%o def A024816（n）：返回和（k代表k in（0..n-1），如果不是k除（n））

%o打印（[A024816（n）for n in srange（1，55）]）#_Peter Luschny_，2023年11月14日

%Y参考A000203（西格玛）、A000217、A004125、A023896、A024916、A066760、A076617、A153485、A173539、A173540、A173541、A244048、A352810、A35281。

%Y参见A342344（用于对称表示）。

%放松，不，很好

%氧1,3

%Paul Jobling（Paul.Jobling（AT）whitecross.com）