当前位置: 代码迷 >> 综合 >> Gym 101061G Repeat it (矩阵快速幂+Dp递推式子)
  详细解决方案

Gym 101061G Repeat it (矩阵快速幂+Dp递推式子)

热度:97   发布时间:2023-11-15 11:41:14.0

题目链接:http://codeforces.com/gym/101061/problem/G

题目大意

给定m和n,要求m个n数字形式拼接起来
其表示的数模1e9+7

题目分析 

矩阵快速幂水题,
递推式子很好列,用二维矩阵快速幂优化下即可.见代码

#include<bits/stdc++.h>
using namespace std;#define debug puts("YES");
#define rep(x,y,z) for(int (x)=(y);(x)<(z);(x)++)
#define ll long long#define lrt int l,int r,int rt
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define root l,r,rt
#define mst(a,b) memset((a),(b),sizeof(a))
#define pii pair<int,int>
#define fi first
#define se second
#define mk(x,y) make_pair(x,y)
const int mod=1e9+7;
const int maxn=1e2+5;
const int ub=1e6;
const int INF=1e9;
ll powmod(ll x,ll y){ll t; for(t=1;y;y>>=1,x=x*x%mod) if(y&1) t=t*x%mod; return t;}
ll gcd(ll x,ll y){return y?gcd(y,x%y):x;}
/*
题目大意:
给定m和n,要求m个n数字形式拼接起来
其表示的数模1e9+7题目分析:
矩阵快速幂水题,
递推式子很好列,用二维矩阵快速幂优化下即可.见代码
*/
int n,m;
struct mat{ll a[3][3];mat(){mst(a,0);}mat operator*(const mat& y)const{mat ret;rep(i,0,3) rep(j,0,3) rep(k,0,3){ret.a[i][j]=(ret.a[i][j]+a[i][k]*y.a[k][j]%mod)%mod;}return ret;}
};
mat powmat(mat x,ll y){mat ret;rep(i,0,3) ret.a[i][i]=1;for(;y;y>>=1,x=x*x) if(y&1)ret=x*ret;return ret;
}
int main(){int t;scanf("%d",&t);ll x=1;while(t--){scanf("%d%d",&m,&n);x=1;while(x<=n) x=x*10;mat t;t.a[0][0]=x%mod,t.a[0][1]=n,t.a[1][1]=1;printf("%lld\n",powmat(t,m).a[0][1]);}return 0;
}

 

  相关解决方案