#line 2 "/home/snow/competitive-programming/competitive-programming-library/snow/template.hpp" #pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include constexpr long long INF = 1LL << 60; constexpr long double PI = 3.141592653589; #define rep(i, n) for (ll i = 0; i < (n); ++i) #define rep1(i, n) for (ll i = 1; i <= (n); ++i) #define rrep(i, n) for (ll i = (n - 1); i >= 0; --i) #define ALL(obj) (obj).begin(), (obj).end() #define RALL(obj) (obj).rbegin(), (obj).rend() #define pb push_back #define to_s to_string #define len(v) (ll) v.size() #define debug(x) cout << #x << ": " << (x) << '\n' using namespace std; typedef long long ll; typedef unsigned long long ull; typedef long double ld; typedef pair P; typedef tuple tpl; template using vec = vector; template using vec2 = vector>; template using vec3 = vector>>; template < typename T > inline bool chmax(T &a, const T &b) { if (a < b) { a = b; return 1; } return 0; } template < typename T > inline bool chmin(T &a, const T &b) { if (b < a) { a = b; return 1; } return 0; } #line 4 "/home/snow/competitive-programming/competitive-programming-library/snow/io/setup.hpp" namespace snow{ struct IoSetup { IoSetup() { std::cin.tie(nullptr); std::ios::sync_with_stdio(false); std::cout << std::fixed << std::setprecision(10); } } iosetup; } #line 5 "/home/snow/competitive-programming/competitive-programming-library/snow/io/helper.hpp" template< typename T1, typename T2 > std::ostream &operator << (std::ostream &os, const std::pair< T1, T2 > &p) { os << p.first << " " << p.second; return os; } template< typename T1, typename T2 > std::istream &operator >> (std::istream &is, std::pair< T1, T2 > &p) { is >> p.first >> p.second; return is; } template< typename T1, typename T2, typename T3 > std::ostream &operator << (std::ostream &os, const std::tuple< T1, T2, T3 > &t) { auto &[a, b, c] = t; os << a << " " << b << " " << c; return os; } template< typename T1, typename T2, typename T3 > std::istream &operator >> (std::istream &is, std::tuple< T1, T2, T3 > &t) { auto &[a, b, c] = t; is >> a >> b >> c; return is; } template< typename T > std::ostream &operator << (std::ostream &os, const std::vector< T > &v){ for (int i = 0; i < (int)v.size(); ++i) { os << v[i] << (i + 1 != v.size() ? " " : ""); } return os; } template< typename T > std::istream &operator >> (std::istream &is, std::vector< T > &v){ for(T &in : v) is >> in; return is; } template< typename T > std::ostream &operator << (std::ostream &os, const std::set< T > &st){ int ct = 0; for(auto& s : st) os << s << (++ct != st.size() ? " " : ""); return os; } template void input(T&... a){ (std::cin >> ... >> a); } void print() { std::cout << '\n'; } template void print(const T& a, const Ts&... b){ std::cout << a; (std::cout << ... << (std::cout << ' ', b)); std::cout << '\n'; } int drop() { std::cout << '\n'; exit(0); } template int drop(const T& a, const Ts&... b){ std::cout << a; (std::cout << ... << (std::cout << ' ', b)); std::cout << '\n'; exit(0); } #line 4 "main.cpp" using namespace snow; // #include "snow/algorithm/binary-search.hpp" // #include "atcoder/modint" // using namespace atcoder; // using mint = modint1000000007; #line 2 "/home/snow/competitive-programming/competitive-programming-library/snow/algorithm/binary-search.hpp" namespace snow { /** * @brief Binary Search * @param EPS Sepcify if using decimal floating point * @return right */ template< typename T, typename F > T binary_search(T left, T right, F f, T EPS = 1) { while(abs(right - left) > EPS){ T mid = (right + left) / 2; if(f(mid)) right = mid; else left = mid; } return right; } } // namespace snow #line 13 "main.cpp" int main() { ll N, M; cin >> N >> M; vec<> B(M); cin >> B; ld rem = B[0] - 1 + N - B[M - 1]; //5.333 auto bi = [&](ld X)->bool { vec f(M); auto dp = [&](auto&& self, ll x)->ld { if(x == M - 1) return 0; ld d = B[x + 1] - B[x]; f[x] = min(self(self, x + 1) + d, X + 1); return f[x]; }; dp(dp, 0); ld e = 0; rep(i, M) e += f[i]; e /= M; return X > e; }; auto ret = snow::binary_search(0, 1e10, bi, 1e-8); ld X = ret; vec f(M); auto dp = [&](auto&& self, ll x)->ld { if(x == M - 1) return 0; ld d = B[x + 1] - B[x]; f[x] = min(self(self, x + 1) + d, X + 1); return f[x]; }; dp(dp, 0); print(f[0] + rem); return 0; }