#include #define M_PI 3.14159265358979323846 using namespace std; //conversion //------------------------------------------ inline int toInt(string s) { int v; istringstream sin(s); sin >> v; return v; } template inline string toString(T x) { ostringstream sout; sout << x; return sout.str(); } inline int readInt() { int x; scanf("%d", &x); return x; } //typedef //------------------------------------------ typedef vector VI; typedef vector VVI; typedef vector VS; typedef pair PII; typedef pair TIII; typedef long long LL; typedef unsigned long long ULL; typedef vector VLL; typedef vector VVLL; //container util //------------------------------------------ #define ALL(a) (a).begin(),(a).end() #define RALL(a) (a).rbegin(), (a).rend() #define PB push_back #define MP make_pair #define SZ(a) int((a).size()) #define SQ(a) ((a)*(a)) #define EACH(i,c) for(typeof((c).begin()) i=(c).begin(); i!=(c).end(); ++i) #define EXIST(s,e) ((s).find(e)!=(s).end()) #define SORT(c) sort((c).begin(),(c).end()) //repetition //------------------------------------------ #define FOR(i,s,n) for(int i=s;i<(int)n;++i) #define REP(i,n) FOR(i,0,n) #define MOD 1000000007 #define rep(i, a, b) for(int i = a; i < (b); ++i) #define trav(a, x) for(auto& a : x) #define all(x) x.begin(), x.end() #define sz(x) (int)(x).size() typedef long long ll; typedef pair pii; typedef vector vi; const double EPS = 1E-8; #define chmin(x,y) x=min(x,y) #define chmax(x,y) x=max(x,y) const int INF = 100000000; struct Edge { int to, from; ll cost; Edge(int from, int to, ll cost): from(from), to(to), cost(cost) {} }; struct UnionFind { vector data; UnionFind(int size) : data(size, -1) { } bool unionSet(int x, int y) { x = root(x); y = root(y); if (x != y) { if (data[y] < data[x]) swap(x, y); data[x] += data[y]; data[y] = x; } return x != y; } bool findSet(int x, int y) { return root(x) == root(y); } int root(int x) { return data[x] < 0 ? x : data[x] = root(data[x]); } int size(int x) { return -data[root(x)]; } }; typedef vector> AdjList; AdjList graph; ll mod_pow(ll x, ll n, ll mod){ ll res = 1; bool c = false; while(n){ if(n&1) res = res * x; if(res > mod){ c = true; res %= mod; } x = x * x %mod; n >>= 1; } if(c) return mod; return res; } #define SIEVE_SIZE 10000+10 bool sieve[SIEVE_SIZE]; void make_sieve(){ for(int i=0; i vector gauss_jordan(const vector>& A, const vector& b){ int n = A.size(); vector> B(n, vector(n+1)); for(int i=0; i abs(B[pivot][i])) pivot = j; } swap(B[i], B[pivot]); if(abs(B[i][i]) < EPS) return vector(); //解なし for(int j=i+1; j<=n; ++j) B[i][j] /= B[i][i]; for(int j=0; j x(n); for(int i=0; i vec; typedef vector mat; mat mul(mat &A, mat &B) { mat C(A.size(), vec((int)B[0].size())); for(int i=0; i 0) { if(n & 1) B = mul(B, A); A = mul(A, A); n >>= 1; } return B; } bool operator<(const pii& a, const pii& b){ if(a.first == b.first) return a.second < b.second; return a.first < b.first; } ll dp[(1<<15)][15]; vector vv; int n; ll rec(int S, int v){ if(dp[S][v]>=0) return dp[S][v]; if(S == (1 << n) -1){ return dp[S][v] = 0; } ll res = INT_MAX; ll waribiki=0; for(int i=0; i>i & 1) waribiki += vv[i]; } for(int u=0; u> u & 1)){ res = min(res, rec(S | (1 << u), u) + max(0ll, vv[u] - waribiki%1000)); } } return dp[S][v] = res; } int main() { cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(20); cin >> n; vv.resize(n); REP(i,n) cin >> vv[i]; ll ans = 1e11; for(int i=0; i<(1<