問題タブ [factorization]

http://projecteuler.netから問題 3 を解決しようとしています。ただし、プログラムを実行しても何も出力されません。私は何を間違っていますか？問題: 数 600851475143 の最大の素因数は?

私が作成しようとしているプログラムは、ユーザー定義のInt値の数を入力し、それらの値の積を生成してから、その積のすべての要素を決定する必要があります。値は、ユーザーが負の数を入力するまでDo-Whileループを介して入力され、その後、すべての正の数が乗算されて、Doubleの形式の積が形成されます。製品の要素をどうやって把握できるかはかなり確信していますが、実際の製品を見つけることになると、私は完全に途方に暮れています。これが私がこれまでに書いたコードですが、製品を作成する試みの失敗を除いたものです。

これが私が理解できたすべてであり、負のIntが入力されるまでValueプロンプトをループします。製品を入手する限り、ループの内側と外側の両方でIf-Statementsを試しましたが、コードが間違っているか、ステートメントが間違っているか、間違ったメソッドを試している可能性があります。誰かにコードを書いてもらうつもりはありませんが、誰かが私を正しい方向に向けることができれば、とてもありがたいです。

不明な数のユーザー定義値から製品を決定するために、どのようなヒントを提案しますか？

このメソッドは int を取ります - しかし、私はエラーメッセージを受け取り続けます.誰かが理由を知っていますか?

構築クラスと、ユーザー入力から数値を読み取り、数値の素因数分解を吐き出すメインメソッドの 2 つのクラスを作成しています。コードは Java を使用しています。

例:
数値を入力: 150
5
5
3
2

ただし、私のプログラムでは、要因のリスト全体を取得しています。

例:
数値を入力: 150
150
75
50
25
5
3
1

素因数を得るためにこれをどのように変更しますか?

主な方法:

これが私のクラスです：

質問はタイトルにあります。ここで、自己生産という用語をより明確にします。自己積とは、数とその桁の積を意味します。例：

自己積(1234) = 1234*1*2*3*4 = 29616.

私は2つの方法を試しました。

強引な

誰のアイデアも、1 と N (範囲の上限) の間のすべての組み合わせをチェックすることです。数値の自己積が範囲内にあるかどうかを確認します。これは比較的低い数値には理想的ですが、範囲が 10^20 と大きくなる可能性があることを念頭に置くと、結果が出力されるまでに時間がかかるため、問題になりつつあります。

因数分解

別のアイデアは、範囲内の数値を因数分解することです。範囲が 60 000 ～ 70 000 の場合、62 688 を確認すると、62 688 は 2*2*2*2*2*3*653 となります。653 を数字にすることはできないことがわかったので、元の数字にする必要があります。次に、正しい答えを得るために 62 688 のすべての因数を組み合わせる必要があり、62 688 = 2612 * 2 * 6 * 1 * 2 であるため、これは 2612 の自己積であることが出力されるはずです。

どちらの状況でも、すべての組み合わせをチェックするという大きな問題に直面しています。

PS私は、数値がn桁の場合、それが何らかの数値の自己積である場合、その数値は少なくともn / 2桁でn以下になることを発見しました。これにより、チェックする必要がある数値のリストが少し小さくなりますが、問題は解決しません。

次元 6x6 の行列 M があり、ランク 1 です。M=A*B となるように、次元 6x1 (A とします) と 1x6 (B とします) の 2 つの行列に因数分解するにはどうすればよいでしょうか。

spptrfMATLAB が LAPACK関数をサポートしているかどうか知っていますか。

この関数は、巨大な正定対称行列のコレスキー分解を計算する必要がある場合に非常に便利です。

1 次元行列として保存されている上三角行列のみを入力として与えることで、因数分解を行うことができます。

または、chol組み込み関数は既にspptrf内部で使用されていますか?

編集

lapackFile Exchange http://www.mathworks.com/matlabcentral/fileexchange/16777-lapackでライブラリを見つけることができ、目的の関数が実装されていspptrfます。

編集2

を呼び出すたびに、自分のマシンで実行されている MATLAB が致命的にクラッシュしますspptrf。

この関数を直接処理する別の方法はありますか?

与えられた式 :

私が試した解決策：

すでに知っている用語 : A、B、C、D、K

p、q、r を指定して用語 s を検索します。

s が < 0 になる まで続行します。すべての項に対して p=0、q=0、r=1...n;

同様に、 s が < 0 になる まで続行します。p=0、q=1..n、r=1...n のすべての項に対して。

と、

最後に、 s が < 0 になる まで続行します。p=1..n、q=1..n、r=1..n のすべての項に対して。

p、q、r を更新するためのコード化された 3 つのループ。

K が 1000...、8145、45000 などのように大きくなると、計算に時間がかかります。

外部ライブラリを提案しないでください...コーディングソリューションを探しています。

スニペットの例

また、次のことに注意してください: preSpecifiedIterations -> 計算前に必要な反復回数を決定することは可能ですか?

上記の問題を解決するためのより良いアルゴリズムはありますか?

たくさん読んでくれてありがとう。

ポラードの P-1 因数分解を Python で実装しようとしています。Rhoメソッドにはいくつかの答えがありますが、このp-1は異なり、p-1についてここで提供できる最良のものはwikiとWolframです。

http://en.wikipedia.org/wiki/Pollardの s_p_%E2%88%92_1_algorithm

http://mathworld.wolfram.com/Pollardp-1FactorizationMethod.html

これは n から何かを因数分解していますが、一貫して p が見つかりません。np と sp はそれぞれ numpy と scipy からのものです。したがって、sp.uint64 の組み込み関数は unsigned long 64 int (予想される整数のサイズのため) であり、np.prod(p) はリスト p の累積積 pi です。

出力で p が見つかりません:

私はPythonを学んでいるので、単純な間違いかもしれません。power2() 関数は 2 乗によるべき乗を使用し、基本的に非常に大きな整数用のスーパーチャージされた pow() です。euc_al_i() は単なる gcd です。好きな gcd() を使用できますが、学習しているので、これらを自分で作成したかったのです。

ここで何がひどく間違っていたのかを突き止めようとしているので、比較的小さな n (20 ビット長ほど) からでも p が見つかりません。