修路（蓝桥杯 2022国赛）（线性dp）

输入

2 2 2
2 1
1 2

输出

5.24

#include <iostream>
#include <bits/stdc++.h>
#include <cstring>
#include <vector>
#include <unordered_map>
#include <queue>
#include <set>
#include <algorithm>
#define x first
#define y second
#define pb emplace_back
#define fu(i,a,b) for(int i=a;i<=b; ++ i)
#define fd(i,a,b) for(int i=a;i>=b;    -- i)
#define endl '\n'
#define ms(x,y) memset(x,y,sizeof x)
#define ios ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
using namespace std;

typedef long long LL;
typedef unsigned long long ULL;
typedef vector<vector<LL>> VVL;
typedef vector<vector<int>> VVI;
typedef vector<LL> VL;
typedef vector<int> VI;
typedef vector<string> VS;
typedef pair<int,int> PII;
typedef vector<PII> VPII;
typedef pair<PII,int> PIII;
typedef pair<double,double> PDD;
typedef pair<double,int> PDI;
typedef pair<char,int> PCI;
typedef pair<string,int> PSI;
typedef pair<int,string> PIS;
typedef pair<int,char> PIC;
typedef pair<LL,LL> PLL;
typedef __int128 i128;
typedef unsigned long long ULL;
const int N =2e3 + 10,M = N * 2,base =400  ,INF = 0x3f3f3f3f,P = 131;
const double eps = 1e-8;
const int mod = 998244353;
const LL LNF=(LL) INF * INF;

double get_dist(int a,int b)
{
    return sqrt((LL)a * a + (LL)b * b);
}
double f[N][N][2];
// i -> 当前走过了 i个上面的 格子   j->当前走过了 j个下面的 格子
// 0-> 当前停留在上面  1->当前停留在下面 
int a[N],b[N];
inline void solve()
{
    int n,m,d;
    cin >> n >> m >> d;
    
    fu(i,1,n) cin >> a[i];
    fu(i,1,m) cin >> b[i];
    
    sort(a+1,a+1+n);
    sort(b+1,b+1+m);
    
    fu(i,1,n)
    f[i][0][0] = a[i], // 只在上面走 不到下面 
    f[i][0][1] = INF; // 无法不经过 下面 就停留在 下面 
    
    fu(i,1,m)
    f[0][i][1]= get_dist(b[1],d) + b[i] - b[1] , // 第一步先走到下面，然后一直在下面走 
    f[0][i][0] = INF;
    
    fu(i,1,n)
    fu(j,1,m)
    {
        f[i][j][0] = min(f[i-1][j][0] + a[i] -a[i-1] ,f[i-1][j][1] + get_dist(a[i]-b[j],d));
//  当前这一步是 在上面走的 ，从下面走上来 b[j] -> a[i]         
        f[i][j][1] = min(f[i][j-1][1] + b[j] - b[j-1],f[i][j-1][0] + get_dist(a[i]-b[j],d));
//  当前这一步是 在下面走的 ，从上面走下来 a[i] -> a[j]        
    }
    // 最后停留在 上面还是 下面 的 花费最小 
    cout << min(f[n][m][0],f[n][m][1]) << endl;
        
}
signed main()
{
//  freopen("1.txt","r",stdin);
    ios
    
    cout << fixed<<setprecision(2);
    int t=1;
//    cin>>t;
    int now = 1;
    while(t -- )
    {
//      cout<<"Case "; 
//      cout<<"Case #"; 
//      cout<<"Scenario #"; 
//      cout<< now ++ <<": ";
//      cout<< now ++ <<": \n";
        solve();
    }


    return 0;
}