HDU6395 多校赛Sequence

最新推荐文章于 2018-08-20 19:54:25 发布

D-DeepWave

最新推荐文章于 2018-08-20 19:54:25 发布

阅读量302

点赞数

CC 4.0 BY-SA版权

分类专栏：数论文章标签：矩阵快速幂整除分块多校赛

本文链接：https://blog.youkuaiyun.com/DEEPWAVE98/article/details/81661448

数论专栏收录该内容

11 篇文章

订阅专栏

本文介绍了解决HDU 6395问题的一种方法，通过结合整除分块思想和矩阵快速幂技术进行优化递推。提供了一段完整的C++代码实现，展示了如何有效地计算大规模数据的递推问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门----http://acm.hdu.edu.cn/showproblem.php?pid=6395

思路：由于有后半部分的下取整，可以联想到整除分块去做。再用矩阵快速幂优化递推即可。

Code：

//2
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#define LL long long
using namespace std;
const LL mod = 1000000007;
LL a, b, c, d, p, n;
struct martix
{
    LL m[4][4];
    martix()
    {
        memset(m, 0, sizeof(m));
    }
    martix operator*(const martix b)const
    {
        LL i, j, k;
        martix ans;
        for(i = 1; i <= 3; i++)
            for(j = 1; j <= 3; j++)
            {
                ans.m[i][j] = 0;
                for(k = 1; k <= 3; k++)
                {
                    ans.m[i][j] += m[i][k] * b.m[k][j];
                    ans.m[i][j] %= mod;
                }
            }
        return ans;
    }
};
void Debug (martix a)
{
    for(int i = 1; i <= 3; i++)
    {
        for(int j = 1; j <= 3; j++)
            cout << a.m[i][j] << ' ';
        cout << endl;
    }
}
martix martix_pow(LL n,LL x)
{
    martix temp, res;
    temp.m[1][1] = d;
    temp.m[1][2] = c;
    temp.m[1][3] = x;
    temp.m[2][1] = 1;
    temp.m[3][3] = 1;
    res.m[1][1] = res.m[2][2] = res.m[3][3] = 1;
    while(n)
    {
        if(n & 1)
            res = res * temp;
        temp = temp * temp;
        n >>= 1;
    }
    return res;
}
int main()
{
    int T;
    scanf("%d", &T);
    while(T--)
    {
        scanf("%lld%lld%lld%lld%lld%lld", &a, &b, &c, &d, &p, &n);
        if(n == 1)
            printf("%lld\n", a);
        else if(n == 2)
            printf("%lld\n", b);
        else
        {
            martix A, B;
            LL l, r;
            A.m[1][1] = A.m[2][2] = A.m[3][3] = 1;
            for(l = 3; l <= n; ++l)
            {
                r = (l <= p) ? min(p / (p / l), n) : n;
                A = (martix_pow(r - l + 1,p/l))*A;
                l=r;
               // Debug(A);
            }
            printf("%lld\n",(A.m[1][1]*b%mod+A.m[1][2]*a%mod+A.m[1][3])%mod);
        }
    }
}