ACM 刘春英老师 课程笔记 1

输入未知次

例：输入两个整数为一组，共有未知组输入

c语法

int a, b;
while(scanf("%d %d", &a, &b) != EOF)

while(scanf("%d %d", &a, &b) == 2)

C++语法

while(cin >> a >> b) //>>读到文件尾返回0

数学基础

约数个数定理

对于一个大于1正整数n可以分解质因数：

则n的正约数的个数就是

其中a1、a2、a3…ak是p1、p2、p3，…pk的指数。

最大公因数GCD

辗转相除法（欧几里得算法）

int gcd(int big, int small) { //big和small参数可以互换
	int t;
	while(small != 0) {
		t = big % small;
		big = small;
		small = t;
	}
	return big;
} 

最小公倍数LCM

LCM(A,B) = A*B / GCD(A,B)

LCM：最小公倍数

GCD：最大公因数

注意：由于A*B容易越界，一般写作

LCM(A,B) = A / GCD(A,B) * B

快速幂（二分加速）

int power(int a, int b) { //递归实现 
	int ans;
	if (n == 0) ans = 1;
	else {
		ans = power(a*a, n/2);
		if (n % 2 == 1) ans *= a;
	}
	return ans;
} 

int power(int a, int b) { //非递归实现 循环 
	int ans = 1;
	while(n) {
		if (n % 2) ans *= a;
		a = a * a;
		n = n /2;
	}
	return ans;
}

二分查找

前提——有序数组，用于单调函数

int BiSearch(int a[], int n, int x) {
	int left = 0, right = n-1;
	while(left <= right) { //注意 = 不能少 
		int middle = (left + right) / 2;
		if (a[middle] == x) 
			return middle;
		if (x > a[middle]) 
			left = middle + 1;
		else
			right = middle - 1; 
	}
	return -1;
} 

int BiSearch(int a[], int x, int left, int right) { //递归
	if (left > right) {   //注意不能有 =
		return -1;
	} else {
		int mid = (left + right) / 2;
		if (a[mid] == x) {
			return mid;
		} else if (x > a[mid]) {
			return BiSearch(a, x, mid+1, right);
		} else if (x < a[mid]) {
			return BiSearch(a, x, left, mid-1);
		}
	}
}

二分求解

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

double y;
double l, r, mid;
double f (double x) {
	return 8*pow(x, 4.0) + 7*pow(x, 3.0) + 2*pow(x, 2.0) + 3.0*x + 6; 
}

int main() {
	int t;
	cin >> t;
	for (int i = 0; i < t; i++) {
		cin >> y;
		if (f(0) <= y && y <= f(100)) {
			l = 0;
			r = 100;
			while (r - l > 1e-7){
				mid = (l + r) / 2;
				double ans = f(mid);
				if (ans > y) {
					r = mid - 1e-8;
				} else {
					l = mid + 1e-8;
				}
			}
			double a = (l+r) / 2;   //相互一一定要先把ans求出，在printf中求会出问题
			printf("%.4f\n", a);
		} else {
			cout << "No solution!" << endl;
		}
	}
	return 0;
}

三分查找

前提——数据的凸性

凸性并不要求可导

double ternarySearch(double l, double r) {
  while (r - l >= EPS) { //EPS 精确度
    double m1 = l + (r - l) / 3.0;
    double m2 = r - (r - l) / 3.0;
    if (f(m1) >= f(m2)) {
      r = m2;
    } else {
      l = m1;
    }
  }
  return l;
}

————————————————

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
                        
原文链接：https://blog.csdn.net/nazeniwaresakini/article/details/104220074