2016 icpc沈阳部分题解

hdu5948 hdu5949 hdu5950 hdu5952 hdu5954

A Thickest Burger

​ 对于每对 A B, 输出 $max(2A+B,2B+A)$

B Relative atomic mass

​ 给出一个字符串，只包含CHO，原子质量分别是12，1，16，求出这个字符串总的质量即可

C Recursive sequence

​ 矩阵快速幂，式子是 $a[i]=2a[i-2]+a[i-1]+i^4$，关键在于根据这个式子推导出转移矩阵，然后直接 $logn$ 实现即可

#include<bits/stdc++.h>
using namespace std;

#define int long long
const int mod =  2147493647;

struct matrix{
	int a[8][8];

	matrix(){
		memset(a, 0, sizeof a);
	}
};

matrix mul(const matrix& x, const matrix& y){
	matrix ans;
	for(int i = 0; i < 7; ++i){
		for(int j = 0; j < 7; ++j){
			for(int k = 0; k < 7; ++k){
				ans.a[i][j] += x.a[i][k] * y.a[k][j];
				ans.a[i][j] %= mod;
			}
		}
	}
	return ans;
}

matrix qpow(matrix a, int b){
	matrix ans;
	memcpy(ans.a, a.a, sizeof a.a);
	if(b == 0) return ans;
	b--;
	while(b){
		if(b & 1LL) ans = mul(ans, a);
		a = mul(a, a);
		b >>= 1LL;
	}	
	return ans;
}

signed main(){
	matrix A, B;
	A.a[0][2] = 16;
	A.a[0][3] = 8;
	A.a[0][4] = 4;
	A.a[0][5] = 2; 
	A.a[0][6] = 1;

	B.a[0][0] = B.a[0][1] = 1;
	B.a[1][0] = 2;
	B.a[2][0] = B.a[2][2] = 1;
	B.a[3][0] = B.a[3][2] = 4;
	B.a[3][3] = 1;
	B.a[4][0] = B.a[4][2] = 6;
	B.a[4][3] = 3;
	B.a[4][4] = 1;
	B.a[5][0] = B.a[5][2] = 4;
	B.a[5][3] = 3;
	B.a[5][4] = 2;
	B.a[5][5] = 1;
	for(int i = 0; i < 7; ++i)
		B.a[6][i] = 1;
	B.a[6][1] = 0;

	int t;
	cin >> t;
	while(t--){
		int n, a, b;
		cin >> n >> a >> b;
		matrix C = A, D = B;
		C.a[0][1] = a;
		C.a[0][0] = b;
		if(n == 1){
			cout << a << endl;
			continue;
		}
		else if(n == 2){
			cout << b << endl;
			continue;
		}
		n -= 2;
		D = qpow(D, n);
		C = mul(C, D);
		cout << C.a[0][0] << endl;
	}
	return 0;
}

D Couting Cliques

​ 给出一个无向图，有n个顶点和m条边，求顶点刚好为S的完全子图。
​ 这里因为n只有100，m只有1000，S不超过10，所以直接爆搜，然后剪枝即可（一开始的时候还在苦苦想着有没有更好的想法，大意了）

#include<bits/stdc++.h>
using namespace std;

const int N = 111;

const int M = 1111;
int e[M], ne[M], h[N], idx;

void add(int a, int b){
	e[idx] = b;
	ne[idx] = h[a];
	h[a] = idx++;
}

int vis[N];
int mp[N][N];
int n, m, s;
int ans;
int a[11];
int at;


inline int read() {
    int X = 0; bool flag = 1; char ch = getchar();
    while (ch < '0' || ch>'9') { if (ch == '-') flag = 0; ch = getchar(); }
    while (ch >= '0' && ch <= '9') { X = (X << 1) + (X << 3) + ch - '0'; ch = getchar(); }
    if (flag) return X;
    return ~(X - 1);
}


void dfs(int now){
	if(at == s){
		ans++;
		return;
	}
	for(int i = h[now]; i != -1; i = ne[i]){
		int y = e[i];
		int flag = 1;
		if(y < now) continue;
		for(int j = 0; j < at; ++j){
			int z = a[j];
			if(mp[y][z] == 0){
				flag = 0;
				break;
			}
		}
		if(flag){
			a[at++] = y;
			dfs(y);
			at--;
		}
	}
}

signed main(){
	int t;
	t = read();
	while(t--){
		memset(mp, 0, sizeof mp);
		memset(h, -1, sizeof h);
		idx = 0;
		ans = 0;
		at = 0;
		n = read(), m = read(), s = read();
		for(int i = 0; i < m; ++i){
			int u, v;
			u = read(), v = read();
			add(u, v);
			mp[u][v] = mp[v][u] = 1;
		}
		for(int i = 1; i <= n - s + 1; ++i){
			a[at++] = i;
			dfs(i);
			at--;
		}
		printf("%d\n", ans);
	}
	return 0;
}

E Do not pour out

​ 当d >= 1 时，上面是一个完整的椭圆，用简单的推算直接求出椭圆的 $a$ 和 $b$，椭圆面积公式 $\pi ab$ 即可。

​ 关键是当 d < 1时，椭圆面是会被截去一部分的，如图

​ 设 $mid,t,y$，根据相似三角形得到 $\frac{t}{mid}=\frac{2-y}{2}$，所以有 $t=\frac{(2-y)mid}{2}$
​ 有效底面积是 $\theta +(t-1)sin(\theta )$，其中 $\theta =arccos(1-t)$，这里不需要特判 $t$ 和 1 的大小关系，神奇的三角函数已经帮我们考虑这种关系了。

​ 所以总结一下公式
$$\begin{gathered} ①t=\frac{(2-y)mid}{2} \\ ②\theta =arccos(1-t) \\ ③S(t)=\theta +(t-1)sin(\theta ) \\ ④V={\int }_0^{2}S(t(y))dy \end{gathered}$$
​ 然后就是代换，③代入④得到
$$⑤V={\int }_0^2(\theta +(t-1)sin\theta )dy$$
​ 由①得 $y=2-\frac{2t}{mid}$，所以得到
$$⑥dy=-\frac{2}{mid}sin\theta d\theta$$
​ 由②
$$⑦t=1-cos\theta$$
​ ⑦代入⑤得到
$$⑧V={\int }_0^2(\theta -sin\theta con\theta )dy$$
​ ⑥代入⑧ 得到
$$⑨V=-\frac{2}{mid}{\int }_a^b((\theta -sin\theta cos\theta )sin\theta )d\theta$$
​ 得到的⑨式中的积分区间是发生了变化的，从原来的对 $y$ 积分变成了对 $\theta$ 积分，这里根据 公式⑥进行新区间范围计算 $[arccos(1-mid),arccos(1)]$，所以最后的积分公式为
$$⑩V=-\frac{2}{mid}{\int }_{arccos(1-mid)}^{arccos(1)}(\theta sin\theta -si{n}^2\theta cos\theta )d\theta$$

​ 然后就是求积分了，的到最后需要实现的计算公式
$$V=-\frac{2}{mid}[sin(x)-xcos(x)-\frac{1}{3}si{n}^3(x){]}_{acos(a-mid)}^{acos(1)}$$
​ 这样就可以通过二分求 $mid$，得到 $mid$ 之后，可以求出斜面和底面的夹角 $\alpha$，根据图形得到以下关系
$$cos\alpha =\frac{mid}{\sqrt{2^2+mi{d}^2}}=\frac{S_{底}}{S_0}$$
​ 这里 $S_0$ 表示斜面面积， $S_{底}$ 表示 $S_0$ 在底面的投影，$S_{底}$ 可以用②③公式求出来，所以答案就是
$$S_0=\frac{S_{底}\sqrt{2^2+mi{d}^2}}{mid}$$

#include<bits/stdc++.h>
using namespace std;
const double PI = acos(-1.0);
const double eps = 1e-10;

double f(double x){
	return sin(x) - x * cos(x) - pow(sin(x), 3) / 3;
}

double getV(double mid){
	double a = acos(1 - mid);
	double b = acos(1);
	return -2.0 / mid * (f(b) - f(a));
}

double getS(double t){
	double theta =  acos(1 - t);
	return theta + (t - 1) * sin(theta);
}

int main(){
	int t;
	cin >> t;
	while(t--){
		double d;
		cin >> d;
		if(fabs(d) < eps){
			printf("%.5f\n", 0);
			continue;
		}
		if(d < 1){
			double l = 0, r = 2;
			while(fabs(r - l) > eps){
				double mid = (l + r) / 2;
				if(getV(mid) < PI * d)
					l = mid;
				else r = mid;
			}
			double mid = (l + r) / 2;
			double S = getS(mid);
			double S1 = S * sqrt(4 + mid * mid) / mid;
			printf("%.5f\n", S1);
		}
		else{
			double a = sqrt(1 + (2 - d) * (2 - d));
			double b = 1;
			double S = PI * a * b;
			printf("%.5f\n", S);
		}
	}
	return 0;
}