【Topcoder SRM 671 Div1 Level3】—BearDestroys（轮廓线dp）

传送门

考虑直接压列是不行的
只能压对角线
发现两条对角线的转移只用一直压$m+1$个状态
然后乱写写就可以了
讨论的情况还是不少

还因为各种傻叉问题调了好久

#include<bits/stdc++.h>
using namespace std;
const int RLEN=(1<<20)|5;
inline char gc(){
	static char ibuf[RLEN],*ib,*ob;
	(ib==ob)&&(ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
	return (ib==ob)?EOF:*ib++;
}
#define gc getchar
inline int read(){
	char ch=gc();
	int res=0,f=1;
	while(!isdigit(ch))f^=ch=='-',ch=gc();
	while(isdigit(ch))res=(res+(res<<2)<<1)+(ch^48),ch=gc();
	return f?res:-res;
}
#define pb push_back
#define re register
#define fi first
#define se second
#define pii pair<int,int>
#define cs const
#define bg begin
template<class tp>inline void chemx(tp &a,tp b){a<b?a=b:0;}
template<class tp>inline void chemn(tp &a,tp b){a>b?a=b:0;}
int mod;
inline int add(int a,int b){return (a+=b)>=mod?a-mod:a;}
inline void Add(int &a,int b){(a+=b)>=mod?a-=mod:0;}
inline int dec(int a,int b){return (a-=b)<0?a+mod:a;}
inline void Dec(int &a,int b){(a-=b)<0?a+=mod:0;}
inline int mul(int a,int b){return 1ll*a*b%mod;}
inline void Mul(int &a,int b){a=1ll*a*b%mod;}
inline int ksm(int a,int b,int res=1){for(;b;b>>=1,Mul(a,a))(b&1)&&(Mul(res,a),1);return res;}
inline int Inv(int x){return ksm(x,mod-2);}
cs int M=(1<<15)|5;
int f[2][M],g[2][M];
int n,m,lim,now;
class BearDestroys{
	public :
		inline int solve(){
			lim=1<<(m+1);
			now=0,f[0][0]=1;
			for(int i=1;i<=n+m-1;i++){
				now^=1;
				memset(f[now],0,sizeof(f[now]));
				memset(g[now],0,sizeof(g[now]));
				for(int s=0;s<lim/2;s++)Add(f[now][s<<1],f[now^1][s]),Add(g[now][s<<1],g[now^1][s]);
				for(int j=1;j<=m;j++){
					if(j>i)continue;now^=1;
					memset(f[now],0,sizeof(f[now]));
					memset(g[now],0,sizeof(g[now]));
					if(i-j+1>n){
						for(int s=0;s<lim;s++)
						Add(f[now][(s&(1<<j))?(s^(1<<j)):s],f[now^1][s]),
						Add(g[now][(s&(1<<j))?(s^(1<<j)):s],g[now^1][s]);
						continue;
					}
					for(int s=0;s<lim;s++)if(f[now^1][s]||g[now^1][s]){
						if(s&(1<<j))Add(f[now][s^(1<<j)],mul(2,f[now^1][s])),Add(g[now][s^(1<<j)],mul(2,g[now^1][s]));
						else{
							if(i-j+1==n){
								if(j==m)Add(f[now][s],mul(2,f[now^1][s])),Add(g[now][s],mul(2,g[now^1][s]));
								else{
									Add(f[now][s|(1<<j)],mul(f[now^1][s],2));
									Add(g[now][s|(1<<j)],mul(add(g[now^1][s],f[now^1][s]),2));
								}
								continue;
							}
							if(j==m){
								if(s&(1<<(j-1))){
									Add(f[now][(s&(1<<j))?(s-(1<<j)):s],mul(2,f[now^1][s]));
									Add(g[now][(s&(1<<j))?(s-(1<<j)):s],mul(2,g[now^1][s]));
								}
								else{
									Add(f[now][((s&(1<<j))?(s-(1<<j)):s)|(1<<(j-1))],mul(2,f[now^1][s]));
									Add(g[now][((s&(1<<j))?(s-(1<<j)):s)|(1<<(j-1))],mul(2,add(f[now^1][s],g[now^1][s])));
								}
								continue;
							}
							if(!(s&(1<<(j-1)))){
								Add(f[now][s|(1<<(j-1))],f[now^1][s]);
								Add(g[now][s|(1<<(j-1))],add(f[now^1][s],g[now^1][s]));
							}
							Add(f[now][s|(1<<j)],mul(f[now^1][s],1+(!!(s&(1<<(j-1))))));
							Add(g[now][s|(1<<j)],mul(add(g[now^1][s],f[now^1][s]),1+(!!(s&(1<<(j-1))))));
						}
					}
				}
			}
			int res=0;
			for(int i=0;i<lim;i++)Add(res,g[now][i]);
			return g[now][0];
		}
		inline int sumUp(int _n,int _m,int _mod){
			n=_n,m=_m,mod=_mod;
			return solve();
		}
};
BearDestroys stargazer;
int main(){
	#ifdef Stargazer
	freopen("lx.cpp","r",stdin);
	#endif
	int n=read(),m=read(),mod=read();
	cout<<stargazer.sumUp(n,m,mod);
}