A Ginkgo Numbers

直接暴力查找(m,n)是否为(p,q)的因子

1

#include   
#include   
#include   
#include   
#include   
#include   
#include

B One-Dimensional Cellular Automaton

一道数论题，直接写TLE无疑，只能YY了一下。各种数组保存。。。写得都晕了，居然1A。

#include   
#include   
#include   
#include   
#include   
#define mem(a,b) memset(a,b,sizeof(a))   
#define FOR(a,b,i) for(i=a;i'9');  
    ret=c-'0';  
    while((c=getchar())>='0'&&c<='9')  
    {  
        ret=ret*10+(c-'0');  
    }  
}  
inline void OT(int a)  
{  
    if(a>=10)  
    {  
        OT(a/10);  
    }  
    putchar(a%10+'0');  
}  
int s[51],f[51][51][51],q[51][51],tx[51];  
int main()  
{  
    int n,m,a,b,c,t,i,j,k,l,p,tt;  
    while(1)  
    {  
        RD(n);  
        RD(m);  
        RD(a);  
        RD(b);  
        RD(c);  
        RD(t);  
        if(n==0&&m==0&&a==0&&b==0&&c==0&&t==0)  
        {  
            break;  
        }  
        For(1,n,i)  
        {  
            RD(s[i]);  
        }  
        if(t==0)//等于0时没有变化，直接输出   
        {  
            FOR(1,n,i)  
            {  
                OT(s[i]);  
                printf(" ");  
            }  
            OT(s[n]);  
            printf("\n");  
        }  
        else  
        {  
            For(1,n,i)//三维数组保存   
            {  
                f[0][i-1][i]=a;  
                f[0][i][i]=b;  
                f[0][i+1][i]=c;  
            }  
            p=0;  
            while(t>1)//转化为2进制进行操作   
            {  
                tx[++p]=t%2;  
                t=t/2;  
            }  
            For(1,p/2,i)  
            {  
                tt=tx[i];  
                tx[i]=tx[p-i+1];  
                tx[p-i+1]=tt;  
            }  
            For(1,p,l)  
            {  
                For(1,n,i)  
                {  
                    For(1,n,j)  
                    {  
                        f[l][i][j]=0;  
                    }  
                }  
                For(1,n,k)  
                {  
                    For(1,n,i)  
                    {  
                        For(1,n,j)  
                        {  
                            f[l][i][j]=(f[l][i][j]+f[l-1][i][k]*f[l-1][k][j])%m;  
                        }  
                    }  
                }  
                if(tx[l]==1)  
                {  
                    For(1,n,i)  
                    {  
                        For(1,n,j)  
                        {  
                            q[i][j]=f[l][i][j];  
                            f[l][i][j]=0;  
                        }  
                    }  
                    For(1,n,k)  
                    {  
                        For(1,n,i)  
                        {  
                            For(1,n,j)  
                            {  
                                f[l][i][j]=(f[l][i][j]+q[i][k]*f[0][k][j])%m;  
                            }  
                        }  
                    }  
                }  
            }  
            mem(q,0);  
            For(1,n,k)  
            {  
                For(1,n,i)  
                {  
                    For(1,n,j)  
                    {  
                        q[i][j]=(q[i][j]+s[k]*f[p][k][j])%m;  
                    }  
                }  
            }  
            FOR(1,n,i)  
            {  
                OT(q[1][i]);  
                printf(" ");  
            }  
            OT(q[1][n]);  
            printf("\n");  
        }  
    }  
    return 0;  
}