bzoj3095--数学题

题目大意：
给定一个长度为n的整数序列x[i]，确定一个二元组(b, k)使得S=Σ(k*i+b- x[i])^2（i∈[0,n-1]）最小

看Claris大神的题解就行了。实际上就是用2次二次函数的性质。

http://www.cnblogs.com/clrs97/p/4703437.html

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
inline char Nc(){
    static char buf[100000],*p1=buf,*p2=buf;
    if(p1==p2){
        p2=(p1=buf)+fread(buf,1,100000,stdin);
        if(p1==p2)return EOF;
    }
    return *p1++;
}
inline void Read(int& x){
    char c=Nc(),b=1;
    for(;c<'0'||c>'9';c=Nc())if(c=='-')b=-1;
    for(x=0;c>='0'&&c<='9';x=x*10+c-48,c=Nc());x*=b;
}
int i,n,x;
double k,b,A,B,C,D;
int main()
{
    Read(n);
    for(i=0;i<n;i++){
        Read(x);
        A+=(double)i*i;B+=(double)i;C+=(double)i*x;D+=x;
    }
    k=(n*C-B*D)/(n*A-B*B);
    b=(D-k*B)/n;
    printf("%.7lf %.7lf",b,k);
    return 0;
}