CH Round #30 摆花[矩阵乘法]

Posted Candy?

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了CH Round #30 摆花[矩阵乘法]相关的知识,希望对你有一定的参考价值。

摆花 CH Round #30 - 清明欢乐赛

背景及描述

艺术馆门前将摆出许多花,一共有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。

 

正解:

这样来看这个问题,我们先定义一种新的花:不摆花,然后我们把所有种类的

花看成一个点,在不互相冲突的种类之间连一条边,其中不摆花不和所有花冲突。我们要求 摆到n个位置,实际上可以认为是求这个图中长度为n的路径的条数,这样把问题转换成了 经典问题,可以运用矩阵乘法求解,得 100 分 

 

关于图中长度为n路径计数:

见图片

技术分享

 

注意:本题必须边从0开始

//
//  main.cpp
//  ch30B
//
//  Created by Candy on 23/10/2016.
//  Copyright © 2016 Candy. All rights reserved.
//

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
const int M=105,MOD=1000000007;
typedef long long ll;
inline int read(){
    char c=getchar();int x=0,f=1;
    while(c<0||c>9){if(c==-)f=-1;c=getchar();}
    while(c>=0&&c<=9){x=x*10+c-0;c=getchar();}
    return x*f;
}
int n,m;
char s[M];
struct mat{
    ll mt[M][M];
    mat(){memset(mt,0,sizeof(mt));}
}g,im;
mat mult(mat &x,mat &y){
    mat t;
    for(int i=0;i<=m;i++)
        for(int k=0;k<=m;k++) if(x.mt[i][k])
            for(int j=0;j<=m;j++)
                t.mt[i][j]=(t.mt[i][j]+x.mt[i][k]*y.mt[k][j]%MOD)%MOD;
    return t;
}
void dp(){
    for(int i=0;i<=m;i++)
            im.mt[i][i]=1;
    
    for(;n;n>>=1,g=mult(g,g))
        if(n&1) im=mult(im,g);
}
int main(int argc, const char * argv[]){
    //freopen("harem.in","r",stdin);
    //freopen("harem.out","w",stdout);
    
    n=read();m=read();
    for(int i=1;i<=m;i++){
        scanf("%s",s+1);
        for(int j=1;j<=m;j++) g.mt[i][j]=(s[j]-0)^1;
    }
    for(int i=0;i<=m;i++) g.mt[i][0]=g.mt[0][i]=1;
    dp();
    ll ans=0;
    for(int i=0;i<=m;i++) ans=(ans+im.mt[0][i])%MOD;
    printf("%lld",ans);
}

 

以上是关于CH Round #30 摆花[矩阵乘法]的主要内容,如果未能解决你的问题,请参考以下文章

算法导论—矩阵链乘法(动态规划)

1571:基础练习 矩阵乘法

Clojure 与 Numpy 中的矩阵乘法

如何在 python 中并行化以下代码片段?

矩阵乘法2 (Codevs No.3147)

试题 基础练习 矩阵乘法暴力矩阵乘法,再次了解一下矩阵