[BZOJ4475][JSOI2015]子集选取（DP+结论）

把集合看成一个$n$位二进制数，第$i$位为$1$表示集合中包含$i$，否则第$i$位为$0$表示集合中不包含$i$。可以想到，如果把二进制数的每一位拆开考虑，那么变成如下模型：
在一个边长为$k$的三角形内，共$\frac{k(k+1)}{2}$个位置，对于一个位置$(i,j)$，如果$(i-1,j)$和$(i,j-1)$都填充了$1$（如果不存在位置$(i-1,j)$或$(i,j-1)$则当作填充$1$处理），那么$(i,j)$可填充$0$或$1$，否则$(i,j)$只能填充$0$。求填充这个三角形的方案数。
可以得到，最后结果等于上面问题结果的$n$次幂。
对于这个问题，考虑一个DP的做法，设$f[i]$为到了第$i$行的方案数。
转移的方法是枚举$j$从$0$到$i-1$，即枚举第$i$行填充$1$的个数，最后加上$1$（前$i$行全部填充$1$）
可以得到，第$i$行中被填充为$1$的位置是此行的前$j$个位置。如果第$i$行的前$j$个值已经被确定为$1$，那么后$i-j$个值也被确定为$0$。以此类推，前$i$行的前$j$列也被定为$1$，这样只剩下一个边长为$i-j-1$的三角形没有被填充，并且这个三角形填充的数字不受已填充的数字的影响。所以得到转移：
$f[i]=1+\sum_{j=0}^{i-1}f[i-j-1]$
也就是$f[i]=1+\sum_{j=0}^{i-1}f[j]$。
由于$f[0]=1$，因此归纳得出$f[k]=2^k$。
所以答案$=2^{nk}$。
代码：

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int PYZ = 1e9 + 7;
int n, K;
int qpow(int a, int b) {
    int res = 1;
    while (b) {
        if (b & 1) res = 1ll * res * a % PYZ;
        a = 1ll * a * a % PYZ;
        b >>= 1;
    }
    return res;
}
int main() {
    cin >> n >> K;
    cout << qpow(2, 1ll * n * K % (PYZ - 1)) << endl;
    return 0;
}