N诺刷题----数学相关

同模余定理

求大数取模的时候可以用这个定理

注意，这里没有除法

1500

例题，过

#include <bits/stdc++.h>
using namespace std;

int main()
{
	int n;
	while(scanf("%d",&n)!=EOF)
	{
		int cnt=1;
		for(int i=0;i<5;i++)
		{
			cnt=(cnt%3)*n%3;
		}
		cout<<cnt<<endl;
	}
	system("pause");
	return 0;
}

最大公约数 GCD

1180 最简真分数

简单，思路就找出所有组合然后看有没有公约数，即最大公约数为1

#include <iostream>
#include <string>
#include <string.h>
#include<algorithm>
#include <vector>
#include <map>
using namespace std;

//求最简真分数代码
int gcd(int a,int b)
{
	if(b==0)
		return a;
	else
		return gcd(b,a%b);
}
int main()
{
	int n;
	int a[600];
	while(cin>>n)
	{
		for(int i=0;i<n;i++)
			cin>>a[i];
		sort(a,a+n);
		int cnt=0;
		for(int i=0;i<n-1;i++)
			for(int j=i+1;j<n;j++)
				if(gcd(a[i],a[j])==1)
					cnt++;
		cout<<cnt<<endl;
	}
	//system("pause");
	return 0;
}

最小公倍数 LCM

1041 求最大公约数和最小公倍数

简单
利用最大公约数求最小公倍数

#include <iostream>
#include <string>
#include <string.h>
#include<algorithm>
#include <vector>
#include <map>
using namespace std;
int gcd(int a,int b)
{
	if(b==0)
		return a;
	else
		return gcd(b,a%b);
}
int main()
{
	int m,n;
	while(cin>>m>>n)
	{
		int yue=gcd(m,n);
		int bei=m*n/yue;
		cout<<yue<<" "<<bei<<endl;
	}
	//system("pause");
	return 0;
}

1353 最大公约数

简单

1426 最大公约数

简单

素数

1013 素数判定

记录：
sqrt（）函数的使用
平实无华的算法

#include <bits/stdc++.h>
using namespace std;

int main()
{
	int n;
	while(cin>>n)
	{
		if(n==1)
			n++;
		for(int i=n;;i++)
		{
			int flag=0;
			for(int j=2;j<=sqrt(i);j++)
				if(i%j==0)
				{
					flag=1;
					break;
				}
			if(flag==0)
			{
				cout<<i<<endl;
				break;
			}
		}
	}
	//system("pause");
	return 0;
}

1355

素数判定 简单

素数筛选

记录
不是很理解标注部分
用的是素数的倍数不是素数
直接背这个代码构造素数

const int maxn = 1000000 + 5;
int prime[maxn];
void getPrime() {
	memset(prime, 0, sizeof(prime));
	for (int i = 2; i <= maxn; i++) //输入数
	{
		if (!prime[i]) //未被动过，每个数赋值i prime[0]存的是素数的个数
			prime[++prime[0]] = i;

		for (int j = 1; j <= prime[0] && prime[j] * i <= maxn; j++) 
		{
			prime[prime[j] * i] = 1;//素数倍数不是素数
			if (i % prime[j] == 0)	break;//这句话不理解，大概可能是去重，去了这行也可构造出这个数组
		}
	}
}

1102

过

1375

简单

分解质因数

1156

这个的难点在于给的数比较大，无法求出大数内的质数
由于我们的素数筛选只能到 1000000， 对于更大的素因子
我们可以不继续分解， 因为不会存在两个大于 1000000 的素因子， 如果存在， 那么这个数的范
围一定大于题目所给的范围 10^9
另一个较大的素数因子直接除法得到

#include <iostream>
#include <string>
#include <string.h>
#include<algorithm>
#include <vector>
#include <map>
using namespace std;
const int maxn = 1000000 + 5;
int prime[maxn];
void getPrime() {
	memset(prime, 0, sizeof(prime));
	for (int i = 2; i <= maxn; i++) //输入数
	{
		if (!prime[i]) //未被动过，每个数赋值i prime[0]存的是素数的个数
			prime[++prime[0]] = i;

		for (int j = 1; j <= prime[0] && prime[j] * i <= maxn; j++) 
		{
			prime[prime[j] * i] = 1;//素数倍数不是素数
			if (i % prime[j] == 0)	break;//这句话不理解，大概可能是去重，去了这行也可构造出这个数组
		}
	}
}
int main()
{
	int n;
	while(cin>>n)
	{
		getPrime();
		int cnt=0;
		for(int i=1;i<=prime[0];i++)
		{
			while(n%prime[i]==0)
			{
				n=n/prime[i];
				cnt++;	
			}
		}
		if(n>1)
			cnt++;
		cout<<cnt<<endl;
	}
	return 0;
}

1284 整除问题

1、对N!分解质因数，所以分别对1、2、3……N分解质因数，然后用一个数组下标记录该质因数的个数

2、对a分解质因数，然后计算相同质因数下，1步质因数个数除以2步质因数个数的最小值即为所求

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1000;
int prime[2][maxn];
void getPrime() {
	memset(prime, 0, sizeof(prime));
	for (int i = 2; i <= maxn; i++) //输入数
	{
		if (!prime[0][i]) //未被动过，每个数赋值i prime[0]存的是素数的个数
			prime[0][++prime[0][0]] = i;

		for (int j = 1; j <= prime[0][0] && prime[0][j] * i <= maxn; j++) 
		{
			prime[0][prime[0][j] * i] = 1;//素数倍数不是素数
			if (i % prime[0][j] == 0)	break;
		}
	}
}
int main()
{
	int n,a;
	int b[1000]={0};
	while(cin>>n>>a)
	{
		getPrime();//做出N内的质因数
		for(int k=2;k<=n;k++)//N个数
		{
			int num=k;
			//对每个数分解质因数
			for(int i=1;i<=prime[0][0];i++)
			{
				while(num%prime[0][i]==0)
				{
					num=num/prime[0][i];
					prime[1][i]++;//记录质因数个数
				}
			}
		}
		//对a分解质因数
		for(int i=1;i<=prime[0][0];i++)
		{
			while(a%prime[0][i]==0)
			{
				a=a/prime[0][i];
				b[i]++;//记录质因数个数
			}
		}
		int value=1000;
		for(int i=0;i<1000;i++)
		{
			if(b[i]!=0)//找a 的质因数部分
			{
				int temp=prime[1][i]/b[i];
				value=min(value,temp);
			}
		}
		cout<<value<<endl;
	}
	return 0;
}

二分幂

1017

幂如果太大，用 for循环复杂度为N
所以考虑将幂二分，1 2 4 8 这样有两倍关系

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

int main()
{
	ll x,n;
	while(cin>>x>>n)
	{
		ll res=1;
		while(n>0)
		{
			if(n&1)
				res=res*x%233333;
			x=x*x%233333;
			n=n>>1;
		}
		cout<<res<<endl;
	}
	return 0;
}

数学规律类型

让我们去找有多少种情况的题目， 第一反应就是他是有规律的， 它的规律
可能简单， 也可能很难。 这时候我们先推出它的前几项。
可借助http://oeis.org/算出
两个规律：
斐波那契数列：
卡特兰数：

高精度问题

1474 大整数加法

1475 大整数乘法

记录：
这道题另一个思路是用数组
先乘法在进位

大整数因子

就是大整数除法
简单