数学建模社区-数学中国

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

作者: lilianjie    时间: 2012-1-11 17:58
标题: 一些初等算术逻缉函数
Z:=IntegerRing() ;- r8 }3 Z, }7 z
Z;
: i  f' h- L0 o' VZ12:=IntegerRing(12);
: X% D5 p; K! w# u$ \. fZ12;
7 _+ p7 c! N2 Z( O1 kZ17:=IntegerRing(19); # X1 c# }  y' e) a: U+ }
Z17;+ c9 y% g3 L0 B# E  l# q0 h
n:=Z!100;) b4 K( a9 j4 @0 M
n;
! ^( V! S' y% o$ o0 g6 j* C7 G7 hm:=Random(Z12);
8 X9 q1 T: [, H) ^  g' _m;7 |5 |; d- ~8 `; a. I# e5 E

# q3 P: n' P& c6 Q+ A* i0 e9 I- p& yn div m;    求余
1 J" G5 u7 P0 E; B) N& _. N& ZZ!n mod Z!m;" V+ t2 c7 q: x* j" V9 ?- p
ExactQuotient(n, 5);求商( k7 c8 f! B0 |# g. e! b
ShiftLeft(n, 1);乘22 t. x" S/ @, Z# V0 @
ShiftLeft(n, 3);乘8  ^/ F: u1 L2 L! r
ShiftLeft(n, 10);乘1024) J( I4 E% X" L9 C5 A+ G7 o0 c
ShiftRight(n, 2);除4* Y+ O4 D4 Y" J& H
ShiftRight(n, 3);除8
2 m5 k- F3 R& {/ Z" d7 l  sShiftRight(n, 10);除1024
! ]0 s- d6 _! n& ?2 QModByPowerOf2(n, 2);模43 @/ }* }' q5 I! `
ModByPowerOf2(n, 3);模8
) L# }6 `7 g- j2 rModByPowerOf2(n, 5);模32
IsEven(n) ;8 E  I" t$ a  `5 X/ y; r/ a% N' J5 y
IsOdd(n);
1 h- b* @2 d$ Z! \% ], c' hIsDivisibleBy(n, 13);
% v6 E, v! N, L) A- K5 L& GIsDivisibleBy(n, 2);' X$ w  t3 I/ a9 w1 ?
IsSquarefree(n) ;
: }3 E0 U; \* V, e( @" SIsPower(n);, h7 |- Z+ D) W8 j5 z: \9 m
IsPower(n, 10) ;
) i0 B1 _& a" q  E; ^IsPower(n, 3) ;
  M( u2 N- |% u# a7 \" JIsSquare(n);- L) M# ^* i8 ?! ^& y, D0 ]% L
IsPrime(n) ;
0 {: [  ~. u' s, [7 h- ^IsIntegral(n);& q8 f1 Z6 ^3 }4 k: r& C
IsSinglePrecision(n);单精度
* m+ J' f+ s: Q0 {( G2 ]; `' @* {, a9 K

: r6 ]! G( _5 q* S7 VInteger Ring
7 _$ e& {0 f- u7 P6 |Residue class ring of integers modulo 12" F/ O- q6 v# o( _' Q* [  r
Residue class ring of integers modulo 19
- ^& _  ]; o1 O  v100
) j( T* K, K; r6 @( {" \0 v9
: A! [9 h- ]0 |; Z$ O11
) x$ S  M. f  a+ }( D1  ~. ~6 n6 w5 p1 g% F5 q  W
20! @' b: d3 M+ e. x2 [- h5 Q4 u
200
- i9 x* h% S  h& T800
' o# r% g5 j7 v: w6 P& f102400
: c( b5 i& \0 m; L6 y# _* h# Z! ~25+ e4 l1 Y4 {6 w( c
128 `) ?) ]* r- e! N- N
0' O* L6 }. I' n1 o5 ?% i
0
$ m, d4 N% M( w5 y. r/ M% l4
- y% F1 D' F- w+ j8 z4
& [7 r7 q# \: e* atrue4 w- N7 E, g1 b
false
* E1 k0 X& {- X" {6 Z9 U7 cfalse# {: ?  @! c* B$ l
true* X4 s7 Q9 v& p( @
false
' l& b" S5 s4 S. F9 ?. o' Ctrue 10 2
3 E' \1 Y0 P# k& [% d- wfalse9 N/ i* v$ M' @, Q5 n  y; x
false
3 @8 {# B# d  N$ z- k! y/ Btrue 10- v; Q' i3 H' G/ a
false
+ |( `% o' w6 x) o$ [true
$ u4 d6 X9 x9 p* l! B) z1 b1 Ztrue
作者: lilianjie    时间: 2012-1-11 18:09
C<c>:=ComplexField(5);C;' C! A6 m/ c- ?1 x7 Y5 Q4 h
- S3 Z# {  V$ i5 x# R
n:=C!(11+22*c);
( e! z+ }; V) N3 ?1 Y' FComplexConjugate(n);( K8 K+ S: I4 A, x4 G3 a
ComplexConjugate(n) ;复共轭
- t6 k9 Q2 z* n, O* ?+ u! eConjugate(n) ;共轭% v! K% p5 w' V  P' N8 t9 S
N:=Norm(n) ;N;范
A:=AbsoluteValue(n);A;, U3 O$ \* S' I3 @" d  L+ r
Abs(n) ;* [6 Q- J6 ]8 I' P3 J5 Y9 G' F2 }
A^2 eq N;) B4 j, w: y) q0 P1 m) w% X6 C
9 W+ ?; c2 G  @% ^3 n# U
Complex field of precision 5
1 ~2 u* i- V) H( {. q; j3 y' q  Q11.000 - 22.000*c
7 X/ q# l! W% N7 C1 A# [4 I, h11.000 - 22.000*c1 T5 o7 j1 Q; z& D( Q( Y
11.000 - 22.000*c  c& j0 ?5 z* n9 A4 ]% G
605.00
6 \0 j, {" D) S* @1 f- f24.597
6 D) n/ d$ _' }5 g" z24.597: w8 _" i0 K" e) ?9 y0 w6 d
true" J" G9 Y  m0 F

作者: lilianjie    时间: 2012-1-11 18:32
本帖最后由 lilianjie 于 2012-1-11 18:49 编辑
$ T$ C7 P5 {( j( G. W  \3 t
; x- y& b5 u3 ~/ y  p3 }8 |5 fC:=RealField();C;1 o$ X: M6 j/ B, r
! _' w0 Z" A/ {( s6 g' i
n:=C!1234/567;' K' M, ?9 }: R8 m. e' {
ComplexConjugate(n);& J" C" g; p* S5 g, }( E8 m) u) P
ComplexConjugate(n) ;
- V* U& S$ S) e. j
( v; k' L7 l$ M8 {2 W8 S  {N:=Norm(n) ;N;范
' o* {2 |, |# yA:=AbsoluteValue(n);A;范
% g, z& S$ Y$ t( ZAbs(n) ;
6 n2 L$ x* O$ A6 T% a/ Q4 V& g$ IA eq N;
3 t& g2 d. U: |( H, Y7 K: p# sIlog2(n);
0 s! X6 I  P0 a# FTruncate(n);取整
( ]5 ]6 [: ^, s$ P: K6 aRound(n);
0 [* f$ B* g/ \! y  JFloor(n);下限整数/ p6 e+ r  c3 K* G) U% u% W
Ceiling(n);上限整数8 Y! I7 p& r2 x% l8 ?6 }
Sign(n);取正/负/0
0 I4 x; W/ v1 N0 `
Real field of precision 30% r: x; z5 f  ]$ I
2.17636684303350970017636684303
' ~. H6 `/ P4 f$ g/ g& M* J; P2.17636684303350970017636684303
4 J" C6 j1 ]: o# f2.17636684303350970017636684303
+ E5 i; x% ~" e' f2 ]  c9 \2.17636684303350970017636684303; B. g# n2 y# d; d- b7 t
2.17636684303350970017636684303
- t+ I) F# d7 D' H1 vTURE6 B& ?% [) X5 `5 O$ W
1
2 p; Q, |! v2 Y' L# y/ l3 r0 y2! @# O- x3 [- p" c8 k
2
; u( k; u6 m/ R' ]4 G6 T3 U) b2( V% T% D" p1 D- a) n. P
31 R' e) B+ ]3 n; I  x( |7 ~- Q
1$ x3 e9 p* Z  t: J  H
& t) V- B: T  ^+ `* e( O( i$ k

+ o* _$ n6 }0 o. GInfinity() ;∞
: V- O( _: w7 h; `; h1 I. VMinusInfinity() ;-∞

) @' o% _: @( W( e; \) xInfinity() ;# q' s3 i/ G  G5 V
MinusInfinity() ;
% `- Z9 V5 z" |NextPrime(n) ;
3 V  r( }  p& I8 Y9 LPreviousPrime(n);& K; F; y- y9 l% f; `: ~* i
NthPrime(n) ;
1 s6 o" Q4 z' q& X
作者: lilianjie    时间: 2012-1-11 18:57
n:=123456;& K* Q' }: b5 K7 F0 t8 y- J  y7 c

* ^; D/ d; T4 x( B3 JCarmichaelLambda(n+1000);卡米歇函数
- d: }8 b9 W. [  w9 k2 S$ bFactoredCarmichaelLambda(n);分解卡米歇函数
2 l+ u8 z( U; h3 L: p- T/ i) f3 p2 h$ e& @
DivisorSigma(4, n) ;因子和2 f) i( i5 z  \# i  F! {- l# U
DivisorSigma(3, n) ;0 f+ ?7 p1 {( p
SumOfDivisors(n);因子和函数
! ]1 L" v# D" m2 `5 b8 s! tNumberOfDivisors(n);因子函数9 W& v, V8 [8 B+ r, E- f3 D6 c% ~
+ n" R& s* X' i8 I
EulerPhi(n) ;欧拉数- S: {' t6 z/ `( U- e1 I# ]4 f9 P
FactoredEulerPhi(n);分解欧拉数, X8 N9 C2 A$ K9 @2 A/ C. y
EulerPhiInverse(n);欧拉函数求逆  S9 J0 V+ H5 f7 f9 D% ?0 O. |0 Y7 S
FactoredEulerPhiInverse(n) ;欧拉函数求逆分解
MoebiusMu(n) ;墨氏函数LegendreSymbol(3, 103);勒让德符号
& a9 U) z' f2 i: j! h( BLegendreSymbol(-3, 103);
2 U0 ^9 k, F, yLegendreSymbol(34, 103);
: k$ K# K3 Q: a; Q+ D  tLegendreSymbol(12, 13);
2 ], |+ G5 b/ h$ \LegendreSymbol(19, 3);3 n* E0 N9 L) u$ h& D0 `( a
JacobiSymbol(15, 13);雅可比符号
( X  z% }' S0 C; G* R9 BJacobiSymbol(15, 13);; |7 C4 J; v9 o( W+ R
JacobiSymbol(150, 103);
  ?1 ]/ E! }3 F: P4 r. T% |KroneckerSymbol(2, 11);
( L- @- X* h2 ~% b& L+ g0 A& qKroneckerSymbol(5, 13);克罗内克符号! d$ l# b1 x1 n! p9 |& l! x4 N
KroneckerSymbol(77, 4);

- [# m3 l. o) N  C+ ^6 g, y* j& G& g
7590
8 O& E( _9 L7 ][ <2, 4>, <3, 1>, <107, 1> ]
. B) X  v7 q. O9 ~1 D* a! V2508455274466997367086 H4 b* v( T0 \- I6 f. k
22300911467196327 k+ V" B. ?7 W. I0 z; q$ Q
327152
# G4 k8 D; O- p/ G6 D5 S28$ D' x% v; h& B$ A& A- l7 x
41088
# K% ]# D4 I& K[ <2, 7>, <3, 1>, <107, 1> ]
+ o  Q0 y( D+ u" N& x[ 123457, 131189, 133757, 185187, 216069, 246914, 246916, 246944, 262378, # H  H9 d) E% v7 u3 x
267514, 288092, 308680, 370374, 370404, 370416, 432138, 463020 ]
  ~% K( U: o: V! |8 v. Z5 A5 L[0 Z2 \. O# {1 [& p
    [ <2, 1>, <3, 1>, <7, 1>, <10289, 1> ],
' u/ j9 D6 h9 z, o/ W; J/ r    [ <2, 1>, <3, 1>, <61729, 1> ],1 O, \. D0 J! R9 ^5 M, C3 o; B/ g
    [ <2, 1>, <13, 1>, <10289, 1> ],7 e7 j0 Q; K1 B$ o# r. \0 [
    [ <2, 1>, <17, 1>, <7717, 1> ],
7 h* N* Z4 {# h; [0 I% v4 |    [ <2, 1>, <123457, 1> ],. e% ^; `( J# k" e% i
    [ <2, 2>, <3, 1>, <5, 1>, <7717, 1> ],
% i; d% Q) V1 s, b    [ <2, 2>, <3, 2>, <10289, 1> ],
/ r, W3 V2 Z, Y/ T! }    [ <2, 2>, <7, 1>, <10289, 1> ],
. K: d: D2 O* q3 S- D    [ <2, 2>, <61729, 1> ],
" D. ?' m: w4 ?- S  P! A! R: f  |' R8 k    [ <2, 3>, <5, 1>, <7717, 1> ],* d( J* {* p; `$ l: o8 p+ n
    [ <2, 4>, <3, 1>, <7717, 1> ],
/ b5 H9 `8 @8 A6 a  `    [ <2, 5>, <7717, 1> ],, h6 k0 h1 @' z) L" ^% ^7 q
    [ <3, 1>, <7, 1>, <10289, 1> ],! ^- S* ]8 P9 x& f
    [ <3, 1>, <61729, 1> ],( j  `3 J7 w! E9 A9 g) V% G' Q& w
    [ <13, 1>, <10289, 1> ],6 U. e' Z  e' b8 l6 b& I* O
    [ <17, 1>, <7717, 1> ],
& [4 P4 Z( ]+ R5 C    [ <123457, 1> ]
: g! [, h# [0 w8 J' b6 r1 k* X]  I9 G' A$ e4 _1 N4 ^& ^+ D5 z
0$ S7 f; }/ F' F- s" Y0 W
-12 B, j4 |9 I9 v- m& c8 ^
1/ x2 \  l# [$ V. R: m. [# [( K
1
  D, p7 c% p  s4 P: a6 D% C1
6 m- n* n- W' e9 Z7 [$ f1
0 G  `2 _4 E: C-17 u5 x5 U% a0 r2 V8 q. Z( l
-16 @1 g; k; J; K  Z7 T! X
-1
% b% ]; e, V% K0 b; J* [, g-1! q9 n& {* ?" q# p* h
-1  ]9 i) ^4 i: H( t7 f$ a
1
作者: lilianjie1    时间: 2012-1-11 19:35
墨氏函数可表示欧拉函数因子和函数等算术函数

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

666.GIF

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

555.GIF

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

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