我们来自五湖四海,不为别的,只因有共同的爱好,为中国互联网发展出一分力!
领航时时彩重庆版

时时彩组选技巧:hdu 1071 The area 高斯消元求二次函数+辛普森积分

2014年07月03日17:40 阅读: 31314 次

领航时时彩重庆版 www.9nwl5.cn 构造系数矩阵,高斯消元求解二次函数,然后两点式求直线函数,带入辛普森积分法无脑AC。。。


?
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
#include<cstdio>
#include<queue>
#include
#include<cstring>
#include<vector>
#include<cmath>
using namespace std;
struct node
{
    double x,y;
}p[4];
double g[10][10];
double f1(double x) //二次函数
{
    return g[0][3]*x*x+g[1][3]*x+g[2][3];
}
double f2(double x) //直线两点式
{
    double y2=p[2].y,x2=p[2].x;
    double y1=p[1].y,x1=p[1].x;
   return (x-x1)/(x2-x1)*(y2-y1)+y1;
}
double f(double x)
{
    return f1(x)-f2(x);
}
double simpson(double a,double b,int n)
{
    double h=(b-a)/n;
    double ans=f(a)+f(b);
    for(int i=1;i<n;i+=2) ans+="4*f(a+i*h);" for(int="" i="2;i<n;i+=2)" return="" ans*h="" 3;="" }="" void="" gauss(int="" cnt)="" {="" int="" i,="" j,="" k;="" double="" tmp,="" big;="" for="" (i="0;" <="" cnt;="" i++)="" (big="0," j="i;" j++)="" if="" (fabs(g[j][i])=""> big)
             {
                 big = fabs(g[j][i]);
                 k = j;
             }
         }
         if (k != i)
         {
             for (j = 0; j <= cnt; j++)
                 swap(g[i][j], g[k][j]);
         }
         for (j = i + 1; j < cnt; j++)
         {
             if (g[j][i])
             {
                 tmp = -g[j][i] / g[i][i];
                 for (k = i; k <= cnt; k++)
                     g[j][k] += tmp * g[i][k];
             }
         }
     }
     for (i = cnt - 1; i >= 0; i--)
     {
         for (j = i + 1; j < cnt; j++)
             g[i][cnt]-=g[j][cnt]*g[i][j];
         g[i][cnt]/=g[i][i];
     }
}
int main()
{
    int cas;
    scanf("%d",&cas);
    while(cas--)
    {
        for(int i=0;i<3;i++) scanf("%lf%lf",&p[i].x,&p[i].y);
        for(int i=0;i<3;i++)
        {
            g[i][0]=p[i].x*p[i].x;  //ax^2
            g[i][1]=p[i].x;         //bx
            g[i][2]=1;              //c
            g[i][3]=p[i].y;         //y
        }
        Gauss(3);
     //   printf("%lfx^2+%lfx+%lf\n",g[0][3],g[1][3],g[2][3]);
        printf("%.2lf\n",simpson(p[1].x,p[2].x,1000));
    }
    return 0;
}
</n;i+=2)></cmath></vector></cstring></algorithm></queue></cstdio>

分享到: 更多
©2019 安全焦点 版权所有.

  • 这些“难民”大部分好吃懒做,无一技之长,犯罪率又高,你要接回去就噹祖宗供着吧。[哈哈] 2019-04-23
  • 铜梁这位30岁小伙抱上吉他骑着摩托去北极过生! 2019-04-23
  • 2018年世界杯防骗宝典!拒绝和骗子一起狂欢! 2019-04-16
  • [微笑]其实很简单就能破这个局:立法禁止通过房地产二次交易获利,炒房就会被杜绝,炒房一旦被杜绝,房价就会受正常供需关系影响波动在合理范围内。 2019-04-16
  • 中国人打仗?被逗了。独生子一代、捧在手里长大的90后00后、站在甲板上飒爽英姿、硝烟战火血肉横飞的场景、会失去行动能力的。这不是耸人听闻。 2019-04-14
  • 乌鲁木齐市水磨沟区开建两座立体停车库 2019-04-14
  • 吕洪蕾:网络信息技术在干部工作中的运用研究 2019-04-04
  • 聚焦中央经济工作会议 2019-03-30
  • 沈杰:在北京打拼的“90后”台湾律师 2019-03-30
  • 法国总统马克龙首次访华 法国居民期待成果 2019-03-25
  • 【新媒体矩阵】长城编小厨 2019-03-25
  • 1949年刚刚进入北平的人民日报校对科夜班工作场面 2019-03-22
  • 江西省“放管服”改革再出实招 2019-03-21
  • 让“毒跑道”绝迹 内蒙古各方联手共同监管校园跑道 2019-03-20
  • 宁夏石嘴山支队隆重举行退伍老兵欢送仪式 2019-03-19
  • 822| 469| 474| 498| 931| 106| 673| 249| 853| 670|