竞赛课第九周（埃式筛法，矩阵乘法）

最新推荐文章于 2024-07-16 18:35:54 发布

做人求其滴

最新推荐文章于 2024-07-16 18:35:54 发布

阅读量344

点赞数 9

分类专栏：算法设计和竞赛实践文章标签：算法 c++ 矩阵乘法斐波那契数列素数埃式筛法

本文链接：https://blog.csdn.net/CH3CH2CH4/article/details/138227306

版权

算法设计和竞赛实践专栏收录该内容

19 篇文章 1 订阅

订阅专栏

本文介绍了如何使用C++编程语言实现埃氏筛法求解区间[2,n]内的所有素数对，以及利用矩阵快速幂方法计算极大数的Fibonacci数，展示了算法的实现和测试结果。

摘要由CSDN通过智能技术生成

1.埃式筛法:求区间[2, n]内所有的素数对

【参考代码】

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

const int N = 1e5;
vector<int> prime;
bool visit[N];

int main()
{
 	int n;
 	cin>>n;
 	memset(visit, false, sizeof(visit));
 	for(int i=2; i<=sqrt(n); i++)
 	{
 		if(visit[i] != true) //i是素数
 		{
 			for(int j=i*i; j<=n; j += i) //原来是2i。优化后i*i。i=5时，2*5 3*5 4*5已经在前面i=2,3,4筛过了
			{
				visit[j] = true; //i*2,3,4,5...都不是素数
			}
		 }
	}
	//通过visit数组记录素数
	for(int i=2; i<=n; i++)
	{
		if(visit[i] != true) //i是素数
		{
			prime.push_back(i);
		}
	}
	for(int i=0; i<prime.size(); i+=2)
	{
		cout << prime[i] << ' ' << prime[i+1] << endl;
	}
 	return 0;
}

【测试结果】

2.求第一亿个Fibonacci数

【参考代码】

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod = 1000000007;

struct matrix{
	ll a[2][2];
};

matrix start;
matrix result;

//(Fn+1) = (1, 1)^n * (F1)
//(Fn  )   (1, 0)  	  (F0)

//矩阵乘法
matrix matrix_multiplication(matrix matrix1, matrix matrix2)
{
	matrix temporary;
	temporary.a[0][0] = (matrix1.a[0][0] * matrix2.a[0][0] + matrix1.a[1][0] * matrix2.a[0][1]) % mod;
	temporary.a[0][1] = (matrix1.a[0][0] * matrix2.a[0][1] + matrix1.a[0][1] * matrix2.a[1][1]) % mod;
	temporary.a[1][0] = (matrix1.a[0][0] * matrix2.a[1][0] + matrix1.a[0][1] * matrix2.a[1][1]) % mod;
	temporary.a[1][1] = (matrix1.a[1][0] * matrix2.a[0][1] + matrix1.a[1][1] * matrix2.a[1][1]) % mod;
	return temporary;
}

//矩阵快速幂
void fastpow(int n)
{
	while(n)
	{
		if(n & 1) //位运算，最后一位是否为1
			result = matrix_multiplication(start, result);
			
		start = matrix_multiplication(start, start);
		n = n >> 1; //右移一位
	}
}

int main() //计算第1亿个Fibonacci数
{
 	int n = 1e8; //一亿
 	start.a[0][0] = start.a[0][1] = start.a[1][0] = 1, start.a[1][1] = 0;
 	//(1, 0)
 	//(0, 1)
 	result.a[0][0] = result.a[1][1] = 1, result.a[1][0] = result.a[0][1] = 0;
 	fastpow(n);
// 	cout << result.a[0][0] % mod << endl;
// 	cout << result.a[0][1] % mod << endl;
 	cout << result.a[1][0] % mod << endl;
// 	cout << result.a[1][1] % mod << endl;
 	return 0;
}
//100000000
//908460138