[Sdoi2010]地精部落

/*
Orz 魏铭。
想出来满分做法。。但是只能70分。
首先，设状态f[i][H][0/1]表示序列posi是H，且为山峰/山谷。然后转移很显然。
不过实现只能强打记忆化搜索。为什么呢？因为他是一个排列。转递推就要状压。就是用bitset也会爆。
抛开空间不讲。时间也过不了。
为什么？转移时O(n)的。如何优化？用前缀和优化变为O(1)。
在处理好空间且保证状态正确的前提下，转成递推，用前缀和转移就好了。当时也想到了，但，写不出来。
当场 魏铭 就用一种很神奇的方法，解决了我的问题，然，我还是看不懂。码力太弱了。。。
此题对我感触很大。
思维：从裸暴力->状压剪枝->换状态->换状态->换状态->优化->正解思路≠>正解代码（码力弱）。
*/
#include<cstdio>
#include<cstring>
#include<algorithm>
#define set(x) freopen(#x".in","r",stdin);freopen(#x".out","w",stdout);
using namespace std;
typedef long long ll;
const int N=555;
int n,vis[N];
ll f[N][N][2];
ll ans,mod;
//d=1 high
//d=0 low
ll dfs(int cur,int h,bool d){
    ll &res=f[cur][h][d];
    if(~res) return res;
    res=0;
    if(cur>=n) return res=1; 
    if(d){
        for(int i=1;i<h;i++){
            if(!vis[i]){
                vis[i]=1;
                res+=dfs(cur+1,i,!d);
                vis[i]=0;
            }
        }
    }
    else{
        for(int i=h+1;i<=n;i++){
            if(!vis[i]){
                vis[i]=1;
                res+=dfs(cur+1,i,!d);
                vis[i]=0;
            }
        }
    }
    return res%=mod;
}
inline void DFS(){
    for(int i=1;i<=n;i++){
        vis[i]=1;
        dfs(1,i,0);
        dfs(1,i,1);
        vis[i]=0;
    }
    for(int i=1;i<=n;i++){
        ans=(ans+f[1][i][0]+f[1][i][1])%mod;
    }
}
int main(){
    set(goblin);
    memset(f,-1,sizeof f);
    scanf("%d%I64d",&n,&mod);
    if(n==1){printf("%I64d",1LL%mod);}
    if(mod==1){printf("0");}
    DFS();
    printf("%I64d",ans);
    return 0;
} 
/*
附魏铭100分代码
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
using namespace std;
#define zbwmqlw
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define abs(x) ((x)>0?(x):-(x))
#define pow(x) (1<<(x))
const int inf=0x3f3f3f3f;
typedef __int64 LL;
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;}

const int N=4300;
int a[20];
bool vis[20],now,last;
int n,ans;
LL f[2][N][2],mo;

void dfs(int depth)
{
    if (depth==n+1) {
        ++ans;
        return;
    }
    int i;
    for (i=1;i<=n;++i) if (!vis[i]) {
        a[depth]=i;
        if (depth>2 && (a[depth]>a[depth-1])==(a[depth-1]>a[depth-2])) continue;
        vis[i]=true;
        dfs(depth+1);
        vis[i]=false;
    }
}

void solvedfs()
{
    ans=0;
    dfs(1);
    printf("%d %d\n",n,ans);
    exit(0);
}

void solvedp()
{
    int i,j,k;
    now=true; last=false;
    memset(f[now],0,sizeof(f[now]));
    for (j=0;j<n;++j) f[now][j][0]=f[now][j][1]=1;
    for (i=2;i<=n;++i) {
        now^=1; last^=1;
        memset(f[now],0,sizeof(f[now]));
        LL s0=0,s1=0;
        for (k=0;k+i-1<=n;++k) s1+=f[last][k][1];
        for (j=0;j+i<=n;++j) {
            s0+=f[last][j][0];
            s1-=f[last][j][1];
            f[now][j][1]=s0%mo;
            f[now][j][0]=s1%mo;
        }
    }
    printf("%I64d\n",(f[now][0][0]+f[now][0][1])%mo);
}

int main()
{
    freopen("goblin.in","r",stdin);
    freopen("goblin.out","w",stdout);
    scanf("%d%I64d",&n,&mo);
    solvedp();
    fclose(stdout);
}

*/