2016summer 训练第一场
A.http://acm.hdu.edu.cn/showproblem.php?pid=5538
求表面积,只需要将所有的1*1的小块扫描一遍。将每一个块与他相邻四周进行比较,如果该快高度大,则将该快高度减去周围块高度然后累加。复杂度O(nm)
#define _CRT_SECURE_NO_DEPRECATE
#include<iostream>
#include<cmath>
#include<algorithm>
typedef long long LL;
const double EPS = 1e-8;
const double PI = acos(-1.0);
using namespace std;
const int MAXN = 55;
int n, m;
int a[MAXN][MAXN];
int dir[4][2] = { { 0, 1 }, { 0, -1 }, { -1, 0 }, { 1, 0 } };
int main(){
int T, i, j, k;
scanf("%d\n", &T);
while (T--){
memset(a, 0, sizeof(a));
int ans = 0;
scanf("%d%d", &n, &m);
for (i = 1; i <= n; i++)
for (j = 1; j <= m; j++){
scanf("%d", &a[i][j]);
if (a[i][j] > 0)
ans++;
}
for (i = 1; i <= n; i++)
for (j = 1; j <= m; j++){
for (k = 0; k < 4; k++){
int x = i + dir[k][0];
int y = j + dir[k][1];
if (a[x][y] < a[i][j]){
ans += a[i][j] - a[x][y];
}
}
}
printf("%lld\n", ans);
}
return 0;
}
2.http://acm.hdu.edu.cn/showproblem.php?pid=5532
原来序列要么按照升序,要么按照降序去判断,首先可以假定原序列是升序(降序同样来处理).首先将原序列从左到右扫描一遍,不断的比较相邻的两个元素,直到遇到某两元素满足
a[i]>a[i+1]或者扫描到末尾时逃出。若扫描到了末尾,则原序列是增序列,满足条件。若是遇到a[i]>a[i+1]跳出,则我们可以断定,我们要删去的元素一定是a[i]或者a[i+1],可以这样来想,如果删去的是其它元素,则该序列中还是存在a[i]和a[i+1]破坏来递增性。于是问题简单了,只需要分别去掉a[i],和a[i+1]来判断一下原序列是否是递增序列即可。同样的,当假设原序列是递减序列时,处理方法相同。
#define _CRT_SECURE_NO_DEPRECATE
#include<iostream>
#include<cmath>
#include<algorithm>
typedef long long LL;
const double EPS = 1e-8;
const double PI = acos(-1.0);
using namespace std;
const int MAXN =1000001;
int a[MAXN];
int n;
bool UpSroted(){
int i;
for (i = 0; i < n - 1; i++)
if (a[i]>a[i + 1]) //若遇到降序的则跳出
break;
if (i>=n - 2) //a[i+1]后面没有数了
return true;
int j = i + 2;
while (j < n - 1){
if (a[j]>a[j + 1])
break;
j++;
}
if (j<n - 1) //再次出现不符合
return false;
if (a[i] <= a[i + 2]) //去掉a[i+1]即可
return true;
if (a[i + 1]>a[i + 2])
return false;
if (i == 0)
return true;
if (a[i - 1] <= a[i + 1])
return true;
return false;
} bool DownSroted(){
int i;
for (i = 0; i < n - 1; i++)
if (a[i]<a[i + 1]) //若遇到降序的则跳出
break;
if (i >= n - 2) //a[i+1]后面没有数了
return true;
int j = i + 2;
while (j < n - 1){
if (a[j]<a[j + 1])
break;
j++;
}
if (j<n - 1) //再次出现不符合
return false;
if (a[i] >= a[i + 2]) //去掉a[i+1]即可
return true;
if (a[i + 1]<a[i + 2])
return false;
if (i == 0)
return true;
if (a[i - 1] >= a[i + 1])
return true;
return false;
}
int main(){
int T, i;
scanf("%d\n", &T);
while (T--){
scanf("%d", &n);
for (i = 0; i < n; i++)
scanf("%d", &a[i]);
bool flag = UpSroted() || DownSroted();
if (flag)
printf("YES\n");
else
printf("NO\n");
}
return 0;
}
3.http://acm.hdu.edu.cn/showproblem.php?pid=5533
原问题给出的都是整数点,这样以来就只有正方形符合要求了,只需要判断一下给定的是否是4个点,这四个点是否可以构成正方形。
#define _CRT_SECURE_NO_DEPRECATE
#include<iostream>
#include<cmath>
#include<algorithm>
using namespace std;
const double PI = acos(-1.0);
const int N = 300;
const double EPS = 1e-8;//实数精度
//点结构类型
struct Point{
double x, y;
Point(double a = 0, double b = 0){ x = a; y = b; }
};
Point operator-(Point a, Point b){
return Point(a.x - b.x, a.y - b.y);
}
//重载==,判断点a,b是否相等
bool operator==(Point a, Point b){
return abs(a.x - b.x) < EPS&&abs(a.y - b.y) < EPS;
}
//比较实数r1与r2的大小关系
int RlCmp(double r1, double r2 = 0){
if (abs(r1 - r2) < EPS)
return 0;
return r1>r2 ? 1 : -1;
}
double Dot(Point p1, Point p2, Point p3, Point p4){
Point a = p2 - p1;
Point b = p4 - p3;
return a.x*b.x + a.y*b.y;
}
//检查以p1p2和p3p4为对角线是否可以构成正方形
bool firstCheck(Point p1, Point p2, Point p3, Point p4){
Point mid0 = Point((p1.x + p2.x) / 2, (p1.y + p2.y) / 2);
Point mid1 = Point((p3.x + p4.x) / 2, (p3.y + p4.y) / 2);
if (!(mid0 == mid1)) //如果中点不重合
return false;
return RlCmp(Dot(p1, p2, p3, p4)) == 0; //对角线垂直
}
bool isSqual(Point P[]){
return firstCheck(P[0], P[1], P[2], P[3]) ||
firstCheck(P[0], P[2], P[1], P[3]) ||
firstCheck(P[0], P[3], P[1], P[2]);
}
int main(){
Point P[N];
int T,i,n;
double x, y;
scanf("%d", &T);
while (T--){
scanf("%d", &n);
for (i = 0; i < n; i++){
scanf("%lf%lf", &x, &y);
P[i] = Point(x, y);
}
if (n != 4){
printf("NO\n");
continue;
}
if (isSqual(P))
printf("YES\n");
else
printf("NO\n");
}
return 0;
}
4.http://acm.hdu.edu.cn/showproblem.php?pid=5536
该问题是贪心+字典树。虽然网上有O(n^3)爆力剪枝可以过,表示不会。将每一个数字的二进制从高位到低位存入到字典树中形成01串。这时候我们需要一个节点变量v,每次插入只需要将v++.接下来,枚举a[i],a[j](j>i),首先在字典树中删掉a[i]和a[j],删除做法很简单,将对应节点位置v--即可,然后用a[i]+a[j]在字典树中进行01串的匹配,匹配完后再插入a[i],a[j]准备下一次的匹配。匹配方法采用贪心思想,依次从高位向低位匹配,若当前数该位为1则去优先匹配0,若没有0匹配,则只能匹配1了。当前位置是0,则优先去匹配1,当没有1匹配,就只能匹配0了。如此下去,直到匹配到最后一位。这样有一个问题,每一个数二进制串长度不一样,给匹配带来不便,做法是,将所有的串高位补0,使得长度为32位。这样就可开心的匹配了。
#define _CRT_SECURE_NO_DEPRECATE
#include<iostream>
#include<algorithm>
typedef long long LL;
const int MAXN = 1000+10;
using namespace std;
struct TrieNode{
TrieNode(){ memset(next, 0, sizeof(next)); v = 0; }
TrieNode* next[2];
LL v;
};
LL MAX;
TrieNode *root;
void Insert(LL x){
TrieNode*p = root;
MAX = 1;
MAX <<= 32;
for (LL i =MAX; i >0; i>>=1){
LL id =(i&x)>0;
if (p->next[id] == NULL)
p->next[id] = new TrieNode;
p = p->next[id];
p->v++;
}
}
void Delete(LL x){
TrieNode*p = root;
MAX = 1;
MAX <<= 32;
for (LL i = MAX; i >0; i >>= 1){
LL id = (i&x)>0;
p = p->next[id];
p->v--;
}
}
LL getMAX(LL x){
TrieNode *p = root;
MAX = 1;
MAX <<= 32;
LL rt,ans=x,i=-1;
for (LL i = MAX; i >0;i>>=1){
LL id = (x&i)>0;
if (p->next[id ^ 1] && p->next[id ^ 1]->v > 0){
if ((id == 0)) //说明x当前位为0,即将变为1
ans +=i ;
p = p->next[id ^ 1];
continue;
}
if (id) //如果x当前位置为1,则即将要变为0
ans -= i;
p = p->next[id];
}
return ans;
}
void Free(TrieNode*T){
if (T){
for (int i = 0; i < 2; i++)
Free(T->next[i]);
free(T);
T = NULL;
}
}
LL a[MAXN];
int main(){
LL n,i,j,T;
scanf("%I64d", &T);
while (T--){
root = new TrieNode;
scanf("%I64d", &n);
for (i = 0; i < n; i++){
scanf("%I64d", &a[i]);
Insert(a[i]);
}
LL ans =0;
for (i = 0; i < n; i++){
Delete(a[i]);
for (j = i + 1; j < n; j++){
Delete(a[j]);
ans = max(ans, getMAX(a[i] + a[j]));
Insert(a[j]);
}
Insert(a[i]);
}
printf("%I64d\n",ans);
Free(root);
}
return 0;
}
5.http://acm.hdu.edu.cn/showproblem.php?pid=5534
思维比较巧妙的一道DP题目。首先n个节点的无向树有2*(n-1)个度。每一节点至少有一个度,所以可以这样考虑,首先给每个节点分配一个度,这样还剩下2*n-2个度,于是问题就转化为了完全背包问题。需要注意的是此时背包总容量是2*n-2,价值转移方程为p[i]-p[0];
#define _CRT_SECURE_NO_DEPRECATE
#include<iostream>
#include<algorithm>
using namespace std;
const int MAXN = 2005 * 2 + 1;
#define INF 0x3f3f3f3f
int n;
int p[MAXN],f[MAXN],V;
//完全背包问题
void ComplePack(int w, int i){
for (int j = i; j <= V; j++)
f[j] = max(f[j],f[j-i] + w);
}
int main(){
int i, T;
scanf("%d", &T);
while (T--){
scanf("%d", &n);
for (i = 0; i < n - 1; i++)
scanf("%d", &p[i]);
V = n - 2;
for (i = 1; i < n - 1; i++)
p[i] -= p[0]; //价值转移
for (i = 1; i < n - 1; i++)
f[i] = -INF;
f[0] = p[0] * n;
for (i = 1; i <= n - 2; i++){
ComplePack(p[i], i);
}
printf("%d\n",f[n-2]);
}
return 0;
}
6.http://acm.hdu.edu.cn/showproblem.php?pid=5531
题目很长,其实题目意思很简单。题意:在一个多边形每一个顶点以顶点为圆心都构造一个园,相邻顶点的园外切,不相邻的顶点的圆无限制。问是否存在这样的构造,存在则还要要求求出所有圆面积和的最小值,以及对应的园的半径。
设多边形的n个点是p1,p2,...pn.n个园半径依次为r1,r2,r3...rn.设d1=p1p2,d2=p2p3,...,dn-1=pn-1pn,dn=pnp1.于是由园相切,我们可以得到n个关系式:
r1+r2=d1
r2+r3=d2
.....
rn-1+rn=dn-1
rn+r1=dn
这时候我们要对n分奇偶数讨论。设数组dis[i]=d[1]-d[2]+d[3]-d[4]...+(-1)i-1*d[i];如果n是奇数,则迭代可知上面n个等式必然有唯一解,求出此解,然后看是否所有的解大于等0即可。如果n是偶数,则上面等式前n项可以变为
r1+r2=dis[1]
r1-r3=dis[2]
r1+r4=dis[3]
....
r1+rn=dis[n-1]
r1+rn=d[n]
显然只有当d[n]==dis[n-1]时才有解。于是接下来只有n-1个方程和n个未知量,一定有解。我们注意到0<=ri<=min(d[i-1],d[i]),(i=2,3..n).而0<=r1<=min(d1,dn)
.利用以上不等式可以求出r1的取值区间。而目标函数area=(r12+r22...+rn2) *π。aera是关于r1的二次函数,对应开口向上,对称轴为(dis[1]+disp[2]...+dis[n-1])/n.
这样最值问题就求解了。
#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);
const int MAXN = 100000 + 10;
struct Point{
double x, y;
};
struct Intervel{
double s, e;
Intervel(double a = 0, double b = 0) :s(a), e(b){}
};
int n;
Point P[MAXN];
double r[MAXN], d[MAXN],dis[MAXN];
double Distance(Point p1, Point p2){
return sqrt((p1.x - p2.x)*(p1.x - p2.x) + (p1.y - p2.y)*(p1.y - p2.y));
}
int RlCmp(double r1, double r2 = 0){
if (abs(r1 - r2) < EPS)
return 0;
return r1>r2?1 : -1;
}
bool Intersect(Intervel x, Intervel y,Intervel &ans){
if (RlCmp(x.e, y.s) == -1)
return false;
if (RlCmp(y.e, x.s) == -1)
return false;
ans =Intervel(max(x.s, y.s), min(x.e, y.e));
return true;
}
bool isPossible(double &area){
if (n & 1){
r[1] = dis[n]/ 2;
if (RlCmp(r[1])==-1)
return false;
area = r[1] * r[1];
for (int i = 2; i <= n; i++){
if (i & 1)
r[i] = r[1] - dis[i-1];
else
r[i] = dis[i-1] - r[1];
if (RlCmp(r[i])==-1)
return false;
area += r[i] * r[i];
}
area *= pi;
return true;
}
else{
if (RlCmp(dis[n]))
return false;
bool flag = 1;
Intervel curr(0,min(d[1],d[n]));
for (int i = 2; i <=n; i++){
if (i & 1)
flag = Intersect(curr, Intervel(dis[i - 1], dis[i - 1] + min(d[i], d[i - 1])), curr);
else
flag = Intersect(curr, Intervel(dis[i - 1] - min(d[i], d[i - 1]),dis[i-1]), curr);
if (!flag)
return false;
}
double t = 0;
for (int i = 1; i <= n - 1; i++)
t += dis[i];
t /= n;
if (t < curr.s)
r[1] = curr.s;
else if (t>curr.e)
r[1] = curr.e;
else
r[1] = t;
area = r[1] * r[1];
for (int i = 2; i <= n; i++){
if (i & 1)
r[i] = r[1] - dis[i - 1];
else
r[i] = dis[i - 1] - r[1];
area += r[i] * r[i];
}
area *= pi;
return true;
}
}
int main(){
int T, i;
scanf("%d", &T);
while (T--){
scanf("%d", &n);
for (i = 1; i <= n; i++)
scanf("%lf%lf", &P[i].x, &P[i].y);
for (i = 1; i < n; i++)
d[i] = Distance(P[i], P[i + 1]);
d[n] = Distance(P[1], P[n]);
dis[0] = 0;
for (i = 1; i <= n; i++){
if (i & 1)
dis[i]=dis[i-1]+d[i];
else
dis[i]=dis[i-1]-d[i];
}
double area;
bool ans = isPossible(area);
if (ans){
printf("%.2lf\n", area);
for (i = 1; i <= n; i++)
printf("%.2lf\n",r[i]);
}
else
printf("IMPOSSIBLE\n");
}
return 0;
}
7.http://acm.hdu.edu.cn/showproblem.php?pid=5527
一道很好的贪心题目。我们要考虑选取数来凑p,但是呢,每个硬币个数是有限制的,顺着贪心是无法下手的。总的硬币个数num和总的钱total是知道的,求凑p得到最大也就是求凑total-p的最小,因为他们的和就是num。于是我们来逆向考虑求最小值。这样我们从最大的硬币选取,一直贪心选下来,基于这样算法成立的条件是,前面硬币面值一定是后面硬币面值的约数,只有这样,前面硬币能达到超过该硬币的值,一定可以用该硬币来替代而使得总硬币个数更小。硬币面值 1 5 10 20 50 100 200 500 1000 2000,我们发现有两对不满足,20 50和200 500,于是会产生这样一个问题。比如我现在有20硬币*3,50硬币*1,当前要选取等于60的硬币,如果贪心来选首选的是50,而实际应该选择是20*3.产生这个原因是由于20不是50的约数。解决的办法就是对50,500进行特判,枚举他们的奇偶情况,因为2个50可以当做100来使用。1. 50奇500偶数。2. 50偶数 500奇数
3. 50奇数 500奇数 4. 50偶数 500偶数。然后把4种情况跑一下就OK了。
#define _CRT_SECURE_NO_DEPRECATE
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
LL p, total,num;
LL nl[11];
LL c[11] = { 0, 1, 5, 10, 20, 50, 100, 200, 500, 1000, 2000 };
LL Compute(LL pp){
LL i, ans = 0;
if (pp<0)
return -1;
for (i = 10; i > 0; i--){
LL k = min(pp/c[i],nl[i]);
if ((k&1)&&(i == 5 || i == 8))
k--; //只可以用掉偶数个
ans += k;
pp-= k*c[i];
}
if (pp!= 0)
return -1;
return num - ans;
}
LL Slove(){
LL ans = -1,pp,ret;
//下面对四中情况进行枚举
if (nl[5] > 0){
nl[5]--;
pp = p - c[5];
ret=Compute(pp);
ans = max(ans, ret-1);
nl[5]++;
}
if (nl[8] > 0){
nl[8]--;
pp = p - c[8];
ret = Compute(pp);
ans = max(ans, ret-1);
nl[8]++;
}
if (nl[5] > 0 && nl[8] > 0){
nl[5]--;nl[8]--;
pp = p - c[5] - c[8];
ret = Compute(pp);
ans = max(ans, ret-2);
nl[5]++; nl[8]++;
}
pp = p;
ret=Compute(pp);
ans = max(ans, ret);
return ans;
}
int main(){
LL i, T;
scanf("%I64d", &T);
while (T--){
total = num = 0;
scanf("%I64d", &p);
for (i = 1; i <= 10; i++)
scanf("%I64d", &nl[i]);
for (i = 1; i <=10; i++){
num += nl[i]; //计算总个数
total += c[i] * nl[i]; //计算总money
}
p = total - p; //从反面考虑 LL ans = Slove();
printf("%I64d\n", ans);
}
return 0;
}
2016summer 训练第一场的更多相关文章
- 大家一起做训练 第一场 E Number With The Given Amount Of Divisors
题目来源:CodeForce #27 E 题目意思和题目标题一样,给一个n,求约数的个数恰好为n个的最小的数.保证答案在1018内. Orz,这题训练的时候没写出来. 这道题目分析一下,1018的不大 ...
- 大家一起做训练 第一场 B Tournament
题目来源:CodeForce #27 B 有n个人比赛,两两之间都有一场比赛,一共 n * (n - 1) / 2 场比赛.每场比赛的记录方式是 a b,表示在a和b的比赛中,a胜出,b失败. 经过研 ...
- 牛客网多校训练第一场 I - Substring(后缀数组 + 重复处理)
链接: https://www.nowcoder.com/acm/contest/139/I 题意: 给出一个n(1≤n≤5e4)个字符的字符串s(si ∈ {a,b,c}),求最多可以从n*(n+1 ...
- 2015多校联合训练第一场Tricks Device(hdu5294)
题意:给一个无向图,给起点s,终点t,求最少拆掉几条边使得s到不了t,最多拆几条边使得s能到t 思路: 先跑一边最短路,记录最短路中最短的边数.总边数-最短边数就是第二个答案 第一个答案就是在最短路里 ...
- Gym-101653:acific Northwest Regional Contest (2019训练第一场)
本套题没有什么数据结构题,图论题,唯一有价值的就是Q题博弈,在最后面,读者可以直接拉到最下面. (还剩下两个,估计每什么价值的题,懒得补了 M .Polyhedra pro:欧拉公式,V-E+F=2: ...
- 大家一起做训练 第一场 G CD
题目来源:UVA 624 题目的意思就是:我现在需要从 t 张CD中拿出一部分来,尽可能的凑出接近 N 这么久的音乐,但是不能超过 N. CD不超过20张,每张长度不超过 N ,不能重复选. 一个很简 ...
- 大家一起做训练 第一场 A Next Test
题目来源:CodeForce #27 A 题目的意思简而言之就是要你输出一个没有出现过的最小的正整数. 题意如此简单明了,做法也很明了. 直接读入所有的数,然后排个序,设置个变量从1开始,出现过+1, ...
- 牛客网多校训练第一场 J - Different Integers(树状数组 + 问题转换)
链接: https://www.nowcoder.com/acm/contest/139/J 题意: 给出n个整数的序列a(1≤ai≤n)和q个询问(1≤n,q≤1e5),每个询问包含两个整数L和R( ...
- 牛客网多校训练第一场 F - Sum of Maximum(容斥原理 + 拉格朗日插值法)
链接: https://www.nowcoder.com/acm/contest/139/F 题意: 分析: 转载自:http://tokitsukaze.live/2018/07/19/2018ni ...
随机推荐
- 图说不为人知的IT传奇故事-2-IBM咬不动的胡桃
此系列文章为“图说不为人知的IT传奇故事”,各位大忙人可以在一分钟甚至几秒内了解把握整个内容,真可谓“大忙人的福利”呀!!希望各位IT界的朋友在钻研技术的同时,也能在文学.历史上有所把握.了解这些故事 ...
- Python之threading多线程
1.threading模块是Python里面常用的线程模块,多线程处理任务对于提升效率非常重要,先说一下线程和进程的各种区别,如图 概括起来就是 IO密集型(不用CPU) 多线程计算密集型(用CPU) ...
- SGX技术初探
一.SGX技术背景 1.1 SGX技术定义 SGX全称Intel Software Guard Extensions,顾名思义,其是对因特尔体系(IA)的一个扩展,用于增强软件的安全性.这种方式并不是 ...
- 团队Alpha版本冲刺(一)
目录 组员情况 组员1(组长):胡绪佩 组员2:胡青元 组员3:庄卉 组员4:家灿 组员5:凯琳 组员6:丹丹 组员7:家伟 组员8:政演 组员9:黄鸿杰 组员10:刘一好 组员11:何宇恒 展示组内 ...
- POJ-2159 最小费用最大流
Going Home 自己写的第一道费用流,图建好一波板子AC.不过还是有几个地方有点迷. 先来 ...
- 历史Linux镜像的问题修复方案
历史Linux镜像创建的ECS云服务器,可能存在NTP没有配置,YUM没有配置,还可能存在最近暴漏较高的安全漏洞,请按照以下步骤进行修复,可以让您的云服务器更加安全,还可以使用阿里云提供的YUM服务进 ...
- vue - computed
computed 的作用主要是对原数据进行改造输出.改造输出:包括格式的编辑,大小写转换,顺序重排,添加符号……. 一.格式化输出结果: 我们先来做个读出价格的例子:我们读书的原始数据是price:1 ...
- FZU Super A^B mod C(欧拉函数降幂)
Problem 1759 Super A^B mod C Accept: 878 Submit: 2870 Time Limit: 1000 mSec Memory Limit : 327 ...
- 【BZOJ1717】[Usaco2006 Dec]Milk Patterns 产奶的模式 (二分+SA)
求重复k次的最长重复子串,解法见罗穗骞大神的后缀数组论文 ; var x,y,rank,sa,h,s,num,c:..maxn] of longint; n,time:longint; functio ...
- [SDOI2008][luogu2463] Sandy的卡片 [kmp]
题面 传送门 思路 这道题里面有三个主要问题: 1.由"数值相等"变成了"加上一个整数以后数值相等"(减去等价于加负数) 2.由"最多匹配多少位(从第 ...