[2017纪中10-22]友谊 DP

60pts考场做法：
设计状态：f[i][j][k][l]为填到第i个友谊颗粒，浪费了j个奇数的友谊颗粒没有匹配，还剩k个偶数圆粒未匹配，还剩l个偶数皱粒未匹配。我们发现k+l=j-(i&1)，于是只需f[i][j][k]即可。
若当前在填奇数：
f[i][j][k]=f[i-1][j][k+1]+f[i-1][j-1][0]* (k==0)+f[i-1][j][k]+f[i-1][j-1][j-1]* (k==j-1)
就是分这个奇数颗粒是圆还是皱，有没有被浪费掉共四种情况。
若当前在填偶数:
f[i][j][k]=f[i-1][j][k-1]* (k!=0)+f[i-1][j][k]* (k!=j)
分圆皱讨论，注意转移要合法。
代码：

#include<iostream>
#include<cstdio>
#define ll long long 
using namespace std;
int n,m,mod;
ll f[2][3100][3100];
int main()
{
    scanf("%d%d%d",&n,&m,&mod);
    f[0][0][0]=1;
    for(int i=1;i<=2*m;i++)
    {
        int v=(i&1);
        for(int j=v;j<=((i+1)>>1);j++)
            for(int k=0;k<=min(j-v,(i>>1));k++)
                if(v) f[v][j][k]=(f[v^1][j][k+1]+f[v^1][j-1][0]*(k==0)+f[v^1][j][k]+f[v^1][j-1][j-1]*(k==j-1))%mod;
                else f[v][j][k]=(((k==0)?0:f[v^1][j][k-1])+f[v^1][j][k]*(k!=j))%mod;
    }
    ll ans=0;
    for(int j=0;j<=n;j++)
        for(int k=0;k<=m;k++)
            ans=(ans+f[0][j][k])%mod;
    printf("%lld",ans);     
    return 0;
}

题解60pts做法：
跟我的本质上应该是相同的。
首先发现第一位和最后一位填什么都不影响匹配，从2转移到n-1，最后答案*4即可。
f[p][i][j]表示用了p个奇数，有i个偶数圆粒未匹配，j个奇数圆粒未匹配，两位两位（偶奇）转移。
圆圆：f[p-1][i][j]
皱皱：f[p-1][i][j]
皱圆：f[p-1][i+1][j-1]，当i==0还可以f[p-1][0][j-1]，可能奇数位的圆粒浪费掉了。
圆皱：f[p-1][i-1][j+1]，当j==0还可以f[p-1][i-1][0]，可能奇数位的皱粒浪费掉了。
代码：

#include<iostream>
#include<cstdio>
#define ll long long
using namespace std;
int n,m,p;
ll f[310][310][310];
int main()
{
    scanf("%d%d%d",&n,&m,&p);
    f[1][0][0]=2;
    for(int i=2;i<=m;i++)
        for(int j=0;j<i;j++)
            for(int k=0;k<i-j;k++)
            {
                f[i][j][k]=f[i-1][j][k]*2+((j==0)?0:f[i-1][j-1][k+1])+((k==0)?0:f[i-1][j+1][k-1]);
                if(j==0&&k!=0) f[i][j][k]+=f[i-1][0][k-1];
                if(j!=0&&k==0) f[i][j][k]+=f[i-1][j-1][0];
                f[i][j][k]%=p;
            }
    ll ans=0;
    for(int i=0;i<n;i++)
        for(int j=0;j<n-i;j++)
            ans=(ans+f[m][i][j]*2)%p;
    printf("%lld",ans);                 
    return 0;
}

题解100pts做法：
因为题解只问了匹配完长度不超过2n的答案和，并不用求出每种长度的答案，我们还可以优化。
我们考虑一开始就先添加n个未匹配的偶数点，这样到最后奇数点应该全部能匹配完才合法（不能存在浪费的情况）。假设一开始我们人为添加的是j个圆粒，n-j个皱粒。
但我们发现，同一个合法序列，可能前面人为添加的j不同，它都能匹配完，这样就算重了。所以我们要建立意义对应关系，我们考虑一个合法序列只能从人为添加最少的j转移过来（如果更少它就匹配不完），也就是在转移途中会至少经历一次j=0（如果没有经历过说明人为添加的j还可以再减）。因为不合法的序列前面怎么分配j都匹配不完，于是这是可行的。
设计状态f[i][j][0/1]表示处理完i组奇偶，还剩j个偶数圆粒（n-j个偶数皱粒），未经历/经历过j=0的方案数。转移见程序。
注意最后要把误转移的f[i][0][0]清空。
代码：

#include<iostream>
#include<cstdio>
#define ll long long
using namespace std;
int n,m,p;
ll f[3100][3100][2];
int main()
{
    scanf("%d%d%d",&n,&m,&p);
    for(int j=1;j<=n;j++)
        f[0][j][0]=1;
    f[0][0][1]=1;   
    for(int i=1;i<=m;i++)
    {
        for(int j=0;j<=n;j++)
        {
            f[i][j][0]=f[i][j][1]=0;
            if(j>0) f[i][j][1]+=f[i-1][j][1];
            if(j==1) f[i][j][1]+=f[i-1][j][0];
            if(j>1) f[i][j][0]+=f[i-1][j][0];
            //圆圆

            if(j<n) {f[i][j][1]+=f[i-1][j][1];f[i][j][0]+=f[i-1][j][0];}
            //皱皱

            if(j>0)f[i][j][1]+=f[i-1][j-1][1];
            if(j>1)f[i][j][0]+=f[i-1][j-1][0];
            if(j==1)f[i][j][1]+=f[i-1][j-1][0];
            //皱圆        

            if(j<n)
            {
                if(j==0) f[i][j][1]+=f[i-1][j+1][0];
                f[i][j][1]+=f[i-1][j+1][1];
                f[i][j][0]+=f[i-1][j+1][0];
            }   
            //圆皱
            f[i][j][0]%=p;f[i][j][1]%=p;
        }
        f[i][0][0]=0;
    }

    ll ans=0;
    for(int i=0;i<=n;i++)
            ans=(ans+f[m][i][1])%p;
    printf("%lld",ans);                 
    return 0;
}

整合后代码：

#include<iostream>
#include<cstdio>
#define ll long long
using namespace std;
int n,m,p;
ll f[3100][3100][2];
int main()
{
    scanf("%d%d%d",&n,&m,&p);
    for(int j=1;j<=n;j++)
        f[0][j][0]=1;
    f[0][0][1]=1;   
    for(int i=1;i<=m;i++)
    {
        for(int j=0;j<=n;j++)
        {
            f[i][j][0]=f[i][j][1]=0;
            if(j==0) f[i][j][1]+=f[i-1][j+1][0];
            if(j>0) f[i][j][1]+=f[i-1][j][1]+f[i-1][j-1][1];
            if(j==1) f[i][j][1]+=f[i-1][j][0]+f[i-1][j-1][0];
            if(j>1) f[i][j][0]+=f[i-1][j][0]+f[i-1][j-1][0];
            if(j<n) 
            {
                f[i][j][1]+=f[i-1][j][1]+f[i-1][j+1][1];
                f[i][j][0]+=f[i-1][j][0]+f[i-1][j+1][0];
            }
            f[i][j][0]%=p;f[i][j][1]%=p;
        }
        f[i][0][0]=0;
    }

    ll ans=0;
    for(int i=0;i<=n;i++)
            ans=(ans+f[m][i][1])%p;
    printf("%lld",ans);                 
    return 0;
}