（未完结）“文远知行杯”gdet第十四届竞赛（网络赛共10题，仅整理出6题）(代码片段)

#include <iostream>
#include <cstdio>
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)>(y)?(y):(x))

using namespace std;

const int mod=192600817;
const int maxn=4e4+5;
long long fio[maxn]=0,1;
long long sum[maxn]=0,1;

void make_table()

    for(int i=2; i<maxn; i++)
    
        fio[i]=(fio[i-1]+fio[i-2])%mod;
        sum[i]=(fio[i]*fio[i])% mod+sum[i-1];
        //后面是要求 sum的差值，这里的sum不要求余，可大致看一下也能发现long long 不会溢出
    


int main()

    make_table();
    /*
    for(int i=0; i<10; i++)
    
        printf("%d %I64d %I64d
", i , fio[i], sum[i]);
    
    */
    int Q;
    while(cin>>Q)
    
        while(Q--)
        
            int a,b,c,d;
            cin>>a>>b>>c>>d;
            a=a*4+b+1;
            c=c*4+d+1;
            b=max(a,c);
            d=min(a,c);
            cout<<(sum[b]-sum[d-1])%mod<<endl;
            //取余别忘记了
        
    

    return 0;

 1 #include <iostream>
 2 using namespace std;
 3 
 4 const int maxn=150005;
 5 long long table[maxn];
 6 
 7 bool myJudge(int n)
 8 
 9     while(n!=1&&n!=4)
10     
11         int sum=0;
12         while(n)
13         
14             sum+=(n%10)*(n%10);
15             n/=10;
16         
17         n=sum;
18     
19     return n==1;
20 
21 
22 void make_table()
23 
24     int cnt=0;
25     for(int i=1; cnt<maxn; i++)
26     
27         if(myJudge(i))
28             table[cnt++]=i;
29     
30 
31 
32 int main()
33 
34     make_table();
35     int Q; cin>>Q;
36     while(Q--)
37     
38         int k; cin>>k;
39         cout<<table[k-1]<<endl;
40     
41 
42     return 0;
43 

 1 #include <iostream>
 2 #include <cstdio>
 3 using namespace std;
 4 
 5 int main ()
 6 
 7     int T; cin>>T;
 8     while(T--)
 9     
10         double x[10],y[10];
11         for(int i=1;i<=6;i++)
12             cin>>x[i]>>y[i];
13             
14         for(int i=2;i<=3;i++)
15             x[i]-=x[1],y[i]-=y[1];
16         for(int i=5;i<=6;i++)
17             x[i]-=x[4],y[i]-=y[4];
18         
19         double a,b,c,d;
20         c=(x[5]*x[3]-x[6]*x[2])/(y[2]*x[3]-x[2]*y[3]);
21         a=(x[5]-y[2]*c)/x[2];
22         d=(y[5]*x[3]-y[6]*x[2])/(y[2]*x[3]-x[2]*y[3]);
23         b=(y[5]-y[2]*d)/x[2];
24         
25         int Q; cin>>Q;
26         while(Q--)
27         
28             cin>>x[2]>>y[2];
29             x[2]-=x[1],y[2]-=y[1];
30             printf("%.2f %.2f
",x[2]*a+y[2]*c+x[4],x[2]*b+y[2]*d+y[4]);
31         
32     
33     return 0;
34 

#include <iostream>
#include <cstdio>
using namespace std;

long long mod=1000000007;

long long poww(long long a, long long b)

    long long ans=1,base=a;
    while(b)
    
        if(b&1)
            ans=(ans*base)%mod;
        base=(base*base)%mod;
        b>>=1;
    
    return ans;


int main ()

    long long n;
    while(cin>>n)
    
        cout<<(poww(2,n)*((n-1)%mod)+1)%mod<<endl;
    
    return 0;

#include <iostream>
#include <cstdio>
using namespace std;

const int mod=123456789;

struct Matrix
    long long m[6][6];
;

Matrix Mul(Matrix A, Matrix B)

    Matrix ans;
    for(int i=0; i<6; i++)
        for(int j=0; j<6; j++)
            ans.m[i][j]=0;
    for(int i=0; i<6; i++)
        for(int j=0; j<6; j++)
            for(int k=0; k<6; k++)
                ans.m[i][j]=((A.m[i][k]*B.m[k][j])%mod+ans.m[i][j])%mod;
    return ans;


Matrix Quick_Pow(Matrix A, long long n)

    Matrix ans;
    for(int i=0; i<6; i++)
        for(int j=0; j<6; j++)
        
            if(i==j)
                ans.m[i][j]=1;
            else
                ans.m[i][j]=0;
        

    while(n)
    
        if(n&1)
            ans=Mul(ans,A);
        A=Mul(A, A);
        n>>=1;
    
    return ans;


int main()

    int T; cin>>T;
    while(T--)
    
        long long n; cin>>n;

        Matrix A;
        for(int i=0; i<6; i++)
            for(int j=0; j<6; j++)
                A.m[i][j]=0;

        A.m[0][0]=1; A.m[0][1]=2; A.m[0][2]=1; A.m[0][3]=3; A.m[0][4]=3; A.m[0][5]=1;
        A.m[1][0]=1;
                                  A.m[2][2]=1;
                                  A.m[3][2]=1; A.m[3][3]=1;
                                  A.m[4][2]=1; A.m[4][3]=2; A.m[4][4]=1;
                                  A.m[5][2]=1; A.m[5][3]=3, A.m[5][4]=3; A.m[5][5]=1;

        Matrix ans=Quick_Pow(A, n-2);
        long long Ans=0;
        long long a[6]=2,1,1,2,4,8;

        /*
        for(int i=0; i<6; i++)
            for(int j=0; j<6; j++)
            
                printf("%I64d ",A.m[i][j]);
                if(j==5)
                    cout<<endl;
            
        cout<<endl;
        */


        for(int i=0; i<6; i++)
        
            Ans=(Ans+(a[i]*ans.m[0][i])%mod)%mod;
        
        cout<<Ans<<endl;

        /*
        for(int i=0; i<6; i++)
            for(int j=0; j<6; j++)
            
                printf("%I64d ",ans.m[i][j]);
                if(j==5)
                    cout<<endl;
            
        */
    
    return 0;