hdu 5411 CRB and Puzzle 矩阵快速幂

最新推荐文章于 2019-07-09 12:53:00 发布
最新推荐文章于 2019-07-09 12:53:00 发布
阅读量555 收藏
点赞数
分类专栏： 快速幂 文章标签： 矩阵乘法 快速幂
版权声明：本文为博主原创文章，遵循 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 4462
什么是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 1657
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 695
题目链接： 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 385
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 65
链接 题解链接: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 1105
题目连接: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 4452
缘起：在学这个之前学了好一阵子矩阵，仍然不太懂。。后来发现解斐波那契数列好像只需要方形矩阵就行了，就学了一下方形矩阵的乘法。 先上代码吧 #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 351
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 596
Problem Description Farmer John有n头奶牛. 某天奶牛想要数一数有多少头奶牛,以一种特殊的方式: 第一头奶牛为1号,第二头奶牛为2号,第三头奶牛之后,假如当前奶牛是第n头,那么他的编号就是2倍的第n-2头奶牛的编号加上第n-1头奶牛的编号再加上自己当前的n的三次方为自己的编号. 现在Farmer John想知道,第n头奶牛的编号是多少,估计答案会很大,你只要输出答案对...
Java程序员面试必备知识点总结
最新发布
09-20
内容概要：本文涵盖了一系列重要的Java基础知识与进阶概念的面试问题及其详细解答，其中包括JVM的工作原理，垃圾回收的实施细节，多线程的解决方案，以及诸如单例设计模式、异常处理和线程安全集合等实际应用中的编程技巧介绍。 适合人群：旨在为准备面试或希望提高Java编程技能的专业人士提供有用指南，无论是对于初级程序员还是资深开发者都有所帮助。 使用场景及目标：适用于Java学习阶段，帮助理解和记忆相关的关键术语和复杂理论；作为快速查阅资料的宝典应对技术评估或职业机会考察。 阅读注意点：强烈建议在研读解答的同时，进行相应的编程实验和实操练习，以便于深化对每个专题的理解和巩固知识点的应用能力。
一个基于健身的社交App，内含跑步轨迹绘制，数据可视化展示，即时通讯，看新闻等模块(毕设&课设&实训&大作业&竞赛&项目)
09-20
项目工程资源经过严格测试可直接运行成功且功能正常的情况才上传，可轻松复刻，拿到资料包后可轻松复现出一样的项目，本人系统开发经验充足（全领域），有任何使用问题欢迎随时与我联系，我会及时为您解惑，提供帮助。 【资源内容】：包含完整源码+工程文件+说明（如有）等。答辩评审平均分达到96分，放心下载使用！可轻松复现，设计报告也可借鉴此项目，该资源内项目代码都经过测试运行成功，功能ok的情况下才上传的。 【提供帮助】：有任何使用问题欢迎随时与我联系，我会及时解答解惑，提供帮助 【附带帮助】：若还需要相关开发工具、学习资料等，我会提供帮助，提供资料，鼓励学习进步 【项目价值】：可用在相关项目设计中，皆可应用在项目、毕业设计、课程设计、期末/期中/大作业、工程实训、大创等学科竞赛比赛、初期项目立项、学习/练手等方面，可借鉴此优质项目实现复刻，设计报告也可借鉴此项目，也可基于此项目来扩展开发出更多功能 下载后请首先打开README文件（如有），项目工程可直接复现复刻，如果基础还行，也可在此程序基础上进行修改，以实现其它功能。供开源学习/技术交流/学习参考，勿用于商业用途。质量优质，放心下载使用。
基于java的考研指导平台设计与实现.docx
09-20
基于java的考研指导平台设计与实现.docx
asp网上办公管理系统设计(源代码+论文).zip
09-20
1、资源项目源码均已通过严格测试验证，保证能够正常运行； 2、项目问题、技术讨论，可以给博主私信或留言，博主看到后会第一时间与您进行沟通； 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用，尤其对于人工智能、计算机科学与技术等相关专业，更为适合；、下 4载使用后，可先查看README.md或论文文件（如有），本项目仅用作交流学习参考，请切勿用于商业用途。 5、资源来自互联网采集，如有侵权，私聊博主删除。 6、可私信博主看论文后选择购买源代码。 1、资源项目源码均已通过严格测试验证，保证能够正常运行； 2、项目问题、技术讨论，可以给博主私信或留言，博主看到后会第一时间与您进行沟通； 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用，尤其对于人工智能、计算机科学与技术等相关专业，更为适合；、下载 4使用后，可先查看README.md或论文文件（如有），本项目仅用作交流学习参考，请切勿用于商业用途。 5、资源来自互联网采集，如有侵权，私聊博主删除。 6、可私信博主看论文后选择购买源代码。 1、资源项目源码均已通过严格测试验证，保证能够正常运行； 2、项目问题、技术讨论，可以给博主私信或留言，博主看到后会第一时间与您进行沟通； 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用，尤其对于人工智能、计算机科学与技术等相关专业，更为适合；、下载 4使用后，可先查看README.md或论文文件（如有），本项目仅用作交流学习参考，请切勿用于商业用途。 5、资源来自互联网采集，如有侵权，私聊博主删除。 6、可私信博主看论文后选择购买源代码。
Aurora情绪树洞是面向大学生为其提供心理健康辅助和情绪宣泄的网站。(毕设&课设&实训&大作业&竞赛&项目)
09-20
项目工程资源经过严格测试可直接运行成功且功能正常的情况才上传，可轻松复刻，拿到资料包后可轻松复现出一样的项目，本人系统开发经验充足（全领域），有任何使用问题欢迎随时与我联系，我会及时为您解惑，提供帮助。 【资源内容】：包含完整源码+工程文件+说明（如有）等。答辩评审平均分达到96分，放心下载使用！可轻松复现，设计报告也可借鉴此项目，该资源内项目代码都经过测试运行成功，功能ok的情况下才上传的。 【提供帮助】：有任何使用问题欢迎随时与我联系，我会及时解答解惑，提供帮助 【附带帮助】：若还需要相关开发工具、学习资料等，我会提供帮助，提供资料，鼓励学习进步 【项目价值】：可用在相关项目设计中，皆可应用在项目、毕业设计、课程设计、期末/期中/大作业、工程实训、大创等学科竞赛比赛、初期项目立项、学习/练手等方面，可借鉴此优质项目实现复刻，设计报告也可借鉴此项目，也可基于此项目来扩展开发出更多功能 下载后请首先打开README文件（如有），项目工程可直接复现复刻，如果基础还行，也可在此程序基础上进行修改，以实现其它功能。供开源学习/技术交流/学习参考，勿用于商业用途。质量优质，放心下载使用。
矩阵的初等变换 hdu
07-27
初等变换是指对矩阵进行一系列基本操作，包括交换两行（列），将某行（列）乘以一个非零标量，将某行（列）的倍数加到另一行（列）。这些操作可以用来简化矩阵，求解线性方程组，计算矩阵的秩等。 在求解线性方程组时，我们可以通过初等变换将增广矩阵化为行阶梯形或简化行阶梯形，从而简化求解过程。以下是一些常见的初等变换： 1. 交换两行（列）：将矩阵中的两行（列）互换位置。 2. 乘以一个非零标量：将矩阵中的某一行（列）上的每个元素都乘以一个非零实数。 3. 将某行（列）的倍数加到另一行（列）：将矩阵中的某一行（列）的每个元素乘以一个实数后，加到另一行（列）对应位置上。 通过这些初等变换，我们可以改变矩阵的形式，使得它更容易进行计算或求解。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值