//include
//------------------------------------------
#include <vector>
#include <list>
#include <map>
#include <set>
#include <deque>
#include <stack>
#include <bitset>
#include <algorithm>
#include <functional>
#include <numeric>
#include <utility>
#include <sstream>
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <string>
#include <cstring>
#include <ctime>

#include <stdexcept>

using namespace std;

//conversion
//------------------------------------------
inline int toInt(string s) { int v; istringstream sin(s); sin >> v; return v; }
template<class T> inline string toString(T x) { ostringstream sout; sout << x; return sout.str(); }

//debug
#define DUMP(x)  cerr << #x << " = " << (x) << endl
#define DEBUG(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ")" << " " << __FILE__ << endl

// type alias
using ll = long long;
using ull = unsigned long long;

// constant value
const ll MOD = 1000000007LL;

template<class ll>
inline ll getCeilExponent(ll a)
{
	if (a <= 0)
	{
		throw invalid_argument("Vaule equal to or less than 0 should not be valid");
	}

	ll index = 0;
	--a;
	for (; a != 0; a >>= 1)
	{
		++index;
	}
	return index;
}

template<class ll>
inline ll getFloorExponent(ll a)
{
	if (a <= 0)
	{
		throw invalid_argument("Vaule equal to or less than 0 should not be valid");
	}

	ll index = 0;
	a >>= 1;
	for (; a != 0; a >>= 1)
	{
		++index;
	}
	return index;
}


map<ll, ull> getPrimeFactorization(ll num)
{
	map<ll, ull> primeToCount;

	int rootNum = sqrt(num);
	vector<bool> isCompositeNumber(rootNum + 1, false);
	for (int p = 2; num != 1 && p <= rootNum; ++p)
	{
		if (!isCompositeNumber[p])
		{
			for (int i = p * 2; i <= rootNum; i += p)
			{
				isCompositeNumber[i] = true;
			}

			int count = 0;
			while (num % p == 0)
			{
				++count;
				num /= p;
			}
			if (count >= 1)
			{
				primeToCount.insert(make_pair(p, count));
			}
		}
	}
	if (num != 1)
	{
		primeToCount.insert(make_pair(num, 1));
	}

	return primeToCount;
}

// Computational Complexity: O(log(max(a, b))) 
template<class ll>
inline ll getGCD(ll a, ll b)
{
	if (b == 0)
	{
		return a;
	}
	return getGCD(b, a % b);
}

// Computational Complexity: O(log(max(a, b))) 
template<class ll>
inline pair<ll, ll> getBezoutsIdentitySoulution(ll a, ll b)
{
	if (b == 0)
	{
		return { 1, 0 };
	}
	auto sol = getBezoutsIdentitySoulution(b, a % b);
	return { sol.second, sol.first - (a / b) * sol.second };
}

inline ll getModValue(ll a, ll m)
{
	return (a % m + m) % m;
}

// Computational Complexity: O(log(max(a, n))) 
inline ll getModInverse(ll a, ll m)
{
	auto sol = getBezoutsIdentitySoulution(a, m);
	if (a * sol.first + m * sol.second != 1)
	{
		return -1;
	}
	return getModValue(sol.first, m);
}

ll getModCombinatorialNumber(ll n, ll r, ll m)
{
	ll res = 1;
	for (; n > 0; --n, --r)
	{
		res = (((res * n) % m) * getModInverse(n, m)) % m;
	}
	return getModValue(res, m);
}

ll getModPower(ll base, unsigned t, ll m)
{
	if (t == 0)
	{
		return 1LL;
	}

	auto res = getModPower(base, t >> 1, m);
	res = getModValue(res * res, m);
	if (t % 2 == 1)
	{
		res = getModValue(res * base, m);
	}
	return res;
}

int main()
{
	cin.tie(0);
	ios::sync_with_stdio(false);

	int n;
	cin >> n;

	vector<ll> as(n);
	vector<ll> sums(n);
	ll sum = 0;
	for (int i = 0; i < n; ++i)
	{
		cin >> as[i];
		sum += as[i];
		sums[i] = sum;
	}

	int set1End = 1;
	int set2End = 2;
	int set3End = 3;
	ll res = 1000000000000000LL;
	for (; set2End + 1 < n; ++set2End)
	{
		ll minLeftDiff = abs(sums[set2End - 1] - sums[set1End - 1] * 2);
		for (; set1End + 1 < set2End; ++set1End)
		{
			ll currLeftDiff = abs(sums[set2End - 1] - sums[set1End] * 2);
			if (minLeftDiff < currLeftDiff)
			{
				break;
			}
			minLeftDiff = currLeftDiff;
		}

		if (set3End == set2End)
		{
			set3End = set2End + 1;
		}
		ll minRightDiff = abs(sums[n - 1] - sums[set3End - 1] * 2 + sums[set2End - 1]);
		for (; set3End + 1 < n; ++set3End)
		{
			ll currRightDiff = abs(sums[n - 1] - sums[set3End] * 2 + sums[set2End - 1]);
			if (minRightDiff < currRightDiff)
			{
				break;
			}
			minRightDiff = currRightDiff;
		}

		vector<ll> cand{
			sums[n - 1] - sums[set3End - 1],
			sums[set3End - 1] - sums[set2End - 1],
			sums[set2End - 1] - sums[set1End - 1],
			sums[set1End - 1]
		};
		res = min(res, *max_element(cand.begin(), cand.end()) - *min_element(cand.begin(), cand.end()));
	}

	cout << res;

	return 0;
}




Measured memory_bytes: 3328