TopCoder SRM 689 Div1 500 MultiplicationTable3

get到了一个新套路

我们先把空集也算进来，直接x++
假设我们现在构造出了一对(n,x)
假设此时的矩阵为
aaa
aaa
aaa

那么可以这么构造出(n+1,x+1)，因为我们发现集合中只要选了n+1，其他所有都得选，所以只多了一个１到n+1的全集。
aaa1
aaa2
aaa3
1230

还可以这么构造出(n+1,x*2)，因为对于之前的每个集合，都可以选择放或不放n+1
aaa3
aaa3
aaa3
3333

然后就把x拆成二进制，随便构造一下就行了。

#include <bits/stdc++.h>
#define ll long long
#define fr(i,x,y) for(int i=x;i<=y;i++)
#define rf(i,x,y) for(int i=x;i>=y;i--)
using namespace std;
const int N=21;
int n,a[N][N];
int m,d[N];

template<class T> void checkmin(T &a,const T &b) { if (b<a) a=b; } 
template<class T> void checkmax(T &a,const T &b) { if (b>a) a=b; }

class MultiplicationTable3 {
public:
    vector <int> construct( int x ) ;
};

void mul2(){
	n++;
	fr(i,1,n) a[i][n]=a[n][i]=i;
}

void plus1(){
	n++;
	for(int i=1;i<n;i++) a[i][n]=a[n][i]=i+1;
	a[n][n]=1;
}

vector <int> MultiplicationTable3::construct(int x) {
    x++;
    n=0;
    while(x){
    	d[++m]=x&1;
    	x>>=1;
    }
    rf(i,m-1,1){
    	mul2();
    	if (d[i]) plus1();
    }
    vector<int> ans;
    fr(i,1,n) fr(j,1,n) ans.push_back(a[i][j]-1);
    return ans;
}