poj2118~Firepersons

最新推荐文章于 2018-07-04 17:53:13 发布

Merc_A

最新推荐文章于 2018-07-04 17:53:13 发布

阅读量1k

点赞数

分类专栏：数论

本文链接：https://blog.youkuaiyun.com/Merc_A/article/details/38050751

版权

数论专栏收录该内容

2 篇文章

订阅专栏

本文记录了作者在解决POJ2118题目时的经历，从最初面对仅有200多人解决的难题感到挫败，到通过研究发现自己的思路正确，使用矩阵乘法和快速幂算法。然而，在实现过程中忽视了n<k的情况，并且在矩阵乘法函数中出现错误，导致WA。经过仔细检查，最终发现并修复了缺少模除的步骤，成功解决问题。

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

poj上的题，比赛结束之后学长跟我说这个题在poj上做出来的只有200多个人，，，，顿时有种不想做了的感觉，晚上回来研究了一会发现自己的思路应该是正确的，就是矩阵乘法加快速幂，时间上肯定没问题，就是不知道错在哪，把代码给学长看了一下，学长说要考虑n<k的情况，我一看自己确实漏了这个，加上之后还是WA。。。。这就蛋疼了，然后就仔仔细细看了好几遍代码，终于发现是在写矩阵乘法的函数时少了一步模除。。。。。

#include<stdio.h>
#include<string.h>
int x,k,a[110],b[110],c[110][110],e[110][110];
void f(int c[110][110],int d[110][110])
{
    int dd[110][110];
    for(int i=1;i<=k;i++)
        for(int j=1;j<=k;j++)
    {
        int t=0;
        for(int l=1;l<=k;l++)
        {
           t+=(c[i][l]*d[l][j])%10000;
        }
        dd[i][j]=t%10000;\\就是在这个地方，一开始漏掉了这个模除。。
    }
    for(int i=1;i<=k;i++)
        for(int j=1;j<=k;j++)
            c[i][j]=dd[i][j];
}
void mi(int n)
{
    while(n!=0)
    {

        if(n%2==1)
        f(e,c);
        f(c,c);
        n/=2;
    }
}
int main()
{
    int sum=0;
    while(scanf("%d",&k)!=EOF)
    {
        sum=0;
        if(k==0) break;
        memset(c,0,sizeof(c));
        memset(e,0,sizeof(e));
        for(int i=1;i<=k;i++)
            for(int j=1;j<=k;j++)
            if(i==j) e[i][j]=1;
        for(int i=0;i<k;i++)
        scanf("%d",&a[i]);
        for(int i=1;i<=k;i++)
            {
                scanf("%d",&b[i]);
                c[1][i]=b[i];
            }
            for(int i=2;i<=k;i++)
            {
                for(int j=1;j<=k;j++)
                {
                    if(i-1==j) c[i][j]=1;
                }
            }

        scanf("%d",&x);
        if(x>=k)
        {
            mi(x-k+1);
       for(int i=1;i<=k;i++)
       sum+=(e[1][i]*a[k-i])%10000;
       sum=sum%10000;
      printf("%d\n",sum);
        }
        else
            printf("%d\n",a[x]);
    }
    return 0;
}

加上之后果然就过了，，，阿门！