Uva 10014 - Simple calculations

额,貌似应该线性方程组消元?

表示果断二分法结果...

看成a[n+1]关于a[1]的函数

然后二分,即可

2Y

第一次看错了n...

#include<stdio.h>
#define eps 1e-3
double abs(double x){
	if(x>0)return x;
	return -x; 
	}
double a[3010],c[3010];
double aim;
int n;
double progress(){
	int i;
	for(i=2;i<=n+1;i++)a[i]=2*a[i-1]+2*c[i-1]-a[i-2];
	return a[n+1]-aim;
	}
int main(){
	int t;
	scanf("%d",&t);
	while(t--){
		scanf("%d",&n);
		int i;
		scanf("%lf%lf",&a[0],&aim);
		double l=-1000.0;
		double r=1000.0;
		for(i=1;i<=n;i++)scanf("%lf",&c[i]);
		while(abs(r-l)>eps){
			double m=(l+r)/2;
			double fl,fr,fm;
			a[1]=l;fl=progress();//printf("f(%.2lf)=%.2lf\n",l,fl);
			a[1]=r;fr=progress();//printf("f(%.2lf)=%.2lf\n",r,fr);
			a[1]=m;fm=progress();//printf("f(%.2lf)=%.2lf\n",m,fm);
			if(fr*fm<0)l=m;
			else r=m;
			}
		printf("%.2lf\n",l);
		if(t)printf("\n");
		}
	return 0;
	}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值