数学建模社区-数学中国

标题: 一些初等算术逻缉函数 [打印本页]

作者: lilianjie    时间: 2012-1-11 17:58
标题: 一些初等算术逻缉函数
Z:=IntegerRing() ;
1 R7 ^; i' U, n# rZ;
" h4 s& `4 j( a7 A4 d9 yZ12:=IntegerRing(12); 4 C' ~7 l0 n! M
Z12;& M$ a6 A/ z. O$ {% ]1 \
Z17:=IntegerRing(19); : Y7 C: W  ^6 B6 ~( K
Z17;* {8 |* m  y5 P
n:=Z!100;+ B; Y- y* k5 W/ \1 \+ A
n;
0 {5 Z! z0 [+ wm:=Random(Z12);
% m$ p9 |* V  m" q2 z9 \& em;& ?- M' k3 \. p0 `+ M9 A

& o. g$ @! I; S. x9 o+ cn div m;    求余+ Z1 _" t* C1 A, k9 K" t
Z!n mod Z!m;
% O5 t- y3 ^3 o( a& nExactQuotient(n, 5);求商- }: a+ J$ H% y/ e, c
ShiftLeft(n, 1);乘2$ I' S2 W* T5 n1 i. A& e: f
ShiftLeft(n, 3);乘8
% h4 ]4 p& L2 Y3 F( eShiftLeft(n, 10);乘1024  P) u- X* L+ ^
ShiftRight(n, 2);除4
# P" ?6 |. e9 t/ m2 v, JShiftRight(n, 3);除8+ o0 Z8 p1 \) {: Q3 e
ShiftRight(n, 10);除1024
' I4 O% {3 z* p5 ?9 i. [ModByPowerOf2(n, 2);模4
: n5 b; D: e. }. w+ Z1 Y9 {ModByPowerOf2(n, 3);模8
0 G' u- J& H* G& L2 GModByPowerOf2(n, 5);模32
IsEven(n) ;9 m6 n9 ?# C# e6 ]1 @' C1 N
IsOdd(n);
( W8 V4 u" Y0 n7 c, hIsDivisibleBy(n, 13);/ K6 p1 [* _( y
IsDivisibleBy(n, 2);
3 x9 A! B# b6 X6 ~& n1 BIsSquarefree(n) ;" o2 C1 Z0 J- u* p7 M
IsPower(n);
5 E7 C% {  \/ a% \+ i: aIsPower(n, 10) ;
: d( y1 U& q+ A- N, d% xIsPower(n, 3) ;
  x9 y+ O5 l5 g# o0 [7 U3 N" SIsSquare(n);9 J1 y) k' f5 i' \' v: x8 q1 @
IsPrime(n) ;
) D0 n* T+ K- b* C" M9 X/ H5 cIsIntegral(n);! e& I$ G/ R7 t2 D0 V- Q$ v% C
IsSinglePrecision(n);单精度% z3 t, ~7 l% o; d% @

0 I/ o6 W5 z  y
! x! P0 R! X/ V! d: }0 P2 k8 {Integer Ring3 p/ P9 Z+ z; T$ F" s2 v
Residue class ring of integers modulo 12
. v( i. o( S1 q! ]/ h6 rResidue class ring of integers modulo 198 \3 T0 O1 X9 l& M! P
100- ^( f3 R. Y; Z9 W  q. h
9
% Z6 f, Q8 N0 d11
& D( Y* N9 Y$ g' Q' M1/ x# F! W& ~! ^/ E5 e0 K. {
20( v% I3 ^. P# K
200% d( s# B( V1 o4 L1 }: \7 `
800
# ?4 P# B" E$ `6 [$ N102400
* T9 e# L' G+ i8 r* V0 c25
$ Z% j5 z9 R/ d, y/ Q5 K1 e* Z12: w4 z% R, [5 t% r
0
1 }/ p" F" s0 K) F  [1 T0
" s8 [0 I, X( R* N/ ]& I4, J$ X: O- s* Y- M8 F
4+ Z/ Q8 F2 H6 w4 s2 J5 \& W9 k
true
7 _7 r$ ~. Y& hfalse  A  I2 ^' |1 I3 R8 u, D
false
. \) a: e+ }* r& M) {) ktrue
& H, ~* i. a& L0 U* O9 t7 Xfalse4 J" g! q& l$ J6 e2 G' u  h
true 10 2
7 t: p: K, S5 U" \5 Y+ r6 ~* kfalse
3 K9 W7 Y6 v( i6 |! T1 l& O0 ^4 Mfalse
$ @2 i$ R* D& etrue 10
7 ^  ]. f6 F; a1 v! Bfalse0 t7 p& A. Q* c# R+ \
true; k) Y1 x" i! y* E$ r6 f& e
true
作者: lilianjie    时间: 2012-1-11 18:09
C<c>:=ComplexField(5);C;
% `5 K2 n  z. X# @, ^7 B4 c+ C3 a; Y
3 i3 ?; ?# G2 \n:=C!(11+22*c);! J" W; @) p7 C! b7 t
ComplexConjugate(n);
0 a# p  c. |# W) z1 @3 p8 qComplexConjugate(n) ;复共轭
& e5 j, M, I" b3 h: YConjugate(n) ;共轭
  t/ a! ~  @- ~/ S9 \N:=Norm(n) ;N;范
A:=AbsoluteValue(n);A;
1 a+ o9 S2 `6 x2 J, kAbs(n) ;
" u4 Q0 h: A& D5 `. U8 ~- u! yA^2 eq N;" C9 e& h# X( U% y, u( e: u3 P+ I
+ G2 @$ e: F" p
Complex field of precision 56 V  _- R+ E* M' g7 R) e
11.000 - 22.000*c
- I  l% B7 c" T8 o  j11.000 - 22.000*c" o' o4 O7 R' ~  f1 d" Y! w% n3 p; v" I
11.000 - 22.000*c
# v7 F" w% @' g) ~1 g605.00" v. }% n* r# A
24.597' x/ j! _( T. P, V+ \  ]
24.5970 \. S, {8 R9 M% z) n2 m
true' W" i4 [# ]- J4 G: h. I  p8 L

作者: lilianjie    时间: 2012-1-11 18:32
本帖最后由 lilianjie 于 2012-1-11 18:49 编辑
! X! k' V/ }4 g0 h9 O
5 G% D& n4 @! T, cC:=RealField();C;
1 Q# C( _5 z+ u- c
' K* |4 ~+ j3 W0 a) u* _n:=C!1234/567;- L* a% C, x: h% K9 q( @0 {7 p+ m
ComplexConjugate(n);5 E6 g! `  @' r3 k8 [
ComplexConjugate(n) ;
( \% @8 B% S% v& u1 Q  G1 z. R9 I, g, |: ]; Q2 _1 n/ G. ~* G
N:=Norm(n) ;N;范% v8 x0 j" G8 X2 ^7 W0 t
A:=AbsoluteValue(n);A;范
: }: V* U8 }( H' DAbs(n) ;
7 h' n$ z! l" k" p9 u7 L$ r- z- i* vA eq N;
. i" z: H; _' o5 X" x! |Ilog2(n);
( ^; V: B- j  `- {! [' ]Truncate(n);取整
) C# F" k0 f/ o3 T0 QRound(n);
% \: A/ V3 R& L( [! d: QFloor(n);下限整数
( n- I: `4 I  A- ^Ceiling(n);上限整数+ K: B6 z8 X( B) y. w6 m+ Q
Sign(n);取正/负/0

; H6 n  u0 c( {0 x$ U$ @+ }Real field of precision 30# F* y; g1 _3 p3 W
2.17636684303350970017636684303
1 T6 @8 ^0 X7 d3 a4 G6 C2.176366843033509700176366843038 S% Y  {4 ^7 m. d* A, v7 m: r
2.17636684303350970017636684303
3 I9 V% G; N2 z2.17636684303350970017636684303
) n" T) h9 g/ A* V8 f! {2.17636684303350970017636684303
2 X* F  H8 \4 ^1 d$ J8 UTURE+ M2 c; n2 ~0 B% ^
1
. F9 q# j% k1 M1 N0 y& ^2) A9 O& M2 L# z5 J' ^. |0 h" j
2; N3 {7 r/ C+ t, o4 E  J6 n
2
# T% }0 u9 @, |7 @4 p3/ {% G7 k; ~$ {5 W& ~7 y) A- a
13 t" _3 z1 S9 j# O# f

  t8 Y9 y+ g6 ^* q3 D% O& t" h; U
8 ~  V) r& [1 K/ Z* pInfinity() ;∞' m7 i7 I$ y$ v# x. w4 h; y
MinusInfinity() ;-∞

& W9 E/ m+ [$ `' I' S# `Infinity() ;' G/ d9 [" Q2 y' M$ |
MinusInfinity() ;6 R( Z' `; I" E3 `
NextPrime(n) ;# [6 _% s3 X0 O5 l
PreviousPrime(n);
2 P' P" D8 `% ?" FNthPrime(n) ;* Z# o; i: k7 ]. v" v. |+ `& Z- Z

作者: lilianjie    时间: 2012-1-11 18:57
n:=123456;! d' x3 S6 K. m# Z3 [8 d
* O2 f' q! o8 J6 L% p
CarmichaelLambda(n+1000);卡米歇函数
5 O/ P6 b! `  S0 M; hFactoredCarmichaelLambda(n);分解卡米歇函数$ L/ R: l, ?+ F# @
; y) ~; M2 \# a0 j% x; W# Y
DivisorSigma(4, n) ;因子和
8 L. s- x  R, ]DivisorSigma(3, n) ;% S8 ]# B; m$ L  `- O$ b
SumOfDivisors(n);因子和函数
8 x1 |2 w7 I# Z' I! v: uNumberOfDivisors(n);因子函数
4 S" I# I* Q/ m; p- t' k1 _3 F9 Q  ^2 j8 i
EulerPhi(n) ;欧拉数, }& f  O7 O# [' u
FactoredEulerPhi(n);分解欧拉数
# ^/ z) K, D! V) }2 }EulerPhiInverse(n);欧拉函数求逆4 \" Y" {9 i0 q
FactoredEulerPhiInverse(n) ;欧拉函数求逆分解
MoebiusMu(n) ;墨氏函数LegendreSymbol(3, 103);勒让德符号; O# h- H* G7 ~: i) v. E
LegendreSymbol(-3, 103);
+ d4 N7 H9 e. _% QLegendreSymbol(34, 103);
; I9 v8 a3 k# n( c7 K& rLegendreSymbol(12, 13);! I6 _5 Z, |2 l' {; d
LegendreSymbol(19, 3);7 Q9 d0 a7 c7 N  w
JacobiSymbol(15, 13);雅可比符号9 t  J+ {- T) l/ ]7 h* V
JacobiSymbol(15, 13);
1 ?+ {* ~/ s& w7 SJacobiSymbol(150, 103);
" k; n% @8 S! c+ w7 ?1 y8 uKroneckerSymbol(2, 11);7 p* q6 C$ ?; s" [
KroneckerSymbol(5, 13);克罗内克符号
* o4 u" k5 @7 l; }& F" A% I5 AKroneckerSymbol(77, 4);

  M" [* A: I: E8 M
; i3 H/ O8 O) o, w7590
7 b: H( L% X$ e/ D4 C- s" F& {[ <2, 4>, <3, 1>, <107, 1> ]
3 U, B8 ?+ Q. {0 N# t250845527446699736708
) u4 R6 |. [9 S3 i9 H2230091146719632
' ^3 }. ~% w# u327152
0 y# t4 V/ a) X' j! v, s3 t7 w28
6 s) }1 C% |7 R  y41088
4 m0 l) m: q; s( E) G[ <2, 7>, <3, 1>, <107, 1> ]0 J) [) [" F& e6 r* N- o
[ 123457, 131189, 133757, 185187, 216069, 246914, 246916, 246944, 262378, : A' \5 I, U" U, J
267514, 288092, 308680, 370374, 370404, 370416, 432138, 463020 ]
3 k4 V5 P; w3 x" J9 \[' A5 Z: K/ I% k& i. ?5 l3 M
    [ <2, 1>, <3, 1>, <7, 1>, <10289, 1> ],3 @  S( [6 ~; j# s% C
    [ <2, 1>, <3, 1>, <61729, 1> ],( C! s0 l3 v* P# U# d7 m6 i+ u, B
    [ <2, 1>, <13, 1>, <10289, 1> ],7 {$ L) c! }, X2 W" ?0 |) Z
    [ <2, 1>, <17, 1>, <7717, 1> ],- s- }( m6 Q; j) O% g; o4 L
    [ <2, 1>, <123457, 1> ],* r! Y4 C! C/ Y* p$ o& ~2 |
    [ <2, 2>, <3, 1>, <5, 1>, <7717, 1> ],% _) {3 ]) j/ }
    [ <2, 2>, <3, 2>, <10289, 1> ],. V" i$ b; x- t
    [ <2, 2>, <7, 1>, <10289, 1> ],
. O* Y  F  W4 h1 A4 k1 O/ f# z    [ <2, 2>, <61729, 1> ],2 s0 s. _# L. d- ]1 c) ?  Q+ [
    [ <2, 3>, <5, 1>, <7717, 1> ],
9 n( g$ n' M. @$ r1 {# B$ v    [ <2, 4>, <3, 1>, <7717, 1> ],
3 d- K' Y. C4 Y4 M) O    [ <2, 5>, <7717, 1> ],% b& V2 V4 w3 l) C  e
    [ <3, 1>, <7, 1>, <10289, 1> ],
$ i8 H  `! @% k9 Z5 G& {    [ <3, 1>, <61729, 1> ],
0 S. x( B2 ?2 Y$ e) U: i    [ <13, 1>, <10289, 1> ],
% ~! d' I- [" [2 V( f    [ <17, 1>, <7717, 1> ],. E7 p! V$ o, G  l1 l
    [ <123457, 1> ]8 M" A1 A9 x' S
]& F3 E4 o; M0 w- h; {' v
0
) G$ G# r& `% [3 A; |. _- _& {( |-1
& g( `, z7 }# M9 B6 S# m15 r$ c5 y) @+ q) U# M
1
5 X( ?& [4 X" c2 r1; C* u, W" ]) D8 n0 }7 {3 K$ C
1
6 W$ u0 S5 B* D-1
3 Y' D) g; m  W) D0 h. W-18 W) G8 M' J5 H
-1
3 A5 r0 q7 \0 I3 e-1
8 V, h2 L$ V; I3 i" z$ J-1* j, u# r/ g3 B
1
作者: lilianjie1    时间: 2012-1-11 19:35
墨氏函数可表示欧拉函数因子和函数等算术函数

666.GIF (13.92 KB, 下载次数: 308)

666.GIF

555.GIF (8.27 KB, 下载次数: 292)

555.GIF

123.GIF (15.34 KB, 下载次数: 310)

123.GIF


作者: 孤寂冷逍遥    时间: 2012-1-12 15:06

作者: 苏惟嫣    时间: 2012-1-12 16:43

作者: 孤寂冷逍遥    时间: 2012-1-12 17:01





欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5