地斗主

地斗主

https://ac.nowcoder.com/acm/problem/19833

题意

的骨牌去填充 的矩阵,问有几种方法。

分析

手算一下前五项,可以得到这么一个序列:1,5,11,36,95

表示拼完前 i 列的方案数。这时考虑 OEIS 。可以得到递推式:

然后利用矩阵快速幂就可以解决了。

#include <bits/stdc++.h>
using namespace std;
#define mem(a,b) memset(a,b,sizeof(a))
#define int long long

struct mat
{
    int a[4][4];
    mat(){mem(a,0);}
};

mat mul(mat a,mat b,int M)
{   
    mat res;
    for(int i = 0; i < 4; i++) 
    {
        for(int j = 0; j < 4; j++) 
        {
            for(int k = 0; k < 4; k++)
            {
                res.a[i][j] += a.a[i][k]*b.a[k][j]%M;
                res.a[i][j] += M;
                res.a[i][j] %= M;
            }
        }
    }
    return res;
}

mat solve(int k,int M)
{
    mat res;
    res.a[0][0] = 1;res.a[0][1] = 5;res.a[0][2] = 1;res.a[0][3] = -1;
    res.a[1][0] = 1;res.a[2][1] = 1;res.a[3][2] = 1;
    mat ans;
    ans.a[0][0] = 36;ans.a[1][0] = 11;ans.a[2][0] = 5;ans.a[3][0] = 1;
    while(k)
    {
        if(k&1) ans = mul(res,ans,M);
        res = mul(res,res,M);
        k >>= 1;
    }
    return ans;
}

signed main()
{
    ios::sync_with_stdio(false);cin.tie(0);
    int T;cin>>T;
    int a[] = {1,1,5,11};
    while(T--)
    {   
        int n,m;
        cin>>n>>m;
        if(n >= 4)
        {
            mat res = solve(n-4,m);
            cout<<res.a[0][0]<<endl;
        }
        else
        {
            cout<<a[n]%m<<endl;
        }
    }
    return 0;
}
全部评论

相关推荐

emo的打工鸭又被画饼了:我看你是外星人,听哥劝,你不应该来地球找工作的,地球要996的
点赞 评论 收藏
分享
2024-12-29 11:08
湖南工业大学 Java
程序员牛肉:简历没什么大问题了。 而且不要再换项目了。三月份就开暑期实习了,现在都一月份了。实在来不及重新开一下项目了。把一个项目写完或许很快,但是把一个项目搞懂吃透并不简单。所以不要换项目了,把你简历上面的两个项目好好挖一挖吧。 具体 体现在:你能不能流利的说出你的项目的每一个功能点代码实现?你能不能说出在这块除了A技术之外,还有其他技术能够实现嘛?如果有其他技术能够实现,那你这块为什么选择了你当前用的这个技术?
投递牛客等公司
点赞 评论 收藏
分享
评论
点赞
收藏
分享

创作者周榜

更多
牛客网
牛客企业服务