0

我在 Maple 中有以下用于二次筛算法的代码：

  QS := proc(n::posint, (c := 1.2, cutoff := 7))
local st, nfb, svprimes, svplogs, nsvprimes, M, sievearray, polroots,
      kn, b, fb, rels, A, deps;
st := Initialization(n,c,cutoff);
print('Using multiplier : ');
print(st[1]);
print('Using smoothness bound: ');
print(st[3]);
fb := st[6];
print('Size of factor base: ');
print(ArrayNumElems(st[6]));
svprimes := st[7];
svplogs := st[8];
nsvprimes := st[9];
M := st[5];
sievearray := st[11];
polroots := st[10];
print('Sieving interval of length : ');
print(2*M+1);
sieve(svprimes, svplogs, nsvprimes, M, sievearray, polroots);
print('Sieving done, searching for smooth values.. ');
kn := st[2];
b := st[4];
rels := FindRelations[kn, b, M, fb, sievearray];
print(LinearAlgebra:-RowDimension(A));
gc();
print('Solving a matrix of size ');
print(LinearAlgebra:-Dimension(A));
deps := Dependencies(A);
print('Number of linear dependien fopund: ');
print(nops(deps));
print('Factors: ');
FindFactors(n, rels, deps)
and proc:

它给出以下错误：错误，“ (”意外。有人可以帮我解决这个问题吗？另外我想知道是否可以将此代码转换为 Java 或 C 代码？我不是要求任何人这样做，而是想知道这是否可能。我会自己试一试，但如果有人愿意，欢迎您。这是一个二次筛算法，我试图用它来分解超过一百位数的长数。

4

1 回答 1

0

改变，

QS := proc(n::posint, (c := 1.2, cutoff := 7))

至，

QS := proc(n::posint, c := 1.2, cutoff := 7)

如果您打算将以下内容作为函数调用，请更改，

FindRelations[kn, b, M, fb, sievearray];

至，

FindRelations(kn, b, M, fb, sievearray);

改变，

and proc:

至，

end proc:

将所有这些打印语句中的所有单右引号（uneval 引号，到 Maple）更改为单左引号（名称引号）或双引号。更好的是，使用userinfo语句。

尝试gc()像那样可靠地调用不一定会按您的预期工作。

于 2013-11-05T04:28:12.320 回答