题解 2017-2018 ACM-ICPC East Central North America Regional Contest (ECNA 2017)
A - Abstract Art
问n个多边形的面积并是多少。
直接粘板子过了…
#define others
#ifdef poj
#include return 0;
return x>0? 1:-1;
}
class Polygon {
private:
int i;
double s;
public:
int n;
point p[maxp];
point& operator[] (int idx) {
return p[idx];
}
void input() {
for (i=0; iscanf("%lf %lf",&p[i].x,&p[i].y);
p[n]=p[0];
}
double Area() {
for (s=0,i=0; i1];
return s/2;
}
};
PDI s[maxN*maxp*2];
Polygon P[maxN];
double S,ts;
int N;
inline double seg(point o,point a,point b) {
if (cmp(b.x-a.x)==0) return (o.y-a.y)/(b.y-a.y);
return (o.x-a.x)/(b.x-a.x);
}
double PolygonUnion() {
int M,c1,c2;
double s1,s2,ret=0;
for (int i=0; ifor (int ii=0; ii0;
s[M++]=mp(0.00,0);
s[M++]=mp(1.00,0);
for (int j=0; jif (i!=j)
for (int jj=0; jj1],P[j][jj]));
c2=cmp(cross(P[i][ii],P[i][ii+1],P[j][jj+1]));
if (c1==0 && c2==0) {
if (((P[i][ii+1]-P[i][ii])^(P[j][jj+1]-P[j][jj]))>0 && i>j) {
s[M++]=mp(seg(P[j][jj],P[i][ii],P[i][ii+1]),1);
s[M++]=mp(seg(P[j][jj+1],P[i][ii],P[i][ii+1]),-1);
}
} else {
s1=cross(P[j][jj],P[j][jj+1],P[i][ii]);
s2=cross(P[j][jj],P[j][jj+1],P[i][ii+1]);
if (c1>=0 && c2<0) s[M++]=mp(s1/(s1-s2),1);
else if (c1<0 && c2>=0) s[M++]=mp(s1/(s1-s2),-1);
}
}
sort(s,s+M);
double pre=min(max(s[0].x,0.0),1.0),now;
double sum=0;
int cov=s[0].y;
for (int j=1; j0.0),1.0);
if (!cov) sum+=now-pre;
cov+=s[j].y;
pre=now;
}
ret+=P[i][ii]*P[i][ii+1]*sum;
}
}
return ret/2;
}
void solve() {
scanf("%d",&N);
for (int i=0; iscanf("%d", &P[i].n);
P[i].input();
ts=P[i].Area();
if (cmp(ts<0)) {
reverse(P[i].p,P[i].p+P[i].n);
P[i][P[i].n]=P[i][0];
ts=-ts;
}
S+=ts;
}
printf("%.9f %.9f\n",S, PolygonUnion());
}
};
int main() {
// file();
Solver::solve();
return 0;
} B - Craters
求一个能保住所有圆的凸包的周长。
将每个圆一万等分成点，然后围一个凸包即可。
//#define backup
#define others
#ifdef poj
#include
// for(int i = 0; i < cnt; i++)
// cout<
double ss = 0, ss2 = 0;
// res[cnt] = res[0];
res[cnt] = res[0];
for(int i = 1; i <= cnt; i++) {
ss += dis(res[i], res[i-1]);
}
// ss = fabs(ss);
printf("%.6f\n", ss);
}
}
int main() {
// file();
solver::solve();
return 0;
} C - DRM Messages
模拟题
//#define backup
#define others
#ifdef poj
#include 2; i++) {
int v = ((str[i]-'A'+v1) + ((str[i+str.size()/2]-'A'+v2))%26+26);
putchar(v%26+'A');
}
}
}
int main() {
// file();
solver::solve();
return 0;
} D - Game of Throwns
拿栈模拟下。
#include return 0;
}
E - Is-A? Has-A? Who Knowz-A?
共4种关系：
例如： a is b ， b is c ， 那么 a is c 。
1. is 和 is ，还是组成 is 。
2. has 和 has ，还是组成 has 。
3. is 和 has ，组成 has 。
4. has 和 is ，组成 has 。
把 is 和 has 看做边的关系，那么过 is 的就是保留原有相对关系，过 has 的就是一定变成 has 。
所以把每个字符串标号成节点，对每个节点进行 bfs ，从而对点与点之间的关系进行修改。
//#define backup
#define others
#ifdef poj
#include F - Keeping On Track
我们考虑枚举每个点作为被删除点，然后来统计第一问的答案。这个过程可以用一遍dfs，一边统计，一遍更新答案。
第二问就只要把第一问的基础上，贪心的把剩下的最大的两块合起来即可。和基本不等式的思想挺像的。
//#define backup
#define others
#ifdef poj
#include
res += ss * V[i];
ss += V[i];
}
if(res > maxv) {
maxv = res;
id = u;
}
}
void dfs2(int u, int fa) {
sz[u] = 1;
int sum = 0;
for(int i = 0; i < G[u].size(); i++) {
int v = G[u][i];
if(v == fa) continue;
dfs2(v, u);
sz[u] += sz[v];
}
}
void solve() {
scanf("%d", &n);
for(int i = 1; i <= n; i++) {
int u, v;
scanf("%d%d", &u, &v);
G[u].push_back(v);
G[v].push_back(u);
}
dfs(0, -1);
dfs2(id, -1);
LL ss = 0;
vectorelse {
sort(all(V));
V[V.size()-2] += V[V.size()-1];
LL ss = 0, res = 0;
for(int i = 0; i < V.size() - 1; i ++) {
res += ss * V[i];
ss += V[i];
}
cout<" "<int main() {
// file();
solver::solve();
return 0;
} G - A Question of Ingestion
:读错题目，成了可以休息一小时从而使得胃容量变回一小时前的状态，写了个 n3 的DP。
令dp[i][j]代表当前准备吃第i个部分，已经连续吃了j次的答案是多少。
枚举休息多长时间转移即可。
//#define backup
#define others
#ifdef poj
#include int main() {
// file();
solver::solve();
return 0;
} H - Sheba’s Amoebas
说了是简单环了，直接dfs或者并查集都可以。
//#define backup
#define others
#ifdef poj
#include int main() {
// file();
solver::solve();
return 0;
} I - Twenty Four, Again
先dfs出交换的结果和交换的代价，然后枚举运算符。
之后手动讨论出所有括号的情况….
用逆波兰表达式计算即可…
真是蠢到不行
#define others
#ifdef poj
#include
char s[111];
for(int i = 0; i <= str.size(); i++) s[i] = str[i];
stack<int> sv;
stack<char> sp;
char c;
int k = 0, flag = 1;
int x, y;
sp.push('\0');
c = s[k];
while (flag) {
if (c >= '0'&&c <= '9') {
sv.push(toNum(s, k));
}
else if (c == '\0'&& sp.top() == '\0')
flag = 0;
else if (c == '(' || (priority(c) > priority(sp.top())))
sp.push(c), k++;
else if (c == ')'&& sp.top() == '(')
sp.pop(), k++;
else if (priority(c) <= priority(sp.top())) {
x = sv.top();
sv.pop();
if(sv.empty()) y = 0;
else {
y = sv.top();
sv.pop();
}
c = sp.top();
sp.pop();
switch (c) {
case '+':y = x + y; break;
case '-':y = y - x; break;
case '*':y = x*y; break;
case '/':
if(x == 0) return 1e9;
if(y % x == 0)
y = y / x;
else
return 1e9;
break;
}
sv.push(y);
}
c = s[k];
}
return sv.top();
}
string x = "+-*/";
map, int> mp;
vector str;
void dfs(int pos, int x, int dep) {
if(dep == x) {
if(mp.count(str)) mp[str] = min(mp[str], dep);
else mp[str] = dep;
return ;
}
if(pos == str.size() - 1) return ;
swap(str[pos], str[pos + 1]);
dfs(0, x + 1, dep);
swap(str[pos], str[pos + 1]);
dfs(pos + 1, x, dep);
}
vector G;
void solve() {
for(int i = 0; i < 4; i++) {
string tmp;
cin >> tmp;
str.push_back(tmp);
}
int res = 1e9;
for(int i = 0; i <= 6; i++)
dfs(0, 0, i);
G.resize(4);
for(auto it : mp) {
int ans = 1e9;
// for(int i = 0; i < G.size(); i++) cout<
// puts("");
// cout<
for(int i = 0; i < 4; i++) {
G[i] = it.first[i];
}
for(int i = 0; i < 4; i++)
for(int j = 0; j < 4; j++)
for(int k = 0; k < 4; k++) {
string str;
//00
str.clear();
str += G[0]; str += x[i]; str += G[1]; str += x[j]; str += G[2]; str += x[k]; str += G[3];
if(cal(str) == 24) ans = min(ans, 0);
//11
str.clear();
str += '('; str += G[0]; str += x[i]; str += G[1]; str += ')'; str += x[j]; str += G[2]; str += x[k]; str += G[3];
if(cal(str) == 24) ans = min(ans, 1);
//22
str.clear();
str += '('; str += G[0]; str += x[i]; str += G[1]; str += x[j]; str += G[2]; str += ')'; str += x[k]; str += G[3];
if(cal(str) == 24) ans = min(ans, 1);
//33
str.clear();
str += '('; str += '('; str += G[0]; str += x[i]; str += G[1]; str += ')'; str += x[j]; str += G[2]; str += ')'; str += x[k]; str += G[3];
if(cal(str) == 24) ans = min(ans, 2);
//44
str.clear();
str += '('; str += G[0]; str += x[i]; str += '('; str += G[1]; str += x[j]; str += G[2]; str += ')'; str += ')'; str += x[k]; str += G[3];
if(cal(str) == 24) ans = min(ans, 2);
//55
str.clear();
str += '('; str += G[0]; str += x[i]; str += G[1]; str += ')'; str += x[j]; str += '('; str += G[2]; str += x[k]; str += G[3]; str += ')';
if(cal(str) == 24) ans = min(ans, 2);
//66
str.clear();
str += G[0]; str += x[i];str += '('; str += G[1]; str += x[j]; str += G[2];str += ')'; str += x[k]; str += G[3];
if(cal(str) == 24) ans = min(ans, 1);
//77
str.clear();
str += G[0]; str += x[i];str += '('; str += G[1]; str += x[j]; str += G[2]; str += x[k]; str += G[3]; str += ')';
if(cal(str) == 24) ans = min(ans, 1);
//88
str.clear();
str += G[0]; str += x[i];str += '('; str += '('; str += G[1]; str += x[j]; str += G[2]; str+=')'; str += x[k]; str += G[3]; str += ')';
if(cal(str) == 24) ans = min(ans, 2);
//99
str.clear();
str += G[0]; str += x[i];str += '('; str += G[1]; str += x[j];str += '('; str += G[2]; str += x[k]; str += G[3]; str += ')';str+=')';
if(cal(str) == 24) ans = min(ans, 2);
//10
str.clear();
str += G[0]; str += x[i]; str += G[1]; str += x[j]; str += '('; str += G[2]; str += x[k]; str += G[3]; str += ')';
if(cal(str) == 24) ans = min(ans, 1);
}
ans += it.second * 2;
res = min(res, ans);
};
if(res == 1e9) puts("impossible");
else cout<int main() {
// file();
Solver::solve();
return 0;
} J - Workout for a Dumbbell
比较繁琐的模拟，要点在于每次到了一个新的器材的位置，根据时间来计算出当前这个器材上有没有人在用。
每次讨论完都要把时间推到当前人准备运动的状态。
//#define backup
#define others
#ifdef poj
#include 10];
}
}
int main() {
// file();
solver::solve();
return 0;
}