hdu 5411 CRB and Puzzle 矩阵快速幂

最新推荐文章于 2019-07-09 12:53:00 发布
最新推荐文章于 2019-07-09 12:53:00 发布
阅读量556 收藏
点赞数
分类专栏： 快速幂 文章标签： 矩阵乘法 快速幂
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/inf_force/article/details/48003077
版权

题意：有n种数量无限的方块，每种方块后可以拼接另一个方块，但是只能拼接给定种类的方块，问最多拼接M次能拼接出多少种不同的样式。

这道题中n只有50，很容易想到矩阵快速幂，毕竟这种题已经烂大街了。然而比赛的时候我居然没有读到这道题！！！！！
重点在于构造转移矩阵，显然如果a后面可以拼接b，那么有 tmp[a][b]=1; 代表a->b有一种转移方式，构造完了之后，对这个矩阵自乘m次，那么矩阵的结果 tmp′[a][b] 代表的就是起点为a，终点为b，经过m次转移的方法数，但是我们要求0次转移 1次转移 2次转移 。。。M次转移的总和，所以把最后一列赋为1，就可以了。
第一次写的时候矩阵构造太挫T了，然后优化了一点勉强过了，
第二次和第三次都是优化过的。
优化过的矩阵构造。
代码1：

//author: CHC
//First Edit Time:  2015-08-25 21:13
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <set>
#include <vector>
#include <map>
#include <queue>
#include <set>
#include <algorithm>
#include <limits>
//#include <time.h>
using namespace std;
typedef long long LL;
const int INF = numeric_limits<int>::max();
const LL LL_INF= numeric_limits<LL>::max();
const int N = 53;
const int mod=2015;
int T,n,m,k,tn;
//C=A*B
void matric_mul(int A[][N],int B[][N],int C[][N]){
    int t[N][N]={0};
    for(int k=0;k<tn;k++){
        for(int i=0;i<tn;i++){
            if(A[i][k]==0)continue;
            for(int j=0;j<tn;j++){
                t[i][j]+=A[i][k]*B[k][j];
                t[i][j]%=mod;
            }
        }
    }
    for(int i=0;i<tn;i++){
        for(int j=0;j<tn;j++){
            C[i][j]=t[i][j];
        }
    }
}
int res[N][N],tmp[N][N];
void print(int A[][N]){
    for(int i=0;i<tn;i++,puts(""))
        for(int j=0;j<tn;j++)printf("%d ",A[i][j]);
    puts("");
}
int q_pow(int nn){
    //printf("res:\n");
    //print(res);
    while(nn>0){
        if(nn&1)matric_mul(res,tmp,res);
        //printf("res:\n");
        //print(res);
        nn>>=1;
        matric_mul(tmp,tmp,tmp);
    }
    int cnt=res[0][n];
    //for(int i=n;i<tn;i++)cnt+=res[0][i];
    return (cnt+1)%mod;
}
void Rd(int &x){
    char ch;
    while((ch=getchar())&&!isdigit(ch));
    x=0;
    do {
        x=(x<<3)+(x<<1)+(ch^48);
    }while((ch=getchar())&&isdigit(ch));
}
int main()
{
    //freopen("1006.in","r",stdin);
    //freopen("10061.out","w",stdout);
    //clock_t now=clock();
    int maxn=0,maxm=0;
    scanf("%d",&T);
    while(T--){
        //scanf("%d%d",&n,&m);
        Rd(n),Rd(m);
        maxn=max(n,maxn),maxm=max(m,maxm);
        tn=n+1;
        memset(tmp,0,sizeof(tmp));
        memset(res,0,sizeof(res));
        for(int i=0,x;i<n;i++){
            //scanf("%d",&k);
            Rd(k);
            for(int j=0;j<k;j++){
                //scanf("%d",&x);
                Rd(x);
                --x;
                tmp[i][x]=1;
                //tmp[i][x+n]=1;
            }
        }
        //for(int i=n;i<tn;i++)tmp[i][i]=1;
        for(int i=0;i<=n;i++)tmp[i][n]=1;
        //print(tmp);
        for(int i=0;i<n;i++)res[0][i]=1;
        res[0][n]=0;
        //print(res);
        printf("%d\n",q_pow(m));
    }
    //printf("%d %d\n",maxn,maxm);
    //printf("%.2lf\n",1.0*(clock()-now)/CLOCKS_PER_SEC);
    return 0;
}
/*
   10
   3 10
   1 2
   1 3
   0
   */

代码2:

//author: CHC
//First Edit Time:  2015-08-25 21:13
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <set>
#include <vector>
#include <map>
#include <queue>
#include <set>
#include <algorithm>
#include <limits>
//#include <time.h>
using namespace std;
typedef long long LL;
const int INF = numeric_limits<int>::max();
const LL LL_INF= numeric_limits<LL>::max();
const int N = 52;
const int mod=2015;
int T,n,m,k,tn;
//C=A*B
void matric_mul(int A[][N],int B[][N],int C[][N]){
    int t[N][N]={0};
    for(int k=0;k<tn;k++){
        for(int i=0;i<tn;i++){
            if(A[i][k]==0)continue;
            for(int j=0;j<tn;j++){
                t[i][j]+=A[i][k]*B[k][j];
                t[i][j]%=mod;
            }
        }
    }
    for(int i=0;i<tn;i++){
        for(int j=0;j<tn;j++){
            C[i][j]=t[i][j];
        }
    }
}
int res[N][N],tmp[N][N];
void print(int A[][N]){
    for(int i=0;i<tn;i++,puts(""))
        for(int j=0;j<tn;j++)printf("%d ",A[i][j]);
    puts("");
}
int q_pow(int nn){
    //printf("res:\n");
    //print(res);
    while(nn>0){
        if(nn&1)matric_mul(res,tmp,res);
        //printf("res:\n");
        //print(res);
        nn>>=1;
        matric_mul(tmp,tmp,tmp);
    }
    int cnt=0;
    for(int i=0;i<tn;i++)
        for(int j=0;j<tn;j++)cnt+=res[i][j];
    return (cnt)%mod;
}
void Rd(int &x){
    char ch;
    while((ch=getchar())&&!isdigit(ch));
    x=0;
    do {
        x=(x<<3)+(x<<1)+(ch^48);
    }while((ch=getchar())&&isdigit(ch));
}
int main()
{
    //freopen("1006.in","r",stdin);
    //freopen("10061.out","w",stdout);
    //clock_t now=clock();
    //int maxn=0,maxm=0;
    scanf("%d",&T);
    while(T--){
        //scanf("%d%d",&n,&m);
        Rd(n),Rd(m);
        //maxn=max(n,maxn),maxm=max(m,maxm);
        tn=n+1;
        memset(tmp,0,sizeof(tmp));
        memset(res,0,sizeof(res));
        for(int i=0;i<=n;i++)res[i][i]=1;
        for(int i=0,x;i<n;i++){
            //scanf("%d",&k);
            Rd(k);
            for(int j=0;j<k;j++){
                //scanf("%d",&x);
                Rd(x);
                --x;
                tmp[i][x]=1;
                //tmp[i][x+n]=1;
            }
        }
        for(int i=0;i<=n;i++)tmp[i][n]=1;
        //for(int i=n;i<tn;i++)tmp[i][i]=1;
        //print(tmp);
        //for(int i=0;i<tn;i++)res[0][i]=1;
        //print(res);
        printf("%d\n",q_pow(m-1));
    }
    //printf("%d %d\n",maxn,maxm);
    //printf("%.2lf\n",1.0*(clock()-now)/CLOCKS_PER_SEC);
    return 0;
}
/*
   10
   3 10
   1 2
   1 3
   0
   */

代码3:

//author: CHC
//First Edit Time:  2015-08-25 21:13
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <set>
#include <vector>
#include <map>
#include <queue>
#include <set>
#include <algorithm>
#include <limits>
//#include <time.h>
using namespace std;
typedef long long LL;
const int INF = numeric_limits<int>::max();
const LL LL_INF= numeric_limits<LL>::max();
const int N = 100;
const int mod=2015;
int T,n,m,k,tn;
//C=A*B
void matric_mul(int A[][N],int B[][N],int C[][N]){
    int t[N][N]={0};
    for(int k=0;k<tn;k++){
        for(int i=0;i<tn;i++){
            if(A[i][k]==0)continue;
            for(int j=0;j<tn;j++){
                t[i][j]+=A[i][k]*B[k][j];
                t[i][j]%=mod;
            }
        }
    }
    for(int i=0;i<tn;i++){
        for(int j=0;j<tn;j++){
            C[i][j]=t[i][j];
        }
    }
}
int res[N][N],tmp[N][N];
void print(int A[][N]){
    for(int i=0;i<tn;i++,puts(""))
        for(int j=0;j<tn;j++)printf("%d ",A[i][j]);
    puts("");
}
int q_pow(int nn){
    //printf("res:\n");
    //print(res);
    while(nn>0){
        if(nn&1)matric_mul(res,tmp,res);
        //printf("res:\n");
        //print(res);
        nn>>=1;
        matric_mul(tmp,tmp,tmp);
    }
    int cnt=0;
    for(int i=n;i<tn;i++)cnt+=res[0][i];
    return (cnt+1)%mod;
}
void Rd(int &x){
    char ch;
    while((ch=getchar())&&!isdigit(ch));
    x=0;
    do {
        x=(x<<3)+(x<<1)+(ch^48);
    }while((ch=getchar())&&isdigit(ch));
}
int main()
{
    //freopen("1006.in","r",stdin);
    //freopen("10061.out","w",stdout);
    //clock_t now=clock();
    int maxn=0,maxm=0;
    scanf("%d",&T);
    while(T--){
        //scanf("%d%d",&n,&m);
        Rd(n),Rd(m);
        maxn=max(n,maxn),maxm=max(m,maxm);
        tn=2*n;
        memset(tmp,0,sizeof(tmp));
        memset(res,0,sizeof(res));
        for(int i=0,x;i<n;i++){
            //scanf("%d",&k);
            Rd(k);
            for(int j=0;j<k;j++){
                //scanf("%d",&x);
                Rd(x);
                --x;
                tmp[i][x]=1;
                tmp[i][x+n]=1;
            }
        }
        for(int i=n;i<tn;i++)tmp[i][i]=1;
        //print(tmp);
        for(int i=0;i<tn;i++)res[0][i]=1;
        //print(res);
        printf("%d\n",q_pow(m-1));
    }
    //printf("%d %d\n",maxn,maxm);
    //printf("%.2lf\n",1.0*(clock()-now)/CLOCKS_PER_SEC);
    return 0;
}
/*
   10
   3 10
   1 2
   1 3
   0
   */
CHCXCHC
关注
专栏目录
中心频点计算公式_5G中的频点计算及实例分析
weixin_39628256的博客
12-21 4465
什么是ARFCNARFCN，英文全称Absolute Radio Frequency Channel Number，即绝对无线频道编号，是指在GSM无线系统中用来鉴别特殊射频通道的编号方案。“ARFCN”一词源自于GSM技术，随着新技术的发展，延伸出其他类似术语，如UMTS / WCDMA的UARFCN，E-UTRAN / LTE的EARFCN，以及现在5G / NR的NR-ARFCN。计算公式：...
HDU 5411 CRB and Puzzle（矩阵快速幂+可达矩阵）
a novice
08-22 1659
HDU 5411题意：Count the number of different patterns by counting the number of different paths of length at most m-1.思路：其实这类问题就是求： S=I+A+A2+...+AmS=I+A+A^2+...+A^m 一般普适的计算方式是： (AmS)=(AI0I)m∗(I0)\left(
hdu5411CRB and Puzzle 矩阵快速幂
ijbuhv的专栏
10-20 330
//n个零件，一个模型最多可由m个零件组成 //在每一个零件后面只能由特定的零件组装 //问能够组成的模型的种类数 //先构造一个n*n的矩阵A，那么对于每一个 //A^(m-1)表示由m个零件组成，对应的零件的 //的种类数 //所以就是求A^1+A^2+...A^(m-1) //快速幂加一个二分可做 #include #include #include using namespace std
hdu 5411 CRB and Puzzle（矩阵快速幂）
piaocoder
08-25 696
题目链接： http://acm.hdu.edu.cn/showproblem.php?pid=5411 解题思路： 题目大意： 给定n个点 常数m 下面n行第i行第一个数字表示i点的出边数，后面给出这些出边。 问：图里存在多少条路径使得路径长度 官方题解： We can count the number of different patterns by count
HDU - 5411 CRB and Puzzle 矩阵快速幂
07-09 76
HDU - 5411 考虑直接dp会T， 用矩阵优化一下就好了。 #include<bits/stdc++.h> #define LL long long #define LD long double #define ull unsigned long long #define fi first #define se second #define mk m...
HDU 5411 CRB and Puzzle （矩阵快速幂水题）
满月照亮的路
11-09 225
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5411 #include&lt;bits/stdc++.h&gt; using namespace std; #define debug puts("YES"); #define rep(x,y,z) for(int (x)=(y);(x)&lt;(z);(x)++) #define ll lon...
hdu 5411 CRB and Puzzle【矩阵快速幂】
mengxiang000000
06-18 386
CRB and Puzzle Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 1013    Accepted Submission(s): 402 Problem Description CRB is now playi
hdu 5411 CRB and Puzzle 矩阵高速幂
weixin_30344795的博客
06-03 66
链接 题解链接:http://www.cygmasot.com/index.php/2015/08/20/hdu_5411/ 给定n个点 常数m 以下n行第i行第一个数字表示i点的出边数。后面给出这些出边。 问：图里存在多少条路径使得路径长度<=m。路径上的点能够反复。 思路: 首先能得到一个m*n*n的dp。dp[i][j]表示路径长度为i 路径的结尾为j的路径个数 。...
HDU5411CRB and Puzzle(矩阵快速幂)
追梦赤子心
08-21 904
HDU5411CRB and Puzzle(矩阵快速幂)2015 Multi-University Training Contest 10
HDU 5015 233 Matrix 【矩阵快速幂】
Tabris的博客
07-29 1112
题目连接:http://acm.hdu.edu.cn/showproblem.php?pid=5015——————————————-.233 MatrixTime Limit: 10000/5000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) Total Submission(s): 1817 Accepted
C++ 矩阵快速幂
Object_的博客
09-16 4459
缘起：在学这个之前学了好一阵子矩阵，仍然不太懂。。后来发现解斐波那契数列好像只需要方形矩阵就行了，就学了一下方形矩阵的乘法。 先上代码吧 #include&amp;amp;amp;amp;lt;cstdio&amp;amp;amp;amp;gt; #include&amp;amp;amp;amp;lt;iostream&amp;amp;amp;amp;gt; #define ll long long #define mod 1000000007 using nam
HDU1757 - A Simple Math Problem - 矩阵快速幂
终有绿洲摇曳在沙漠
03-05 355
1.题目描述 A Simple Math Problem Time Limit: 3000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 4507    Accepted Submission(s): 2716 Problem Descripti
HDU-6460 矩阵快速幂 含幂项 二项展开
刷题记录
03-17 599
Problem Description Farmer John有n头奶牛. 某天奶牛想要数一数有多少头奶牛,以一种特殊的方式: 第一头奶牛为1号,第二头奶牛为2号,第三头奶牛之后,假如当前奶牛是第n头,那么他的编号就是2倍的第n-2头奶牛的编号加上第n-1头奶牛的编号再加上自己当前的n的三次方为自己的编号. 现在Farmer John想知道,第n头奶牛的编号是多少,估计答案会很大,你只要输出答案对...
基于libev和mongoose的集成多种协议的异步嵌入式网络库Evmongoose设计源码
最新发布
09-25
该项目是采用libev和mongoose构建的，支持多协议集成的异步嵌入式网络库Evmongoose的设计源码，包含47个文件，其中包括15个Lua脚本、13个C语言源代码、4个文本文件、4个Markdown文件、2个密钥文件、2个PEM文件、2个头文件、1个Git忽略文件、1个授权文件、1个SHTML文件。Evmongoose库支持TCP、HTTP、WebSocket和MQTT等多种协议，并提供了Lua API接口，适用于嵌入式网络开发领域。
计算机软件设计-Axure RP8-账号登录页面
09-25
实现动态面板【账号密码登录|短信登录】切换 短信验证码【60-0秒】登录倒计时
bitComet，比特彗星多线程下载器
09-25
bitComet，比特彗星多线程下载器
基于GitCGB2104项目的Java与HTML综合设计源码
09-25
本项目为GitCGB2104项目，集成Java与HTML技术，总计包含85个文件，涵盖36个Java源文件、16个Java类文件、9个XML配置文件、8个HTML文件、7个YAML文件、6个Idea项目配置文件、1个SQL脚本文件、1个PNG图片文件及1个文本文件。该项目设计精巧，旨在提供全面的技术支持。
基于Java语言的校园疫情防控管理后端设计源码
09-25
该项目是一款以Java语言开发的校园疫情防控管理后端系统源码，总计包含126个文件，其中包括84个Java源代码文件、21个PNG图片文件、11个XML配置文件、2个JAR包文件、1个Git忽略文件、1个属性文件和1个Markdown文件。系统设计旨在为校园疫情防控提供高效的后台管理支持。
矩阵的初等变换 hdu
07-27
初等变换是指对矩阵进行一系列基本操作，包括交换两行（列），将某行（列）乘以一个非零标量，将某行（列）的倍数加到另一行（列）。这些操作可以用来简化矩阵，求解线性方程组，计算矩阵的秩等。 在求解线性方程组时，我们可以通过初等变换将增广矩阵化为行阶梯形或简化行阶梯形，从而简化求解过程。以下是一些常见的初等变换： 1. 交换两行（列）：将矩阵中的两行（列）互换位置。 2. 乘以一个非零标量：将矩阵中的某一行（列）上的每个元素都乘以一个非零实数。 3. 将某行（列）的倍数加到另一行（列）：将矩阵中的某一行（列）的每个元素乘以一个实数后，加到另一行（列）对应位置上。 通过这些初等变换，我们可以改变矩阵的形式，使得它更容易进行计算或求解。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值