hdu1290 献给杭电五十周年校庆的礼物

http://acm.hdu.edu.cn/showproblem.php?pid=1290

神奇的方法

类似的分割的题目，不管是曲线也好折线也好

在二维空间中，都可以用二次式表达式表达

在三维空间中，都可以用三次表达式表达

自己硬算出前几个，然后用待定系数法求出一个表达式即可

#include<iostream>
#include<cstdio>
#include<cstring>
#include<fstream>
#include<algorithm>
#include<map>
#include<vector>
#include<cmath>
using namespace std;
#define ll long long

int main()
{
    cin.sync_with_stdio(false);
    int n;
    while (cin>>n)
    {
        cout<<(n*n*n+5*n)/6+1<<endl;
    }
}

另外附上列主元高斯消元法求方程组233，拒绝手算

#include <iostream>
#include <fstream>
#include <iomanip>
#include <cmath>
#define MAXN 500
#define MB n
#define PRECISION 0.001
#define DOSWAP
using namespace std;
double Matrix[MAXN][MAXN];
double x[MAXN];
int n;

void ShowMat()
{
    cout<<"Mat:"<<endl;
    for (int i=0;i<n;i++)
    {
        for (int j=0;j<n+1;j++)
        {
            cout<<setw(10)<<Matrix[i][j];
        }
        cout<<endl;
    }
    cout<<endl;
}

void SwapRow(int a,int b)
{
    double tmp;
    for (int i=0;i<=n;i++)
    {
        tmp=Matrix[a][i];
        Matrix[a][i]=Matrix[b][i];
        Matrix[b][i]=tmp;
    }
}
int main()
{
    fstream in("test2.txt");
    in>>n;
    if (!in)
    {
        cout<<"No File Found"<<endl;
        return 0;
    }
    for (int i=0;i<n;i++)
    {
        for (int j=0;j<n+1;j++)
        {
            in>>Matrix[i][j];
        }
    }
    for (int i=0;i<n;i++)
    {
        in>>x[i];
    }
    //RESOLVE
    for (int i=0;i<n;i++)
    {
        #ifdef DOSWAP
        double maxone(0);
        int mmax(i);
        for (int j=i;j<n;j++)
        {
            if (abs(Matrix[j][i])>abs(maxone))
            {
                maxone=abs(Matrix[j][i]);
                mmax=j;
            }
        }
        cout<<i<<mmax<<endl;
        if (i!=mmax) SwapRow(i,mmax);
        #endif // DOSWAP
        for (int j=i+1;j<=n;j++)
            Matrix[i][j]/=Matrix[i][i];
        Matrix[i][i]=1;
        for (int j=i+1;j<n;j++)
        {
            for (int k=i+1;k<=n;k++)
                Matrix[j][k]-=Matrix[i][k]*Matrix[j][i];
            Matrix[j][i]=0;
        }
        ShowMat();
    }

    for (int i=n-1;i>=0;i--)
    {
        double sum(0);
        for (int j=i+1;j<n;j++)
        {
            sum+=Matrix[i][j]*x[j];
        }
        x[i]=Matrix[i][MB]-sum;
    }

    for (int i=0;i<n;i++)
    {
        cout<<"x["<<i<<"]="<<x[i]<<endl;
    }
    return 0;
}