2016summer 训练第二场

1.http://acm.hdu.edu.cn/showproblem.php?pid=5112

#define _CRT_SECURE_NO_DEPRECATE
#include<iostream>
#include<algorithm>
using namespace std;
const int MAXN = 100000+10;
struct Node{
	int t;
	int x;
	bool operator<(Node&a){
		return t < a.t;
	}
};
Node S[MAXN];
int main(){
	int i, T,n;
	int iCase = 1;
	scanf("%d", &T);
	while(T--){
		scanf("%d", &n);
		for (i = 0; i < n; i++)
			scanf("%d%d", &S[i].t, &S[i].x);
		sort(S, S + n);
		double Max = 0;
		for (i = 0; i < n - 1; i++){
			double v = abs(S[i + 1].x - S[i].x)*1.0 / (S[i + 1].t - S[i].t);
			Max = max(Max, v);
		}
		printf("Case #%d: %.2lf\n",iCase++, Max);
	}
	return 0;
}

　　2.http://acm.hdu.edu.cn/showproblem.php?pid=5122

一开始一看到以为求逆序来搞，一看题就写了个数状数组，nlogn超时了。当然一定是太性急了，这个题其实就是判断下它后面是否有数字比它小，有的话需要将其往后移一次，否则不变。这样只需要判断下每个数后面是否有有比它小的数，于是将数组倒置，即为判断每个数前面是否有比它大的数，只需要从前向后扫描一遍，需要一个数Min记录当前已扫描结束的序列中最小的数。这样只要碰到后迷死你右数大于这个最小的数就计数一次，否则更新这个Min。

#define _CRT_SECURE_NO_DEPRECATE
#include<iostream>
#include<algorithm>
using namespace std;
const int MAXN = 10000001;
int a[MAXN];
int main(){
	int i, T, n;
	int iCase = 1;
	scanf("%d", &T);
	while (T--){
		scanf("%d", &n);
		for (i = n; i > 0; i--)
			scanf("%d", &a[i]);
		int Min = a[1];
		int ans = 0;
		for (i = 2; i <= n; i++){
			if (a[i] > Min)
				ans++;
			else
				Min = a[i];
		}
		printf("Case #%d: %d\n",iCase++,ans);
	}
	return 0;
}

　　3.http://acm.hdu.edu.cn/showproblem.php?pid=5120

模板很重要啊。这里主要的是求两个圆相交的面积，有了模板直接搞。最终面积=大圆交大圆-2*大圆交小圆+小圆交小圆。

#define _CRT_SECURE_NO_DEPRECATE
#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
const double EPS = 1e-8;
const double pi = acos(-1.0);
struct Point{
	double x, y;
	Point(double a=0, double b=0) :x(a), y(b){}
};
struct Circle{
	Circle(Point a, double x) :o(a), r(x){}
	Point o;
	double r;
};
int RlCmp(double r1, double r2){
	if (abs(r1 - r2) < EPS)
		return 0;
	return r1 < r2 ? -1 : 1;
}
Point operator-(Point a, Point b){
	return Point(a.x - b.x, a.y - b.y);
}
Point operator+(Point a, Point b){
	return Point(a.x + b.x, a.y + b.y);
}
double Mold(Point a){
	return sqrt(a.x*a.x + a.y*a.y);
}
double Dis(Point a, Point b){
	return Mold(a - b);
}
//园与圆相交面积模板
double CircleCrossArea(Circle A,Circle B){
	double r1 = A.r, r2 = B.r;
	double d = Dis(A.o, B.o),r=min(r1,r2);
	if (RlCmp(d, r1 + r2) >= 0)
		return 0;        //相离或者外切
	if (RlCmp(d, abs(r1 - r2)) <= 0)
		return pi*r*r;       //内含
	//将r1放在圆心
	double x1 = (d*d + r1*r1 - r2*r2) / (2 * d);
	double s1 = x1*sqrt(r1*r1 - x1*x1) - r1*r1*acos(x1/r1);
	//将r2放在圆心
	double x2 = (d*d + r2*r2 - r1*r1) / (2 * d);
	double s2 = x2*sqrt(r2*r2 - x2*x2) - r2*r2*acos(x2 / r2);
	return abs(s1 + s2);
}
int main(){
	int T, iCase = 1;
	Point q,o;
	double r1, r2;
	scanf("%d", &T);
	while (T--){
		scanf("%lf%lf",&r1,&r2);
		scanf("%lf%lf",&q.x,&q.y);
		scanf("%lf%lf", &o.x,&o.y);
		Circle A1(q, r1);
		Circle A2(q, r2);
		Circle B1(o, r1);
		Circle B2(o, r2);
		double a = CircleCrossArea(A2, B2);
		double b = CircleCrossArea(A2, B1);
		double c = CircleCrossArea(A1, B1);
		printf("Case #%d: %.6lf\n", iCase++, a - 2 * b + c);
	}
	return 0;
}

　　4.http://acm.hdu.edu.cn/showproblem.php?pid=5113

一个很好的搜索题。对矩阵从左到右，一行一行的进行搜索(不要把回溯写错或者写漏掉.....老毛病).注意剪枝，一个剪枝条件是,当某一种颜色的数量超过(格子总数+1)/2时，剩下颜色拼出的一定是不合法的。

#define _CRT_SECURE_NO_DEPRECATE
#include<iostream>
#include<algorithm>
using namespace std;
const int MAXN = 7;
int a[7][7];
int b[MAXN*MAXN];
int k, n, m,iCase;
int dir[4][2] = { { 1, 0 }, { -1, 0 }, { 0, 1 }, { 0, -1 } };
bool ans;
bool isLeg(int x, int y,int val){
	//判断在a[x][y]处填颜色val是否合法
	return x - 1 >=0 && a[x - 1][y] != val
		&&y - 1 >=0 && a[x][y - 1] != val;
}
void DFS(int r,int c,int cur){
	for (int i = 1; i <= k; i++){
		if (b[i]>(cur+1)/2)   //剪枝
			return;
	}
	if (c > m){
		c = 1;
		r++;
	}
	if (r > n){   
		//r>n时，表明已经搜索完成
		ans = true;
		printf("Case #%d:\nYES\n", iCase++);
		for (int i = 1; i <= n; i++)
		for (int j = 1; j <= m; j++){
			printf(j == m ? "%d\n" : "%d ", a[i][j]);
		}
	}
	if (ans)
		return;
	for (int i = 1; i <= k; i++){
		if (b[i] > 0 && isLeg(r, c, i)){
			a[r][c] = i;
			b[i]--;
			DFS(r, c+1,cur-1);
			//记得回溯洛......
			b[i]++;
			a[r][c] = 0;
		}
	}
}
int main(){
	int T, i,j;
	scanf("%d", &T);
	iCase = 1;
	while (T--){
		scanf("%d%d%d", &n, &m,&k);
		for (i = 1; i <= k; i++)
			scanf("%d", &b[i]);
		memset(a, 0, sizeof(a));
		ans = false;
		DFS(1, 1,n*m);
		if (!ans)
			printf("Case #%d:\nNO\n", iCase++);
	}
	return 0;
}

　　5.http://acm.hdu.edu.cn/showproblem.php?pid=5119

背包问题，真难看出来....。用dp[i][j]表示前i个数中异或值为j的选法数。则得到状态转移方程dp[i][j]=dp[i-1][j]+dp[i-1][j^a[i]];   这里有一个性质:若a^b=c，则b^c=a,a^c=b;题目还有一个坑点，Mi<=10^6，但是他们的异或值可能会超过10^6,于是数组开在10^6大一点不行，但是开在10^7又会爆掉。考虑到10^6的二进制数异或值不会大于2^20次方，从而可以拿这个数来开数组。另外，第i行只与第i-1行有关，可以使用滚动数组。

#define _CRT_SECURE_NO_DEPRECATE
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
const int MAXN = 1<<20;
LL n, m;
LL a[45];
LL f[2][MAXN];
int main(){
	LL i, iCase, T,j;
	iCase = 1;
	scanf("%I64d", &T);
	while (T--){
		scanf("%I64d%I64d", &n, &m);
		for (i = 1; i <= n; i++)
			scanf("%I64d", &a[i]);
		memset(f, 0, sizeof(f));
		//从前0个数中选取0个数，异或值等于0,方法有一种，就是什么也不选
		f[0][0] = 1;
		for (i = 1; i <= n; i++){
			for (j=0;j<MAXN; j++){
					f[i%2][j] = f[(i-1)%2][j] + f[(i-1)%2][j^a[i]];
			}
		}
		LL ans = 1;
		ans <<= n;
		for (i =0; i <m;i++)
			ans-= f[n%2][i];
		printf("Case #%I64d: %I64d\n",iCase++,ans);
	}
	return 0;
}