1. C++ / Говнокод #17879

    +54

    1. 01
    2. 02
    3. 03
    4. 04
    5. 05
    6. 06
    7. 07
    8. 08
    9. 09
    10. 10
    11. 11
    12. 12
    13. 13
    14. 14
    15. 15
    16. 16
    17. 17
    18. 18
    19. 19
    20. 20
    21. 21
    22. 22
    23. 23
    24. 24
    25. 25
    26. 26
    27. 27
    28. 28
    29. 29
    30. 30
    31. 31
    32. 32
    33. 33
    34. 34
    35. 35
    36. 36
    37. 37
    38. 38
    39. 39
    40. 40
    41. 41
    42. 42
    43. 43
    44. 44
    45. 45
    46. 46
    47. 47
    48. 48
    49. 49
    50. 50
    51. 51
    52. 52
    53. 53
    54. 54
    55. 55
    56. 56
    57. 57
    58. 58
    59. 59
    60. 60
    61. 61
    62. 62
    63. 63
    64. 64
    65. 65
    66. 66
    67. 67
    68. 68
    69. 69
    70. 70
    71. 71
    72. 72
    73. 73
    74. 74
    75. 75
    76. 76
    77. 77
    78. 78
    79. 79
    80. 80
    81. 81
    82. 82
    83. 83
    84. 84
    85. 85
    86. 86
    87. 87
    88. 88
    89. 89
    90. 90
    91. 91
    92. 92
    93. 93
    94. 94
    95. 95
    96. 96
    97. 97
    98. 98
    std::pair<long double, long double> Calculator::leastsquares(const QVector<double> &x, const QVector<double> &yy) const
{
    QVector<double> y = yy;
    for (int i = 0; i < y.size(); ++i) {
        y[i] += 1.0;
    }
    long double A, B, a = 0, b = 0, at = 0, tt = 0, bt = 0, tmp = 0;
    for (int i = 0; i < x.size(); ++i) {
        tmp += x[i] * x[i] * y[i];
    }
    
    a = tmp;
    tmp = 0;
    
    for (int i = 0; i < x.size(); ++i) {
        tmp += y[i] * log(y[i]);
    }
    
    a *= tmp;
    tmp = 0;
    
    for (int i = 0; i < x.size(); ++i) {
        tmp += y[i] * x[i];
    }
    
    at = tmp;
    tmp = 0;
    
    for (int i = 0; i < x.size(); ++i) {
        tmp += y[i] * x[i] * log(y[i]);
    }
    
    at *= tmp;
    tmp = 0;
    
    a -= at;
    
    for (int i = 0; i < x.size(); ++i) {
        tmp += y[i];
    }
    
    at = tmp;
    tmp = 0;
    
    for (int i = 0; i < x.size(); ++i) {
        tmp += x[i] * x[i] * y[i];
    }
    
    tt = at * tmp;
    tmp = 0;
    
    for (int i = 0; i < x.size(); ++i) {
        tmp += x[i] * y[i];
    }
    
    tt -= tmp * tmp;
    tmp = 0;
    
    a /= tt;
    
    A = exp(a);
    
    for (int i = 0; i < x.size(); ++i) {
        tmp += y[i];
    }
    
    b = tmp;
    tmp = 0;
    
    for (int i = 0; i < x.size(); ++i) {
        tmp += x[i] * y[i] * log(y[i]);
    }
    
    b *= tmp;
    tmp = 0;
    
    for (int i = 0; i < x.size(); ++i) {
        tmp += y[i] * x[i];
    }
    
    bt = tmp;
    tmp = 0;
    
    for (int i = 0; i < x.size(); ++i) {
        tmp += y[i] * log(y[i]);
    }
    
    bt *= tmp;
    tmp = 0;
    
    b -= bt;
    
    b /= tt;
    
    B = b;
    
    return std::make_pair(A, B);
}

    Запостил: Abbath, 28 Марта 2015

    Комментарии (0) RSS

    Добавить комментарий

    Ошибка компиляции комментария:
    1. Гости могут высказаться только в понедельник, среду, четверг или воскресение
    ava Я, guest, находясь в здравом уме и твердой памяти, торжественно заявляю:
    А не использовать ли нам bbcode?
    • [b]жирный[/b] — жирный
    • [i]курсив[/i] — курсив
    • [u]подчеркнутый[/u] — подчеркнутый
    • [s]перечеркнутый[/s] — перечеркнутый
    • [blink]мигающий[/blink] — мигающий
    • [color=red]цвет[/color] — цвет (подробнее)
    • [size=20]размер[/size] — размер (подробнее)
    • [code=<language>]some code[/code] (подробнее)
    Проверочный код