[DP]CH Round#30 T2摆花

头一次参加比赛..还是中途发现有这个的...就加了进来....真是一个悲伤的故事...

T1一开始敲了一个，交上去才发现，是错的..没考虑从句套从句...果断cout<<"No";

T2敲了个dfs，结果还没写对，叫hzwer帮忙敲了一个矩阵乘法，结果他还把01看反了.....

T3纯模拟暴力，爆0，T4T5放弃治疗，真是一个悲伤的故事

原题：

背景及描述

艺术馆门前将摆出许多花，一共有n个位置排成一排，每个位置可以摆花也可以不摆花。有些花如果摆在相邻的位置（隔着一个空的位置不算相邻），就不好看了。假定每种花数量无限，求摆花的方案数。

输入格式

输入有1+m行，第一行有两个用空格隔开的正整数n、m，m表示花的种类数。接下来的m行，每行有m个字符1或0,若第i行第j列为1，则表示第i种花和第j种花不能排在相邻的位置，输入保证对称。（提示：同一种花可能不能排在相邻位置）。

输出格式

输出只有一个整数，为方案数（这个数字可能很大，请输出方案数除以1000000007的余数）。

样例输入

2 2
01
10

样例输出

7

样例说明

七种方案为(空，空)、（空，1）、（1、空）、（2、空）、（空、2）、（1,1）、（2,2）。

数据范围与约定

• 20%的数据，1＜n≤5,0＜m≤10。
• 60%的数据，1＜n≤200,0＜m≤100。
• 100%的数据，1＜n≤1000000000，0＜m≤100。

本蒟蒻当时dfs暴力爆0，dp暴力都不会...后来也只是看懂了dp暴力，特此献上

#include<iostream>
#include<cstdio>
#define ll long long 
#define mod 1000000007
using namespace std;
int n,m;char ch[101];
ll a[101][101],f[10001][101],ans;
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)
	{
		scanf("%s",ch);
		for(int j=1;j<=m;j++)
		   if(ch[j-1]=='0')
		      a[i][j]=1;
	}
	f[0][0]=1;
	for(int i=1;i<=n;i++)
		for(int j=0;j<=m;j++)
		   for(int k=0;k<=m;k++)
		      if(a[j][k]||j==0||k==0)
                         f[i][j]=(f[i][j]+f[i-1][k])%mod;
    for(int i=0;i<=m;i++)ans=(ans+f[n][i])%mod;
    printf("%d",ans);
	return 0;
}

再附上hzwer的矩阵乘法Orz

#include<iostream>
#include<cstdio>
#define ll long long 
#define mod 1000000007
using namespace std;
int n,m;char ch[101];
ll c[101],a[101][101],b[101][101],ans;
void mul(ll a[101][101],ll b[101][101],ll s[101][101])
{
    ll t[101][101];
    for(int i=0;i<=m;i++)
       for(int j=0;j<=m;j++)
       {
          t[i][j]=0;
          for(int k=0;k<=m;k++)
             t[i][j]=(t[i][j]+a[i][k]*b[k][j])%mod;
         }
    for(int i=0;i<=m;i++)
       for(int j=0;j<=m;j++)
          s[i][j]=t[i][j];
}
void mul2(ll a[101][101],ll b[101],ll s[101])
{
    ll t[101];
    for(int i=0;i<=m;i++)
       {
          for(int k=0;k<=m;k++)
             t[i]=(t[i]+a[i][k]*b[k])%mod;
         }
    for(int i=0;i<=m;i++)s[i]=t[i];
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)
	{
		scanf("%s",ch);
		for(int j=1;j<=m;j++)
		   if(ch[j-1]=='0')
		      a[i][j]=1;
	}
	for(int i=0;i<=m;i++)a[i][0]=a[0][i]=1;
	for(int i=0;i<=m;i++)b[i][i]=1;
	while(n)
       {
              if(n&1)mul(a,b,b);
              n>>=1;
              mul(a,a,a);
       }
       c[0]=1;
       mul2(b,c,c);
       for(int i=0;i<=m;i++){ans+=c[i];ans%=mod;}
       printf("%d\n",ans);
	return 0;
}