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

poj 2118 Firepersons

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

poj 2118 Firepersons

#include <cstdio>#include <cstring>#include <cstdlib>#include <iostream>#define SIZE 110#define mod 10000using namespace std;int n,m,a[SIZE],b[SIZE];struct MT{    int x,y;    int mt[SIZE][SIZE];}ans,def;inline MT operator *(MT a,MT b){    MT c;    memset(c.mt,0,sizeof c.mt);    c.x=a.x; c.y=b.y;    for(int i=1;i<=a.x;i++)        for(int j=1;j<=b.y;j++) for(int k=1;k<=b.x;k++)     c.mt[i][j]=(c.mt[i][j]+(a.mt[i][k]%mod)*(b.mt[k][j]%mod))%mod;    return c;}bool read(){    for(int i=0;i<n;i++) scanf("%d",&a[i]);    for(int i=0;i<n;i++) scanf("%d",&b[i]);    scanf("%d",&m);    if(m<=n-1) {printf("%dn",a[m]);return false;}    m=m-(n-1);    for(int i=1;i<=n;i++) ans.mt[i][1]=a[i-1];    ans.x=n; ans.y=1;    memset(def.mt,0,sizeof def.mt);    for(int i=1;i<n;i++) def.mt[i][i+1]=1;    for(int i=1;i<=n;i++) def.mt[n][i]=b[n-i];    def.x=def.y=n;    return true;}void go(){    while(m)    {        if(m&1) ans=def*ans;        def=def*def;        m>>=1;    }    printf("%dn",ans.mt[n][1]);}int main(){    while(scanf("%d",&n),n)    {        if(!read()) continue;        else go();    }    return 0;}
转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/374283.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

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

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