Skip to main content
9 / 12
タイトルを編集

ある条件を満たす、matrixクラスのオブジェクトを生成したい(C++)

C++の初心者です。以下の通り、質問します。わかる方、いらっしゃれば、ご教示ください。

聞きたいこと

C++でmatrixクラスのオブジェクトで、以下の条件を満たすものを作成したいです。

  • オブジェクトの列数はSである
  • オブジェクトは0か1を成分に持つ
  • オブジェクトのそれぞれの行の要素の和はNである

ここで、S個の成分のうち、N個、1を選ぶ組み合わせは、M = S C N (Cはコンビネーション)通り、存在します。したがって、

  • オブジェクトの行数はMである

これら4つの条件を満たすオブジェクトを作成したいです。

具体例を挙げます。例えば、S=4、N=2であれば、オブジェクトの列要素(S=)4つのうち、(N=)2つのみを1にすることが許されます。ですので、取りうる場合の数(オブジェクトの行数)は(M=)6通りです。

上の条件を満たすオブジェクトは、それぞれの行が
(0,0,1,1), (0,1,0,1), (0,1,1,0),
(1,0,0,1), (1,0,1,0), (1,1,0,0)
となることに注意すると、

(0,0,1,1
0,1,0,1
0,1,1,0
1,0,0,1
1,0,1,0
1,1,0,0)

というM(=6)×S(=4)のオブジェクトになっていることがわかります。
このようなオブジェクトを任意の自然数S, Nに対して作成したいというのが、やりたいことです。

試したこと

N=2の場合には、以下のようなコードを用いることで、一般のSに対して、オブジェクトを作成することができました。

#include <vector>
#include <complex>

using namespace std;
typedef complex<double> Complex;
typedef vector<Complex> cvector;
typedef vector<cvector> cmatrix;

int main(void)
{
    int S = 4;
    int N = 2;

    //Mを求める
    int a1 = 1; // S*(S-1)*・・・( S - (N-1)) //N回の掛け算
    int a2 = 1; // N!



    for (int i = 0; i < N; i++) {
        a1 = a1 * (S - i);
    }

    for (int i = 0; i < N; i++) {
        a2 = a2 * (N - i);
    }

    int M = a1 / a2;  //行列の行数を決定 M = S C N
    
    //行列Bに値を代入
    cmatrix B(M, cvector((S), 0));
    for (int i = 0; i < S - 1; i++) {
        for (int j = 0; j < S - i - 1; j++) {
            B[j + S * i - 0.5 * i * (i + 1)][(S - 1) - i] = 1;
            B[j + S * i - 0.5 * i * (i + 1)][(S - 1) - (i + j + 1)] = 1;
        }
    }

    return 0;
}

これにより、cmatrixのオブジェクト"B"が

B =
(0,0,1,1
0,1,0,1
0,1,1,0
1,0,0,1
1,0,1,0
1,1,0,0)

という形で作成されました。
上のコードでは、S = 4としていますが、任意のSでこのコードは動きます。ただ、Nに対してはN=2しか許されません。(というのも、for文を使って、Bに直接、N回値を書き込むという操作を行っているからです。)
ネットで調べていると、ビット演算というものを使えばできると書いてあるものもありましたが、如何せん初心者であるため、よくわかっていません。

どなたか、ご教示をお願い致します。