Pagini recente » Cod sursa (job #2830329) | Cod sursa (job #557838) | Cod sursa (job #2199968) | Cod sursa (job #2582920) | Cod sursa (job #355138)
Cod sursa(job #355138)
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
#define FIN "bibel.in"
#define FOUT "bibel.out"
#define N 18
#define MAXST 131072
#define oo 1000000000
#define minim(a,b) ((a) > (b)) ? (b) : (a)
int d1[N], d2[N][MAXST], P[N], d[N][N];
pair <short int, short int> v1[N], v2[N];
inline int dist(pair <short int, short int> a, pair <short int, short int> b)
{
return ((b.first - a.first) * (b.first - a.first) + (b.second - a.second) * (b.second - a.second));
}
void solve()
{
int i, st, k;
for (i = 0; i <= 17; ++ i)
{
d1[i] = d2[i][P[v1[0].first] - 1];
//memcpy(d1[i], d2[i], sizeof(d2[i]));
for (k = 0; k < P[17]; ++ k)
d2[i][k] = oo;
}
for (i = 1; i <= v2[0].first; ++ i)
for (k = 1; k <= v2[0].first; ++ k)
d[i][k] = dist(v2[i], v2[k]);
for (i = 0; i < v2[0].first; ++ i)
for (k = 0; k < v1[0].first; ++ k)
d2[i][P[i]] = minim(d2[i][P[i]], d1[k] + dist(v1[k + 1], v2[i + 1]));
for (st = 0; st < P[v2[0].first]; ++ st)
for (i = 0; i < v2[0].first; ++ i)
if (st & P[i])
for (k = 0; k < v2[0].first; ++ k)
if (st & P[k])
d2[i][st] = minim(d2[i][st], d2[k][st - P[i]] + d[k + 1][i + 1]);
}
int main()
{
int i, t = 0, r;
freopen(FIN, "r", stdin);
freopen(FOUT, "w", stdout);
v1[1] = make_pair(0, 0);
v1[0].first = 1;
for (i = 0; i < 18; ++ i)
P[i] = 1 << i;
while (t != 2)
{
scanf("%u", &t);
if (!t) {
++ v2[0].first;
scanf("%u%u", &v2[v2[0].first].first, &v2[v2[0].first].second);
}
if (t == 1)
{
solve();
for (i = 0, r = oo; i < v2[0].first; ++ i)
r = minim(r, d2[i][P[v2[0].first] - 1]);
printf("%u\n", r);
memcpy(v1, v2, sizeof(v2));
v2[0].first = 0;
}
}
}