CF-300D-Painting Square（dp or NTT）

题意：http://codeforces.com/problemset/problem/300/D

四分树 分k次的所有种类数，层数有限制

思路：

可以根据二叉树推广，不停往上加点的思想

（https://www.cnblogs.com/CJLHY/p/11099739.html）

首先讲dp：

从方程看出 i-j 和 k-l 的地位其实是等价的，所以我们可以先搞 i+j 的所有情况，也就是 k+l的所有情况，记为temp【sum】

而我们要的答案也就是 ans【i+j】== temp【i】【j】（就是叠两次）

而NTT方法就是加快卷积https://www.luogu.com.cn/problemnew/solution/CF300D

#define IOS ios_base::sync_with_stdio(0); cin.tie(0);
#include <cstdio>//sprintf islower isupper
#include <cstdlib>//malloc  exit strcat itoa system("cls")
#include <iostream>//pair
#include <fstream>//freopen("C:\Users\13606\Desktop\Input.txt","r",stdin);
#include <bitset>
//#include <map>
//#include<unordered_map>
#include <vector>
#include <stack>
#include <set>
#include <string.h>//strstr substr strcat
#include <string>
#include <time.h>// srand(((unsigned)time(NULL))); Seed n=rand()%10 - 0~9;
#include <cmath>
#include <deque>
#include <queue>//priority_queue<int, vector<int>, greater<int> > q;//less
#include <vector>//emplace_back
//#include <math.h>
#include <cassert>
#include <iomanip>
//#include <windows.h>//reverse(a,a+len);// ~ ! ~ ! floor
#include <algorithm>//sort + unique : sz=unique(b+1,b+n+1)-(b+1);+nth_element(first, nth, last, compare)
using namespace std;//next_permutation(a+1,a+1+n);//prev_permutation
//******************
clock_t __START,__END;
double __TOTALTIME;
void _MS(){__START=clock();}
void _ME(){__END=clock();__TOTALTIME=(double)(__END-__START)/CLOCKS_PER_SEC;cout<<"Time: "<<__TOTALTIME<<" s"<<endl;}
//***********************
#define rint register int
#define fo(a,b,c) for(rint a=b;a<=c;++a)
#define fr(a,b,c) for(rint a=b;a>=c;--a)
#define mem(a,b) memset(a,b,sizeof(a))
#define pr printf
#define sc scanf
#define ls rt<<1
#define rs rt<<1|1
typedef pair<int,int> PII;
typedef vector<int> VI;
typedef unsigned long long ull;
typedef long long ll;
typedef double db;
const db E=2.718281828;
const db PI=acos(-1.0);
const ll INF=(1LL<<60);
const int inf=(1<<30);
const db ESP=1e-9;
const int mod=7340033;
const int N=(int)1e3+10;

ll dp[50][N],temp[50][N];

int main()
{
    for(int i=0;i<=31;++i)dp[i][0]=temp[i][0]=1;
    for(int c=1;c<=31;++c)
    {
        for(int i=1;i<=1000;++i)
            for(int j=0;j<i;++j)
                dp[c][i]=(dp[c][i]+temp[c-1][j]*temp[c-1][i-1-j]%mod)%mod;
        for(int i=1;i<=1000;++i)
            for(int j=0;j<=i;++j)
                temp[c][i]=(temp[c][i]+dp[c][j]*dp[c][i-j]%mod)%mod;
    }
    int q;
    sc("%d",&q);
    while(q--)
    {
        int n,k;
        sc("%d%d",&n,&k);
        int cnt=0;
        while(n&1&&n>1)
        {
            cnt++;
            n/=2;
        }
        pr("%lld
",dp[cnt][k]);
    }
    return 0;
}

/**************************************************************************************/