为了防止后效性,应该将两组数做一个sort,这样就可以直接进行dp,方程十分显然。
注意浮点数绝对值要使用fabs
= = swap也可以交换数组 方便了很多。。。
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<stdio.h>
using namespace std;
double a[510],b[510],f[510][510];
int n1,n2;
int main()
{
<span style="white-space:pre"> </span>int T;
<span style="white-space:pre"> </span>cin>>T;
<span style="white-space:pre"> </span>while(T--)
<span style="white-space:pre"> </span>{
<span style="white-space:pre"> </span>memset(f,63,sizeof(f));
<span style="white-space:pre"> </span>f[0][0]=0;
<span style="white-space:pre"> </span>cin>>n1>>n2;
<span style="white-space:pre"> </span>for(int i=1;i<=n1;i++) cin>>a[i];
<span style="white-space:pre"> </span>for(int i=1;i<=n2;i++) cin>>b[i];
<span style="white-space:pre"> </span>if(n1>n2) { swap(n1,n2); swap(a,b); }
<span style="white-space:pre"> </span>sort(a+1,a+n1+1);
<span style="white-space:pre"> </span>sort(b+1,b+n2+1);
<span style="white-space:pre"> </span>for(int i=1;i<=n1;i++)
<span style="white-space:pre"> </span>{
<span style="white-space:pre"> </span>f[i][i]=f[i-1][i-1]+fabs(a[i]-b[i]);
<span style="white-space:pre"> </span>for(int j=i+1;j<=n2;j++)
<span style="white-space:pre"> </span> f[i][j]=min(f[i-1][j-1]+fabs(a[i]-b[j]),f[i][j-1]);
<span style="white-space:pre"> </span>}
<span style="white-space:pre"> </span>printf("%.1lf\n",f[n1][n2]);
<span style="white-space:pre"> </span>}<span style="white-space:pre"> </span>
<span style="white-space:pre"> </span>while(1);
<span style="white-space:pre"> </span>return 0;
}