6.18

1. gcd
   讲解链接：https://blog.youkuaiyun.com/qq_34992845/article/details/53043339
   题目：

• hdu 6545
  
  代码：

#include <cstdio>
#include <iostream>;
#include <stdio.h>
using namespace std;
int n;
int minn = -100000;
int gcd(int x,int y){//来求最大公约数
	return (x == 0)?y:gcd(y%x,x);
	
}
int main(){
	while(~scanf("%d",&n)){
		if(n < 2){
			break; 
		}
		int sum = (1+n)*n/2;
		int ans = minn;
		for(int i = 2;i < sum;i ++){
			if(sum %i == 0){
				ans = i > ans?i:ans;
			}
			
		}
		printf("%d\n",ans);
		
	}
	return 0;
	
} 

• hdu 1575
  wa 代码：

#include<bits/stdc++.h>
using namespace std;
const int mod = 9973;
const int maxn= 10;
int n,k;
int matrix[maxn][maxn];
int ans[maxn][maxn];
int temp[maxn][maxn];
void qpow(int k,int n){
	if(k == 0){//ans = matrix
		for(int i = 0;i < n;i ++){
			for(int j = 0;j < n;j ++){
				ans[i][j] = (i == j)?1:0;
			}
		}
		return; 
		
	}
	else if(k%2 == 1){
		qpow(k-1,n);
		 //做一次乘
		for(int y = 0;y < n;y ++){
			for(int t = 0;t < n;t ++){
				for(int i = 0;i < n;i ++){
					temp[y][t] += ans[y][i]*matrix[i][t]%mod;
				}
			}
		}
		for(int i = 0;i < n;i ++){
			for(int j = 0;j < n;j ++){
				ans[i][j] = temp[i][j];
			}
		}
		memset(temp,0,sizeof(temp));
		return;
	}
	else{
		qpow(k/2,n);
		for(int y = 0;y < n;y ++){
			for(int t = 0;t < n;t ++){
				for(int i = 0;i < n;i ++){
					temp[y][t] += ans[y][i]*ans[i][t]%mod;
				}
			
			}
		}
		
		for(int i = 0;i < n;i ++){
			for(int j = 0;j < n;j ++){
				ans[i][j] = temp[i][j];
			}
		}
		memset(temp,0,sizeof(temp));
		return;
		
	}
}
int main(){
	int t;
	scanf("%d",&t);
	while(t--){
		scanf("%d %d",&n,&k);
		for(int i = 0;i < n;i ++){
			for(int j = 0;j < n;j ++){
				scanf("%d",&matrix[i][j]);
			}
			
		}
		
		qpow(k,n);
		
		int ans_num = 0;
		for(int i = 0;i < n;i ++){
			ans_num += ans[i][i]%mod;
		}
		ans_num = ans_num%mod;
		printf("%d\n",ans_num);
	}
	
}