CodeForces 621E Wet Shark and Blocks (矩阵快速幂+水题)

题目链接:http://codeforces.com/problemset/problem/621/E

题目大意:

给定b个带有n个数的砖块,
每个砖块可以选则一个数字作为要拼凑的数字的下一位,
问有多少种选择方案可以使得凑出的数modx等于k.

题目分析: 

比较裸的矩阵快速幂。
我们不难发现矩阵的维度就是x 的限度,
100*100的矩阵条件是允许的。
然后对b进行幂次级别运算(完全符合题意好吧T_T。。。)
最重要的就是状态转移矩阵了,这个
也特别明显了,dp.a[i][j]就是表示j通过一个砖块的各种val
最后能变成i 的方案数。
比较水比较舒服了。

#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=100;
const int ub=1e6;
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;}
/*
题目大意:
给定b个带有n个数的砖块,
每个砖块可以选则一个数字作为要拼凑的数字的下一位,
问有多少种选择方案可以使得凑出的数modx等于k.

题目分析:
比较裸的矩阵快速幂。
我们不难发现矩阵的维度就是x 的限度,
100*100的矩阵条件是允许的。
然后对b进行幂次级别运算(完全符合题意好吧T_T。。。)
最重要的就是状态转移矩阵了,这个
也特别明显了,dp.a[i][j]就是表示j通过一个砖块的各种val
最后能变成i 的方案数。
比较水比较舒服了。
时间复杂度:(x*x*logb)。
*/

int n,b,k,x;
struct mat{
    ll a[maxn][maxn];
    mat(){
        mst(a,0);
    }
    mat operator*(const mat& y)const{
        mat ret;
        rep(i,0,x) rep(j,0,x) rep(p,0,x)
            ret.a[i][j]=(ret.a[i][j]+a[i][p]*y.a[p][j]%mod)%mod;
        return ret;
    }
};
mat matpow(mat y,int n){
    mat ret;rep(i,0,x) ret.a[i][i]=1;
    for(;n;y=y*y,n>>=1) if(n&1) ret=y*ret;
    return ret;
}
int val;
int main(){
    cin>>n>>b>>k>>x;
    k%=x;
    mat dp;
    rep(i,0,n){
        cin>>val;
        rep(j,0,x) dp.a[(j*10+val)%x][j]++;
    }
    dp=matpow(dp,b);
    cout<<dp.a[k][0]<<endl;
    return 0;
}

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值