第十七天题解

A:

struct mat{
    int m[maxn][maxn];
}unit;//矩阵的数据结构

**重载矩阵*强调内容*乘法**
mat operator * (mat a,mat b)
{
    mat ret;
    ll x;
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
    {
        x=0;
        for(int k=0;k<n;k++)
            x+=mod((ll)a.m[i][k]*b.m[k][j]);
        ret.m[i][j]=mod(x);
    }
    return ret;
}

#include<iostream>
#include<algorithm>
using namespace std;
#define mod(x) ((x)%9973)
typedef long long ll;
const int maxn =11;
int n;
struct mat{
    int m[maxn][maxn];
}unit;
//矩阵乘法
mat operator * (mat a,mat b)
{
    mat ret;
    ll x;
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
    {
        x=0;
        for(int k=0;k<n;k++)
            x+=mod((ll)a.m[i][k]*b.m[k][j]);
        ret.m[i][j]=mod(x);
    }
    return ret;
}
void init_unit()
{
    for(int i=0;i<maxn;i++)
        unit.m[i][i]=1;
    return ;
}
mat pow_mat(mat a,ll n)
{
    mat ret=unit;
    while(n)
    {
        if(n&1)
        {
            ret= ret*a;
        }
        n>>=1;
        a=a*a;
    }
    return ret;
}
int main()
{
    ios::sync_with_stdio(false);
    int x;
    int T;
    cin>>T;
    init_unit();
    while(T--&&cin>>n>>x)
    {
        mat a;
        for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
        cin>>a.m[i][j];
        a=pow_mat(a,x);
        ll ans=0;
        for(int i=0;i<n;i++)
            ans+=a.m[i][i];
        ans%=9973;
        cout<<ans<<endl;
    }
    return 0;
}

B:
题意：给定n和m n表示1……n的数字排列。求第m个排列。

STL中的神器：next_permutation（后一个）和prev_permutation（前一个）函数

依照STL文档的描写叙述，next_permutation函数将按字母表顺序生成给定序列的下一个较大的序列。直到整个序列为减序为止。prev_permutation函数与之相反。是生成给定序列的上一个较小的序列。

二者原理同样。仅遍例顺序相反

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
using namespace std;
#define N 1047
int num[N];
int main()
{
	int n , m , i ,k;
	while(~scanf("%d%d",&n,&m))
	{
		memset(num,0,sizeof(num));
		for(i = 0 ; i < n ; i++)
		{
			num[i] = i+1;
		}
		for(i = 1 ; i < m ; i++)
		{
			next_permutation(num,num+n);
		}
		for(i = 0 ; i < n-1 ; i++)
			printf("%d ",num[i]);
		printf("%d\n",num[n-1]);
	}
	return 0;
}

C:
思路：简单数论问题
gcd(a,b) 为a,b最大公因数
最小公倍数=a*b/gcd(a,b)
gcd(a,b)最简单的写法就是用递归写，好理解也很好记

#include<stdio.h>
int Gcd(int a, int b){
    return b == 0? a: Gcd(b, a%b);
}
int main(){
    int x,y;
    while(scanf("%d%d",&x,&y)!=EOF) printf("%d\n",x*y/Gcd(x,y));
} 

D:

#include "cstring"
#include "iostream"
#include "string.h"
#include "cstdio"
using namespace std;
int ans = 0;
int n = 0;



int main()
{
    while (scanf("%d", &n) != EOF)
    {
        long long a[50];
        a[1] = 3;
        a[2] = 8;
        if (n < 3)
        {
            printf("%lld\n", a[n]);
            continue;
        }
        for (int cnt = 3; cnt <= n; cnt++)
        {
            a[cnt] = 2 * (a[cnt - 1] + a[cnt - 2]);
        }
        printf("%lld\n", a[n]);
    }
}