【题解】HDU Homework(倍增)

矩阵题一定要多多检查一下是否行列反了...

一百个递推项一定要存101个

说多了都是泪啊

一下午就做了这一道题因为实在是太菜了太久没写这种矩阵的题目...

设一个行向量\(e\),和一个增逛矩阵\(A\),他们咋定义的见我那篇讲线性递推博客

现在我们再预处理\(st\)矩阵数组,其中\(st_i=A^{2^i}\)。

考虑这样一种做法,我们考虑让\(e\)总共有101个值,然后当第一个值被增逛为\(f_{q.n-100}\)时,暴力将\(e\)中的第\(101\)项修改。中间的转移直接利用这个倍增数组。中间的转移由于是倍增,而且是行乘以一个矩阵,这样的话复杂度就是\(O(m^2)\)。考虑一个询问就要处理一次最终复杂度\(O(m^3\log n+qm^2\log n)\)

细节要处理一下

//@winlere
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector> using namespace std; typedef long long ll;
inline int qr(){
register int ret=0,f=0;
register char c=getchar();
while(c<48||c>57)f|=c==45,c=getchar();
while(c>=48&&c<=57) ret=ret*10+c-48,c=getchar();
return f?-ret:ret;
}
const int maxn=1e2+3;
const int mod=1e9+7;
int n,m,q,T;
const int M=101;
struct MAT{
int data[maxn][maxn];
MAT(){memset(data,0,sizeof data);}
inline int*operator [](int x){return data[x];}
inline MAT operator *(MAT&f){
MAT ret;
for(int k=1;k<=M;++k)
for(int t=1;t<=M;++t)
for(int i=1;i<=M;++i)
ret[t][i]=(ret[t][i]+1ll*data[t][k]*f[k][i])%mod;
return ret;
}
}st[31]; struct line{
int data[maxn];
line(){memset(data,0,sizeof data);}
inline int&operator[](int x){return data[x];}
inline line operator *(MAT&f){
line ret;
for(int i=1;i<=M;++i)
for(int t=1;t<=M;++t)
ret[t]=(ret[t]+1ll*data[i]*f[i][t])%mod;
return ret;
}
}e; int init[maxn];
pair<pair<int,int>,vector<int> > Q[maxn]; inline void pow(const int&num){
for(int t=0;t<30;++t)
if(num>>t&1)
e=e*st[t];
} int main(){
int cnt=0;
while(~scanf("%d%d%d",&n,&m,&q)){
memset(st[0].data,0,sizeof st[0].data);
memset(e.data,0,sizeof e.data);
for(int t=1;t<=m;++t) e[t]=init[t]=qr()%mod;
for(int t=1;t<M;++t) st[0][t+1][t]=1;
T=qr();
for(int t=1;t<=T;++t) st[0][M-t+1][M]=qr()%mod;
for(int t=1;t<30;++t) st[t]=st[t-1]*st[t-1];
for(int t=1;t<=q;++t){
Q[t].first.first=qr();
Q[t].first.second=qr();
Q[t].second.clear();
Q[t].second.push_back(0);
for(int i=1;i<=Q[t].first.second;++i) Q[t].second.push_back(qr());
}
sort(Q+1,Q+q+1);
int l=1;
for(int t=m+1;t<=M;++t){
if(t==Q[l].first.first&&l<=q){
int ret=0;
for(int i=1;i<=Q[l].first.second;++i)
ret=(ret+1ll*init[t-i]*Q[l].second[i])%mod;
init[t]=ret;
++l;
}
else {
int ret=0;
for(int i=1;i<=T;++i)
ret=(ret+1ll*init[t-i]*st[0][M-i+1][M])%mod;
init[t]=ret;
}
}
for(int t=1;t<=M;++t) e[t]=init[t];
if(n<=M) {printf("Case %d: %d\n",++cnt,init[n]); continue;}
int cur=1;
for(;l<=q;++l){
if(Q[l].first.first>n) break;
pow(Q[l].first.first-(cur+M-1));
cur+=Q[l].first.first-(cur+M-1);
e[M]=0;
for(int t=1;t<=Q[l].first.second;++t)
e[M]=(e[M]+1ll*Q[l].second[t]*e[M-t])%mod;
}
pow(n-cur);
printf("Case %d: %d\n",++cnt,e[1]);
}
return 0;
}

【题解】HDU Homework(倍增)的更多相关文章

  1. [题解]hdu 1009 FatMouse' Trade(贪心基础题)

    Problem Description FatMouse prepared M pounds of cat food, ready to trade with the cats guarding th ...

  2. HDU 2586 倍增法求lca

    How far away ? Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)To ...

  3. hdu 6107--Typesetting(倍增)

    题目链接 Problem Description Yellowstar is writing an article that contains N words and 1 picture, and t ...

  4. 题解 HDU 3698 Let the light guide us Dp + 线段树优化

    http://acm.hdu.edu.cn/showproblem.php?pid=3698 Let the light guide us Time Limit: 5000/2000 MS (Java ...

  5. 题解 HDU 5279 YJC plays Minecraft

    题目传送门 题目大意 给出\(n\)以及\(a_{1,2,...,n}\),表示有\(n\)个完全图,第\(i\)个完全图大小为\(a_i\),这些完全图之间第\(i\)个完全图的点\(a_i\)与\ ...

  6. 题解-bzoj4320 Homework

    Problem bzoj4320 Solution 前置技能:分块+线段树+卡常+一点小小的数学知识 考试时A的 这种题无论怎么处理总有瓶颈,套路分块,设\(k\)以下的插入时直接暴力预处理,查询时直 ...

  7. 题解——HDU 4734 F(x) (数位DP)

    这道题还是关于数位DP的板子题 数位DP有一个显著的特征,就是求的东西大概率与输入关系不大,理论上一般都是数的构成规律 然后这题就是算一个\( F(A) \)的公式值,然后求\( \left [ 0 ...

  8. 题解——HDU 2089 不要62(数位DP)

    最近在学数位DP 应该是入门题吧 设\( dp[i][0/1] \)表示到第\( i \)位时,前一位是否是6的满足条件的数的个数 然后就是套路 注意\( limit \)的限制条件以及转移时候信息的 ...

  9. 题解——HDU 1848 Fibonacci again and again

    一道组合游戏的题目 SG函数的板子题 预处理出SG函数的值然后回答询问即可 代码 #include <cstdio> #include <algorithm> #include ...

随机推荐

  1. @NOIP2018 - D1T2@ 货币系统

    目录 @题目描述@ @题解@ @代码@ @题目描述@ 在网友的国度中共有 n 种不同面额的货币,第 i 种货币的面额为 a[i],你可以假设每一种货币都有无穷多张.为了方便,我们把货币种数为 n.面额 ...

  2. python基础---集合类型(Sets)

    集合类型(Sets) 集合对象是不同的(不可重复)hashable对象的无序集合.常见用法包括:成员关系测试.移除序列中的重复.以及科学计算,例如交集.并集.差分和对称差分.通俗点来说,集合是一个无序 ...

  3. oracle函数 nls_charset_id(c1)

    [功能]返回字符集名称参应id值 [参数]c1,字符型 [返回]数值型 sql> select nls_charset_id('zhs16gbk') from dual; nls_charset ...

  4. @bzoj - 2388@ 旅行规划

    目录 @description@ @solution@ @accepted code@ @details@ @description@ 请你维护一个序列,支持两种操作: (1)某个区间 [x, y] ...

  5. HZOJ 题

    首先对于n<=100的点,直接暴力dp,f[i][j][k]表示时间为i,在i,j位置的方案数,枚举转移即可,期望得分40. ) { ) { f[][][]=; ;i<=n;i++) ;x ...

  6. Android ListView显示底部的分割线

    有些时候,我们会提出这样的需求,希望ListView显示底部(顶部)的分割线,这样做,会使得UI效果更加精致,如下图所示: 如果搜索资料,大家会搜到一堆相关的方法,最多的莫过于设置listview的f ...

  7. 如何用JS和HTML 做一个桌面炒股小插件【原创】

    首先,使用node-webkit 做环境,废话不多说,直接贴HTML <!DOCTYPE html> <html xmlns="http://www.w3.org/1999 ...

  8. hdu 1277 全文检索 (直接映射查找 || 自动机)

    Problem - 1277 无聊做水题的时候发现的一道题目.这道题第一反应可以用自动机来解决.当然,条件是各种限制,从而导致可以用直接映射标记的方法来搜索.具体的做法就像RK算法一样,将字符串has ...

  9. DispatcherTimer 应用实例

    public partial class MainWindow : Window { public MainWindow() { InitializeComponent(); //实例化 Dispat ...

  10. oracle用WHERE替代ORDER BY

    ORDER BY 子句只在两种严格的条件下使用索引. ORDER BY中所有的列必须包含在相同的索引中并保持在索引中的排列顺序. ORDER BY中所有的列必须定义为非空. WHERE子句使用的索引和 ...