线性回归——以一道等差数列的题为例

Shanhj

已于 2022-07-23 21:13:12 修改

阅读量206

点赞数

文章标签：线性回归算法机器学习

于 2022-07-23 21:12:36 首次发布

本文链接：https://blog.csdn.net/Shanhj/article/details/125952802

版权

线性回归等差数列代价计算 C++代码数值优化

关键词由CSDN通过智能技术生成

数论专栏收录该内容

5 篇文章 0 订阅

订阅专栏

对于一组 x_i 和 y_i 计算线性回归方程的公式如下：
在这里插入图片描述
题目大意：给定一个数列，将其变成等差数列，修改的代价是
求最小代价

直接套用上面的线性回归方程公式即可。

#include <bits/stdc++.h>
#define inf 0x3f3f3f3f
const long long N = 3e5 + 5;
using namespace std;

namespace GTI //快读模板
{
    char gc(void)
    {
        const long long S = 1 << 16;
        static char buf[S], *s = buf, *t = buf;
        if (s == t) t = buf + fread(s = buf, 1, S, stdin);
        if (s == t) return EOF;
        return *s++;
    }
    long long gti(void)
    {
        long long a = 0, b = 1, c = gc();
        for (; !isdigit(c); c = gc())
            b ^= (c == '-');
        for (; isdigit(c); c = gc())
            a = a * 10 + c - '0';
        return b ? a : -a;
    }
}
using GTI::gti;

long long num[N];

signed main()
{
    long long T;
    T = gti();
    while (T--)
    {
        long long n = gti();
        long long sxy = 0, sx = 0, sy = 0, sx2 = 0;
        long double b, a, ans = 0;
        for (long long i = 1; i <= n; i++)
        {
            num[i] = gti();
            sxy += i * num[i];
            sx += i;
            sy += num[i];
            sx2 += i * i;
        }
        b = ((long double)n * sxy - (long double)sx * sy) / ((long double)n * sx2 - (long double)sx * sx);
        a = (long double)sy / n - b * sx / n;
        long double ai = a;
        for (long long i = 1; i <= n; i++)
        {
            ai += b;
            ans += (num[i] - ai) * (num[i] - ai);
        }
        printf("%.10lf\n", ans);
    }
    return 0;
}