摘要：
很久没更新东西，这个学期比较忙~~~基于最近对大作业的探讨，写了一些比较有用的代码，拿出来分享一下。直接进主题：SolveDiotN[A_List, M_] :=Module[{T = {}, Y = {}, X, x, y, i, d, n, u}, n = Length[A]; Array[d, n]; X = Array[x, n]; For[i = 1, i <= n, i++, AppendTo[Y, x[i] >= 0]]; For[i = 1, i <= n, i++, AppendTo[T, x[i]]]; d[2] = GCD[A[[1]], A[[2]]] 阅读全文
摘要：
终于考完Mathematica，趁有时间发个BLOG先。一直自学着潘承洞的 《初等数论》 \，里面关于原根、指标的计算实在繁杂，所以编了几个Packages在做题时用，希望对学习数论的朋友有用。（已经假定你具备已有的知识） 先介绍几个Mathematica的本身的命令： EulerPhi[n] // 欧拉函数MultiplicativeOrder[k, n, 1] //求mod n，k的指数 PrimitiveRoot[n] // 求n最小正原根 PrimePowerQ[n] // 判别n是否素数的方幂 MultiplicativeOrder[k, n, r] //(字符编辑原因，不能识别La 阅读全文
摘要：
实践了一个晚上，终于把大概的思路弄明白~~~~为此还编了个模板。在Mathematica 7.0模板如下： BeginPackage["自定义函数名`",{"需要用到的PACKAGE（通常在AddOns可以找到）`"," `"}] 自定义函数名::usage="对函数的注释"Begin["`Private`"] 程序段 End[]EndPackage[]想了想，比较空乏，给个原型吧BeginPackage["PrimitiveQ`", "NumberTheory`& 阅读全文