栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 面试经验 > 面试问答

zoj 3772 Calculate the Function

面试问答 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

zoj 3772 Calculate the Function

#include <iostream>#include <cstring>#include <stdio.h>using namespace std;const int MAX_N = 100000 + 100;const int MOD = 1000000007;struct Matrix{    long long m[2][2];    Matrix(){memset(m, 0, sizeof(m));}    Matrix operator * (const Matrix &b)    {        Matrix temp;        for(int i = 0; i < 2; i++) for(int j = 0; j < 2; j++)     for(int k = 0; k < 2; k++)         temp.m[i][j] = ((m[i][k] * b.m[k][j]) + temp.m[i][j]) % MOD;        return temp;    }};struct Node{    int l, r;    Matrix mat;};Node node[MAX_N << 2];int a[MAX_N];void Build(int rt, int l, int r){    int m;    node[rt].l = l;    node[rt].r = r;    if(l == r)    {        node[rt].mat.m[0][0] = 1;        node[rt].mat.m[0][1] = a[l];        node[rt].mat.m[1][0] = 1;        node[rt].mat.m[1][1] = 0;    }    else    {        m = (l + r) >> 1;        Build(rt << 1, l, m);        Build(rt << 1 | 1, m + 1, r);        node[rt].mat = node[rt << 1 | 1].mat * node[rt << 1].mat;    }}Matrix Query(int rt, int ql, int qr){    int m;    if(ql == node[rt].l && node[rt].r == qr)        return node[rt].mat;    else    {        m = (node[rt].l + node[rt].r) >> 1;        if(qr <= m) return Query(rt << 1, ql, qr);        else if(ql > m) return Query(rt << 1 | 1, ql, qr);        else return Query(rt << 1 | 1, m + 1, qr) * Query(rt << 1, ql, m);    }}int T, n, m, ql, qr;int main(){    scanf("%d", &T);    while(T--)    {        Matrix res, f;        scanf("%d%d",&n, &m);        for(int i = 1; i <= n; i++) scanf("%d", &a[i]);        Build(1, 1, n);        for(int i = 1; i<= m; i++)        { scanf("%d%d", &ql, &qr); if(qr - ql >= 2) {     f.m[0][0] = a[ql + 1];     f.m[1][0] = a[ql];     res = Query(1, ql + 2, qr) * f;     printf("%dn", res.m[0][0]); } else     printf("%dn", a[qr]);        }    }    return 0;}
转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/380093.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号