牛顿法求方程近似解
牛顿法是一种求解方程近似根的方法。小蒜的老师给他出了一道题目,让他用牛顿法求下面这个方程的近似解。
2x^3 - 4x^2 + 3x - 6 = 0
根据牛顿法的定义,小蒜对 f(x) = 2x^{3} - 4x^{2} + 3x - 6进行了求导,得到了f'(x) = 6x^{2} - 8x + 3
如给出的代码所示,小蒜已经将 f(x) 和f′(x) 写成了两个 C 语言的函数。并且完成一个newton
函数,实现了牛顿法求解的目的,并在main
函数中进行了调用。newton
函数中,x 的初值被设为了 1.5,按照牛顿法进行近似解迭代会到 10^{-6} 的误差范围时停止。你可以运行一下程序,看看小蒜实现的程序的计算结果(带入方程看看他计算得是否正确呢)。
#include <stdio.h>
#include <math.h>
#define EPSILON 1e-6
double f(double x) {
return 2 * pow(x, 3) - 4 * pow(x, 2) + 3 * x - 6;
}
double f_prime(double x) {
return 6 * pow(x, 2) - 8 * x + 3;
}
double h(double x) {
return pow(x, 3) - 4 * pow(x, 2) + 3 * x - 6;
}
double h_prime(double x) {
return 3 * pow(x, 2) - 8 * x + 3;
}
double newton(double (*fp)(double), double (*fp_prime)(double)) {
double x = 1.5;
while (fabs(fp(x)) > EPSILON) {
x = x - fp(x) / fp_prime(x);
}
return x;
}
int main() {
printf("%g\n", newton(f,f_prime));
printf("%g\n", newton(h, h_prime));
return 0;
}
二分法求方程近似解
二分法是一种求解方程近似根的方法。对于一个函数 f(x),使用二分法求 f(x) 近似解的时候,我们先设定一个迭代区间(在这个题目上,我们之后给出了的两个初值决定的区间 [−20,20]),区间两端自变量 x 的值对应的 f(x) 值是异号的,之后我们会计算出两端 x 的中点位置 x′ 所对应的 f(x') ,然后更新我们的迭代区间,确保对应的迭代区间的两端 x 的值对应的f(x) 值还会是异号的。
重复这个过程直到我们某一次中点值 x' 对应的 f(x′)<ϵ (题目中可以直接用EPSILON
)就可以将这个 x'作为近似解返回给 main
函数了。
例如:
上面所示的一个迭代过程的第一次的迭代区间是[a1,b1],取中点 b2,然后第二次的迭代区间是 [a1,b2],再取中点a2,然后第三次的迭代区间是[a2,b2],然后取 a3,然后第四次的迭代区间是[a3,b2],再取红色中点 c,我们得到发现f(c) 的值已经小于ϵ,输出 c作为近似解。
在这里,我们将用它实现对形如 px + q = 0 的一元一次方程的求解。
在这里,你完成的程序将被输入两个正整数 p 和 q(你可以认为测评机给出的 0<∣p∣≤1000 且 0<∣q∣≤1000),程序需要用二分法求出 px + q = 0 的近似解。
输入格式
测评机会反复运行你的程序。每次程序运行时,输入为一行,包括一组被空格分隔开的符合描述的正整数 p 和 q。你可以认为输入数据构成的方程 px+q=0 都是有解且解在[−20,20] 的区间内。
输出格式
输出为一行,包括一个数字。为方程 px+q=0 的近似解。请使用四舍五入的方式保留小数点后 4 位小数。
样例输入1
55 9
样例输出1
-0.1636
方法一:
#include <stdio.h>
#include <math.h>
#define EPSILON 1e-7
double bisection(int p, int q, double (*func)(int, int, double));
double f(int p, int q, double x);
int main() {
int p;
int q;
scanf("%d%d", &p, &q);
printf("%.4f\n", bisection(p, q, f));
return 0;
}
double bisection(int p, int q, double (*func)(int, int, double)) {
double a = -20;
double b = 20;
double c;
while (1) {
c = (a + b) / 2;
if (f(p, q, c) == 0)
break;
else if (f(p, q,a) * f( p, q,c) < 0)
b = c;
else
a = c;
}
return c;
}
double f(int p, int q, double x) {
return p * x + q;
}
方法二:
#include <stdio.h>
#include <math.h>
#define EPSILON 1e-7
double bisection(int p, int q, double (*func)(int, int, double));
double f(int p, int q, double x);
int main() {
int p;
int q;
scanf("%d%d", &p, &q);
printf("%.4f\n", bisection(p, q, f));
return 0;
}
double bisection(int p, int q, double (*func)(int, int, double)) {
double left = -20; // 初值设定迭代区间的左端点为 -20
double right = 20; // 初值设定迭代区间的右端点为 20
double x = 0; // x 用于保存每次迭代的中点位置值
double fx = 0; // fx 用于保存 x 对应的 f(x) 的值
while (1) {
x = (left + right) / 2; // 计算迭代区间的中点位置
fx = func(p, q, x); // 计算中点位置 x 对应的 f(x) 的值
if (fabs(fx) < EPSILON) {
break; // f(x) 的值已经小于给定的 EPSILON,结束迭代
} else if (fx * func(p, q, left) < 0) {
right = x; // 中点位置 x 对应的 f(x) 与左端点对应的 f(left) 异号,更新迭代区间的右端点为 x
} else {
left = x; // 中点位置 x 对应的 f(x) 与左端点对应的 f(left) 同号,更新迭代区间的左端点为 x
}
}
return x;
}
double f(int p, int q, double x) {
return p * x + q;
}
斐波那契数列
用 fn 表示新数列的第 n 项,则有:f1=f2=1,f_n =a f_{n-1} + bf_{n-2} (n>2)。
输入格式
输入每行包含 4 个整数 n(1≤n≤100),a(1≤a≤10),b(1≤b≤10),p(1≤p≤2000)。
输出格式
输出fn 对 p 取模的值。
格式说明
输出时每行末尾的多余空格,不影响答案正确性
样例输入
3 1 1 1000
样例输出
2
#include<stdio.h>
#include<math.h>
int Fib(int a, int b, int n,int p){
if(n==1||n==2) return 1;
else return (a * Fib(a, b, n - 1, p) + b * Fib(a, b, n - 2, p)) % p;//注意p取余位置
}
int main(){
int n,a,b,p;
scanf("%d %d %d %d",&n,&a,&b,&p);
int result = Fib(a,b,n,p);
printf("%d\n",result);
return 0;
}
字符串弱等于
如果两字符串 aa, bb 满足下面两个条件之一,我们称 aa 弱等于 bb,或者 bb 弱等于 aa。
- aa 等于 bb;
- aa,bb 的长度相同,并且长度都为偶数,把字符串 aa 串从正中间拆开,拆成 a1,a2,同样把 bb 串拆成 b1,b2。a1 弱等于 b1 并且 a2 弱等于b2,或者 a1 弱等于 b2 并且a2 弱等于b1。
输入格式
输入共两行,每行一个字符串,字符串只包含小写字母。字符串的长度小于2x10^5。
输出格式
如果输入两个字符串是弱等于的,输出YES
,否则输出NO
。
格式说明
输出时每行末尾的多余空格,不影响答案正确性
样例输入
abbaabbabaababab
样例输出
YES
#include <stdio.h>
#include <string.h>
int isWeakEqual(char *str1, char *str2) {
int len1 = strlen(str1);
int len2 = strlen(str2);
// 条件1: 字符串相等
if (strcmp(str1, str2) == 0) {
return 1;
}
// 条件2: 长度相等且都为偶数
if (len1 == len2 && len1 % 2 == 0) {
// 将字符串拆分成两部分
int halfLen = len1 / 2;
char a1[halfLen + 1];
char a2[halfLen + 1];
char b1[halfLen + 1];
char b2[halfLen + 1];
strncpy(a1, str1, halfLen); // 字符串str1的前半部分
strncpy(a2, str1 + halfLen, halfLen); // 字符串str1的后半部分
strncpy(b1, str2, halfLen); // 字符串str2的前半部分
strncpy(b2, str2 + halfLen, halfLen); // 字符串str2的后半部分
a1[halfLen] = '\0';
a2[halfLen] = '\0';
b1[halfLen] = '\0';
b2[halfLen] = '\0';
// a1 弱等于 b1 并且 a2 弱等于 b2,或者 a1 弱等于 b2 并且 a2 弱等于 b1
if ((isWeakEqual(a1, b1) && isWeakEqual(a2, b2)) ||
(isWeakEqual(a1, b2) && isWeakEqual(a2, b1))) {
return 1;
}
}
return 0; // 不满足弱等于条件
}
int main() {
char str1[200001]; // 根据题目要求的字符串长度上限定义字符数组
char str2[200001];
scanf("%s", str1);
scanf("%s", str2);
if (isWeakEqual(str1, str2)) {
printf("YES\n");
} else {
printf("NO\n");
}
return 0;
}
对数方程
#include <stdio.h>
#include <math.h>
double equation(double x, double a) {
return x + log(x) - a;
}
double solveEquation(double a) {
double left = 0.1; // 左边界 (可以根据具体问题调整)
double right = 1e9; // 右边界 (可以根据具体问题调整)
double precision = 1e-4; // 精度
while (right - left > precision) {
double middle = (left + right) / 2;
double result = equation(middle, a);
if (result >= 0) {
right = middle;
} else {
left = middle;
}
}
return right; // 取区间中点作为近似解
}
int main() {
double a;
scanf("%lf", &a);
double solution = solveEquation(a);
printf("%.6lf\n", solution);
return 0;
}