这是一个什么用也没有的模板
貌似牛顿迭代快如闪电,吊打二分。
虽然二次函数可以直接求根公式,,,但是我只是想练练牛顿迭代。
(要是以后考数学的时候我能直接用这个程序求零点就爽了Hhhh)
#include<bits/stdc++.h> #define ll long long #define newton 100 using namespace std; double a,b,c; double x,y;inline bool zero(double k){return k>-0.00000001&&k<0.00000001; }inline void solve(int t){double k,d;//t==1时是求左零点,否则是右零点 if(t==1) x=-b/(2*a)-1;else x=-b/(2*a)+1;//迭代次数 int tmp=newton;while(tmp--){//算出当前点的切线方程,a*x^2+b*x+c的在(x0,y0)处的斜率就是2*a*x0+b k=x*a*2+b,d=y-k*x;//求切线与x轴交点 x=-d/k;y=x*x*a+x*b+c;} }int main(){scanf("%lf%lf%lf",&a,&b,&c);if(b*b-4*a*c<-0.00000001){puts("-1");return 0;}solve(1);if(zero(x)) x=0;if(zero(y)) y=0;printf("%.6lf %.6lf\n",x,y);if(!zero(b*b-4*a*c)){solve(2);if(zero(x)) x=0;if(zero(y)) y=0;printf("%.6lf %.6lf\n",x,y);}return 0; }
然后这个是求A的k次方根的
#include<bits/stdc++.h> #define ll long long using namespace std; double a,b,c,d; double x,y,A; int k;inline void solve(){x=1.00,y=1-A;int tmp=1000;while(tmp--){printf("%lf %lf\n",x,y);//求一下导 b=(double)k*pow(x,k-1);//确定点和斜率就可以确定切线方程了 d=y-b*x;//然后就可以找到和x轴的交点了 x=-d/b;//再找到y=x^k上横坐标是x的点 y=pow(x,k)-A;} }int main(){//求A的k次方根 scanf("%d%lf",&k,&A);solve();printf("%.6lf\n",x);return 0; }
转载于:https://www.cnblogs.com/JYYHH/p/8334312.html
这是一个什么用也没有的模板相关推荐
- 心怡哇1.0.版本一个女生的网站xinyiwa模板 ,WordPress主题1.0女性资讯模板主题
介绍 心怡哇1.0.版本一个女生的网站xinyiwa模板,WordPress<心怡哇>女性资讯模板主题,本模板专注女性资讯,界面简单大方,栏目分明. 目前是1.1基础版(原1.0版作废). ...
- 一个注册页面的前端模板(html+css+javascript)可自适应屏幕
一个注册页面的前端模板(html+css+javascript)可自适应屏幕,修改样式即可用. 截图 代码 <!DOCTYPE html> <html><head> ...
- pc wap模板 html,快速创建一个pc或者wap站点模板
whatsns模板命名规则移动端模板命名格式:以wap结尾自动当作wap模板处理,如 fronzewap pc端模板以default为例,复制default模板文件夹重命名其他英文名称即可,如:bai ...
- java抽组件_GitHub - ysc/HtmlExtractor: HtmlExtractor是一个Java实现的基于模板的网页结构化信息精准抽取组件。...
HtmlExtractor是一个Java实现的基于模板的网页结构化信息精准抽取组件,本身并不包含爬虫功能,但可被爬虫或其他程序调用以便更精准地对网页结构化信息进行抽取. HtmlExtractor是为 ...
- 一个简单的SHELL脚本模板
一个简单的SHELL脚本模板,适用于Java项目的启动,停止,重启. 如果在windows上编辑,请使用Notepad++ 设置编码格式为UTF-8无BOM编码,脚本的文档格式为UNIX格式,否则会因 ...
- 手把手教你使用nodejs编写一个【使用远程仓库模板,快速创建项目模块】的cli(命令行)
目录 实现步骤 初始化cli项目 项目目录 创建交互式命令 拉取远程仓库代码,读取仓库中的模板 拉取远程仓库代码 ora 终端 loading 读取仓库中的模板 将选择的模板复制写入目标项目 Comm ...
- Mustache模板技术,一个比freemarker轻量级的模板引擎
一.初识Mustache 同样也是看Dropwizard才知道这个东西的,以前一直是知道诸如freemarker这样的模板引擎,这个是头一次听说,但是听周围的朋友说最早这个东西是出自于JS的,Drop ...
- 一个自己写的PHP模板引擎
这是一个自己写的编译型的模板引擎(不包括缓存部分)贴上代码.有问题可以随时跟帖. 单文件版:JTemplate.class.php 1 2 3 4 5 6 7 8 9 10 11 12 13 14 1 ...
- 一个使用多年的Makefile模板
layout: post title: category : linux系统 tags : [Makefile] 许久以前,一直使用自己总结的Makefile模板,这个模板也是基于现有资料整理而成的. ...
最新文章
- linux下锐捷拨号的问题
- python国产_Python勒索软件来袭,国产杀软集体失身
- ML:自然语言处理NLP面试题
- gl.vertexAtteib3f P42 讲数据传给location参数指定的attribute变量
- 巧用1个GPIO控制2个LED显示4种状态
- Yii和ThinkPHP对比心得
- python queue 生产者 消费者_生产者、消费者模型---Queue类
- 记忆模糊知识点5-11
- 深入浅出设计模式之状态模式、代理模式
- 南京应用计算机学校排名,南京软件工程专业比较好的17所大学推荐
- 360篡改html文件打开方式,Chrome浏览器被360导航篡改怎么办
- 【论文笔记】气道树分割:A 3D UNet-Graph Neural Network for Airway Segmentation
- 知云软件登陆但不显示登陆二维码
- 一枚中级网络工程师的工作日常,能引起多少同行的共鸣啊。
- 授权服务是如何颁发授权码和访问令牌的?
- oracle创建用户saler,Linux环境下sqlldr一个csv文件
- “¥” 人民币符号及其他币种键盘快捷键
- 民营中小企业融资难的原因
- UltraISOPE 9.6.2.3059简体中文注册版/单文件版+软碟通
- 渡一教育公开课web前端开发JavaScript精英课学习笔记(三)条件语句,循环语句