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;
}