最猥琐的一道题。
分两块做:
①前50%
很显然是状压dp,我们设f[i][j]表示前i行,第i行状态为j的方案数。
状态中,1表示下凸,0表示其他。
然后我们预处理出所有搭配起来合法的状态,就可以直接转移。
什么叫合法呢?就是这两个状态i和j,i&j=0因为不可能两个都下凸,这样就重叠了。
而且i|j以后,连续的0必须是偶数个,不然就不会合法,这很显然。
#include<cstdio>
#include<iostream>
#include<cstring>
#define mo 1000000007
using namespace std;
long long n,m;
long long f[105][2050];
int b[2050][1005];
bool bz[2050];
char s[110];
int main()
{
scanf("%s%lld",s+1,&m);
int i,j,x,k;
memset(bz,true,sizeof(bz));
for (i=0;i<=(1<<m)-1;++i)
{
x=i;
j=0;
k=1;
while (x||k<=m+1)
{
if ((x&1)==1||k>m)
{
if(j%2)
{
bz[i]=false;
break;
}
j=0;
}
if ((x&1)==0)
++j;
++k;
x>>=1;
}
}
for (i=0;i<=(1<<m)-1;++i)
for (j=0;j<=(1<<m)-1;++j)
if((i&j)==0&&bz[i|j]) b[i][++b[i][0]]=j;
if(strlen(s+1)<=3)
{
for (i=1;i<=strlen(s+1);++i)
n=n*10+(s[i]-48);
f[0][0]=1;
for (i=1;i<=n;++i)
for (j=0;j<=(1<<m)-1;++j)
for (k=1;k<=b[j][0];++k)
(f[i][j]+=f[i-1][b[j][k]])%=mo;
printf("%lld\n",f[n][0]%mo);
}
}
②后50%
我们发现每次我们预处理的东西p[i][j]=1表示从i这个状态转移到j这个状态合法。
那么每次转移一行都是这个所以直接就可以用一个矩阵乘法,也就是你预处理那东西,直接快速幂上n次,答案就是p[0][0]了。
#include<cstdio>
#include<iostream>
#include<cstring>
#define mo 1000000007
using namespace std;
long long n,m;
long long f[105][2050];
int b[2050][1005];
bool bz[2050];
char s[210];
long long matr[35][35],ans[35][35],matrix[35][35],ans2[35][35];
int n1[205];
int main()
{
scanf("%s%lld",s+1,&m);
int i,j,x,k;
memset(bz,true,sizeof(bz));
for (i=0;i<=(1<<m)-1;++i)
{
x=i;
j=0;
k=1;
while (x||k<=m+1)
{
if ((x&1)==1||k>m)
{
if(j%2)
{
bz[i]=false;
break;
}
j=0;
}
if ((x&1)==0)
++j;
++k;
x>>=1;
}
}
int l=0;
if(m<=5) l=1;
for (i=0;i<=(1<<m)-1;++i)
for (j=0;j<=(1<<m)-1;++j)
if((i&j)==0&&bz[i|j])
{
b[i][++b[i][0]]=j;
if(l) matr[i][j]=1;
}
if(strlen(s+1)<=3)
{
for (i=1;i<=strlen(s+1);++i)
n=n*10+(s[i]-48);
f[0][0]=1;
for (i=1;i<=n;++i)
for (j=0;j<=(1<<m)-1;++j)
for (k=1;k<=b[j][0];++k)
(f[i][j]+=f[i-1][b[j][k]])%=mo;
printf("%lld\n",f[n][0]%mo);
}
else
{
for (i=strlen(s+1);i>=1;--i)
n1[i]=s[strlen(s+1)-i+1]-48;
n1[0]=strlen(s+1);
ans[0][0]=1;
while(n1[0]>1||(n1[0]==1&&n1[1]>0))
{
if(n1[1]&1)
{
memset(ans2,0,sizeof(ans2));
for (k=0;k<=(1<<m)-1;++k)
for (i=0;i<=(1<<m)-1;++i)
for (j=0;j<=(1<<m)-1;++j)
ans2[i][j]=(ans2[i][j]+ans[i][k]*matr[k][j])%mo;
memcpy(ans,ans2,sizeof(ans2));
}
memset(matrix,0,sizeof(matrix));
for (k=0;k<=(1<<m)-1;++k)
for (i=0;i<=(1<<m)-1;++i)
for (j=0;j<=(1<<m)-1;++j)
(matrix[i][j]+=matr[i][k]*matr[k][j])%=mo;
memcpy(matr,matrix,sizeof(matrix));
if (n1[1]%2==1) n1[1]--;
for (i=1;i<n1[0];i++)
{
if (n1[i+1]%2==1) n1[i+1]--,n1[i]+=10;
n1[i]/=2;
}
n1[n1[0]]/=2;
if (n1[n1[0]]==0) n1[0]--;
}
printf("%lld\n",ans[0][0]);
}
}