글 읽기 - 안녕하세요! 수치해석 뉴튼메소드 관련해서 질문있습니다!

hd96dud 6년 전 0

따로 수치해석 뉴튼법으로 고차원 다항식을 풀려고 하는데요

다음 해를 어떻게 찾아야 할 지가 고민입니다.

예를 들어 -10~ 10 사이의 해 값을 구한다 하면

최저값인 -10의 값을 넣고 해를 찾았으면 다음 해를 찾아야 하는데 거기서 최소값을 어떻게 잡아야할까요?

막히는 부분에서의 다음해를 찾기위해 넘어가는 ans값을 정해야 하는데 식에 따라 값을 정하기가 애매합니다...

부탁드립니다!

뉴튼법에 빠삭하신분 계시면 알려주세요!

#include<cstdio>
#include<stdio.h>
#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;


vector<double> v;
vector<double> a;
int num = 0;

double f(vector<double> &x, double inp) 
{
	double ret = 0;
	int sz = x.size()-1;
	for (int i = 0; i <= sz; i++)
		ret += x[i] * pow(inp, sz - i);
	return ret;
}

double df(vector<double> &x, double inp) 
{
	double ret = 0;
	int sz = x.size()-1;
	for (int i = 0; i < sz; i++)
		ret += x[i]*(sz - i)*pow(inp, sz-i-1);
	return ret;
}

void cal(vector<double> &x, double init, double limit)
{
	if (init > limit)
		return;
	num++;
	double ans = init;
	for (int i = 0; i < 10000; i++)
		ans = ans - (f(x,ans) / df(x,ans));
	a.push_back(ans);
	if (num < v.size() - 1)
		cal(x, ans+3, limit); //<< 막히는 부분
	return;
}


int main()
{
	int c;
	scanf("%d", &c);
	while (c--)
	{
		int n;
		scanf("%d", &n);
		for (int i = 0; i <= n; i++)
		{
			double xx;
			scanf("%lf", &xx);
			v.push_back(xx);
		}
		cal(v, -10, 10);
		sort(a.begin(), a.end());
		int sz = a.size();
		for (int i = 0; i < sz; i++)
			printf("%lf ", a[i]);
		printf("\n");
		a.clear();
		v.clear();
		num = 0;
	}
}

댓글을 작성하려면 로그인해야 합니다.

문제

문제

출처

ICPC

안녕하세요! 수치해석 뉴튼메소드 관련해서 질문있습니다!

hd96dud 6년 전 0

Baekjoon Online Judge

채점 현황

문제

유저 대회

고등학교 대회

출처

대학교 대회

도움말