組合問題的算法實現

/*
問題描述：對於一組各不相同的數字，從中任意抽取1-n個數字，構成一個新的集合。
求出所有的可能的集合。例如，對於集合{1,2,3}，其所有子集爲{1}，{2},{3},{1,2},{1,3},{2,3}{1,2,3}， 
給定一個數組（元素各不相同），求出數組的元素的所有非空組合（即數組的所有非空子集）
解法：位向量法。用一個輔助數組表示各個元素的狀態。1表示在集合中，0表示不在數組中。遞歸地求解所有的子集。
算法描述如下：//這裏的算法對空集也輸出了，可修改之使得只輸出非空集合。
*/

#include <iostream>
using namespace std;
void getSubset(int list[],bool v[],int a,int b){
	if(a == b){
		for(int i = 0; i < b; i++){
			if(v[i])
			cout<<list[i]<<"  ";
		}
		cout<<endl;
		return;
	}
	v[a] = true;
	getSubset(list,v,a+1,b);
	v[a] = false;
	getSubset(list,v,a+1,b);

}
int main(){
    int li[] = {1,2,3,4};
    bool v[] = { false,false,false,false};
    getSubset(li,v,0,4);
}

/*
USER_ID: icpc#2017pre2015140537
PROBLEM: 1471
SUBMISSION_TIME: 2017-04-02 14:37:55

#include <iostream>
#include <vector>
#define MOD 1000000007
using namespace std;
vector<vector<long long> > muliMatrix(const vector<vector<long long> > &m1, const vector<vector<long long> > &m2) {
    vector<vector<long long> > res(m1.size() ,vector<long long>(m2[0].size(), 0));
    for (int i = 0; i < m1.size(); i++) {
        for (int j = 0; j < m2[0].size(); j++) {
            for (int k = 0; k < m2.size(); k++) {
                res[i][j] += (m1[i][k] * m2[k][j]);
            }
            res[i][j] = res[i][j] % MOD;
        }
    }
    return res;
}
     
vector<vector<long long> > matrixPower(const vector<vector<long long> > &m, long long p) {
    vector<vector<long long> > res(2 ,vector<long long>(2, 0));
     
    for (int i = 0; i < res.size(); i++) res[i][i] = 1;
    vector<vector<long long> > tmp(m);
    for(; p != 0; p >>= 1) {
        if ((p & 1) != 0) res = muliMatrix(res, tmp);
        tmp = muliMatrix(tmp, tmp);
    }
    return res;
}
int main()
{
    long long n;
    cin >> n;
    if (n == 0) cout << 1 << endl;
    else {
        vector<vector<long long> > v(2, vector<long long>(2, 0));
        v[0][0] = 2;
        v[0][1] = 1;
        v[1][0] = 3;
        v[1][1] = 0;   
        vector<vector<long long> > res = matrixPower(v, n - 1);
        cout << ((res[0][0] + res[1][0]) % MOD)<< endl;
    }
    return 0;
}

*/