BZOJ刷题记录---普及组难度
总目录详见https://blog.csdn.net/mrcrack/article/details/90228694
| 序号 | 题号 | 算法 | 思想难度 | 实现难度 | 总难度 | 推荐指数 |
| 1 | 2659 | 数学 | 20 | 7 | 27 | 2 |
| 2 | 1968 | 枚举 | 15 | 13 | 28 | 2 |
| 3 | 2761 | 去重 | 12 | 18 | 30 | 6 |
| 4 | 3767 | 高精度A+B | 10 | 20 | 30 | 2 |
| 5 | 1432 | 找规律 | 27 | 5 | 32 | 3 |
| 6 | 3098 | 随机卡HASH | 23 | 10 | 33 | 2 |
| 7 | 1800 | 统计 | 20 | 16 | 36 | 3 |
| 8 | 1218 | 枚举 | 20 | 17 | 37 | 2 |
| 9 | 4001 | 找规律 | 28 | 9 | 37 | 3 |
| 10 | 2760 | 字符串处理 | 12 | 26 | 38 | 1 |
| 11 | 1088 | DP/判断 | 22 | 17 | 39 | 6 |
| 12 | 1034 | 贪心 | 24 | 16 | 40 | 6 |
| 13 | 1207 | DP | 23 | 17 | 40 | 2 |
| 14 | 2748 | DP | 21 | 19 | 40 | 4 |
| 15 | 4302 | 分类讨论 | 29 | 11 | 40 | |
| 16 | 1083 | 最小生成树 | 17 | 25 | 42 | 9 |
| 17 | 2563 | 转换后排序 | 27 | 15 | 42 | 7 |
| 18 | 3097 | 构造卡HASH | 27 | 15 | 42 | 2 |
| 19 | 1022 | 简单博弈 | 30 | 13 | 43 | 3 |
| 20 | 2465 | 贪心 | 25 | 20 | 45 | 4 |
| 21 | 2122 | 暴力 | 28 | 18 | 46 | 2 |
| 22 | 1197 | DP | 31 | 19 | 50 | 2 |
| 23 | 1295 | 搜索 | 20 | 30 | 50 | 5 |
| 24 | 1821 | 最小生成树 | 25 | 25 | 50 | 4 |
| 25 | 2134 | 期望DP | 30 | 20 | 50 | 4 |
测评地址：
1.bzoj 2659 //2659: [Beijing wc2012]算不出的算式 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2659
2.bzoj 1968 //1968: [Ahoi2005]COMMON 约数研究 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1968
3.bzoj 2761 //2761: [JLOI2011]不重复数字 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2761
4.bzoj 3767 此题存在版权，故不再支持提交，保留在此只供大家参考题面！ 望见谅！
5.bzoj 1432 //1432: [ZJOI2009]Function //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1432
6.bzoj 3098 //3098: Hash Killer II //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=3098
7.bzoj 1800 //1800: [Ahoi2009]fly 飞行棋 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1800
8.bzoj 1218 //1218: [HNOI2003]激光炸弹 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1218
9.bzoj 4001 //4001: [TJOI2015]概率论 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=4001
10.bzoj 2760 //2760: [JLOI2011]小A的烦恼 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2760
11.bzoj 1088 //1088: [SCOI2005]扫雷Mine //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1088
12.bzoj 1034 //1034: [ZJOI2008]泡泡堂BNB //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1034
13.bzoj 1207 //1207: [HNOI2004]打鼹鼠 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1207
14.bzoj 2748 //2748: [HAOI2012]音量调节 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2748
15.bzoj 4302 //4302: Hdu 5301 Buildings //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=4302
16.bzoj 1083 //1083: [SCOI2005]繁忙的都市 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1083
17.bzoj 2563 //2563: 阿狸和桃子的游戏 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2563
18.bzoj 3097 //3097: Hash Killer I //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=3097
19.bzoj 1022 //1022: [SHOI2008]小约翰的游戏John //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1022
20.bzoj 2465 //2465: [中山市选2009]小球 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2465
21.bzoj 2122 //2122: 工作评估 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2122
22.bzoj 1197 //1197: [HNOI2006]花仙子的魔法 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1197
23.bzoj 1295 //1295: [SCOI2009]最长距离 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1295
24.bzoj 1821 //1821: [JSOI2010]Group 部落划分 Group //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1821
25.bzoj 2134 //2134: 单选错位 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2134
题解：
1.bzoj 2659 //2659: [Beijing wc2012]算不出的算式 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2659
//2659: [Beijing wc2012]算不出的算式
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2659
//可能遇到int溢出，采用unsigned long long
//样例通过，提交Time_Limit_Exceed，顿时傻眼。
//看来要用等差数列求和公式，被公式中的字母迷惑了。
//马上修改,发现，遇到了取整的问题。难办啊2019-5-19 22:21
//https://www.cnblogs.com/zhber/p/4116413.html此文思路写得不错
//先画出个坐标系，标出坐标为(p,q)的点，然后从(0,0)到(p,q)连一条线，画出以这条线为对角线的长方形。显然直线的斜率是p/q。
//你会发现Σ[kp/q]其实就是直线下方的整点数，而Σ[kq/p]就是直线上方的整点数。
//计算p==q的时候还要考虑直线上有整点的情况
//说实话，打表，探得公式，其实很困难。
//样例解释如下：
//样例通过，提交Wrong_Answer，阅读代码，发现
//else printf("%llu\n",(1+(p-1)/2)/2*(p-1)/2+(1+(q-1)/2)/2*(q-1)/2);//此处写成else printf("%llu\n",(1+(p-1)/2)*(p-1)/2+(1+(q-1)/2)*(q-1)/2);//等差数列，求和
//提交AC。2019-5-20 22:09
//bzoj是ACM提交模式。
#include <stdio.h>
#define ULL unsigned long long
int main(){
ULL p,q,k,L,sum=0;
scanf("%llu%llu",&p,&q);
if(p!=q)printf("%llu\n",(p-1)*(q-1)/4);
else printf("%llu\n",(1+(p-1)/2)/2*(p-1)/2+(1+(q-1)/2)/2*(q-1)/2);//此处写成else printf("%llu\n",(1+(p-1)/2)*(p-1)/2+(1+(q-1)/2)*(q-1)/2);//等差数列，求和
return 0;
}
2.bzoj 1968 //1968: [Ahoi2005]COMMON 约数研究 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1968
//1968: [Ahoi2005]COMMON 约数研究
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1968
//0 < N < 1000000要么O(logn)，要么O(nlogn)
//想到了 素数筛，想到了 约数成对出现 还想到了质因数分解6=2^1*3^1 约数个数(1+1)*(1+1)=4
//先写个素数筛 还不错，能独立写出。2019-5-21
//质数的约数为2，和数的约数分解质因数来处理
//2^20=1024*1024，故具体到某个数，阶数不超过20，故在筛的过程中，再计算 阶数 ,不超时，极限是20*10^6=2*10^7
//n=1无需特判，可读程序，即明白。
//测试了n=1000000,输出13970034，发现，还是要超时
//将以下两句进行优化
//for(j=1;j<=cnt;j++)
// if(x%prime[j]==0)break;
//改成 while(x%prime[j])j++;
//测试了，n=1000000,输出13970034 发现还是超时
//n=100000,秒出，如何才能不超时？
//提交，果不其然Time_Limit_Exceed。2019-5-21
//以下为超时算法。2019-5-21
//https://blog.csdn.net/aarongzk/article/details/50647495此文思路奇特，没有经过训练，确实难以想到，摘抄如下
//直接求每个数的因数个数是比较麻烦的，我们可以换一个角度考虑。
//对于每个数i，在1到n的范围内有多少个数是它的倍数？答案显然是n/i。于是最终的答案就是∑(n/i)。
//测试了n=1000000，真是秒出啊。
//样例通过，提交AC。2019-5-21
//以下为AC代码。
#include <stdio.h>
#define LL long long
LL ans=0;
int main(){
int n,i;
scanf("%d",&n);
for(i=1;i<=n;i++)ans+=n/i;
printf("%lld\n",ans);
return 0;
}
//1968: [Ahoi2005]COMMON 约数研究
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1968
//0 < N < 1000000要么O(logn)，要么O(nlogn)
//想到了 素数筛，想到了 约数成对出现 还想到了质因数分解6=2^1*3^1 约数个数(1+1)*(1+1)=4
//先写个素数筛 还不错，能独立写出。2019-5-21
//质数的约数为2，和数的约数分解质因数来处理
//2^20=1024*1024，故具体到某个数，阶数不超过20，故在筛的过程中，再计算 阶数 ,不超时，极限是20*10^6=2*10^7
//n=1无需特判，可读程序，即明白。
//测试了n=1000000,输出13970034，发现，还是要超时
//将以下两句进行优化
//for(j=1;j<=cnt;j++)
// if(x%prime[j]==0)break;
//改成 while(x%prime[j])j++;
//测试了，n=1000000,输出13970034 发现还是超时
//n=100000,秒出，如何才能不超时？
//提交，果不其然Time_Limit_Exceed。2019-5-21
//以下为超时算法。2019-5-21
#include <stdio.h>
#include <string.h>
#define LL long long
int prime[10000],not_prime[1000100],cnt=0;
LL sum=1;
int f(int x){//计算x的约数的个数
int ans=1,c,j=1;
while(x>1){
c=1;
while(x%prime[j])j++;
while(x%prime[j]==0)c++,x/=prime[j];
ans*=c;
}
return ans;
}
void linear_shaker(int n){
int i,j;
memset(not_prime,0,sizeof(not_prime));
for(i=2;i<=n;i++){
if(!not_prime[i])prime[++cnt]=i;
for(j=1;prime[j]*i<=n;j++){//忘了加次层次的｛｝
not_prime[prime[j]*i]=1;
if(i%prime[j]==0)break;
}
if(!not_prime[i])sum+=2;//此处写成 if(not_prime[i])sum+=2;
else sum+=f(i);//一开始将此if else 写到for(i=2循环之外了
}
}
int main(){
int i,n;
scanf("%d",&n);
linear_shaker(n);
printf("%lld\n",sum);
return 0;
}
3.bzoj 2761 //2761: [JLOI2011]不重复数字 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2761
//2761: [JLOI2011]不重复数字
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2761
//思路：结构体记录位置，值
//按值快排，
//结构题复制，去重
//按位置快排，输出结果
//样例通过，提交Wrong_Answer，真叫没想法。
//仔细想了想，快排是不稳定排序，问题强调"只保留第一次出现的数"
//马上对cmp1进行修改
//int cmp1(node a,node b){
// return a.val<b.val;
//}
//改成
//int cmp1(node a,node b){
// return a.val==b.val?a.pos<b.pos:a.val<b.val;
//}
//提交，等了10分钟，AC了，没想到网站卡壳被我遇到了。
//快排是不稳定排序，这个坑不是第一次遇到了。2019-5-21
#include <cstdio>
#include <algorithm>
#include <cstring>
#define maxn 50100
using namespace std;
struct node{
int pos,val;
}a[maxn],b[maxn];
int cmp1(node a,node b){
return a.val==b.val?a.pos<b.pos:a.val<b.val;
}
int cmp2(node a,node b){
return a.pos<b.pos;
}
int main(){
int T,n,i,cnt;
scanf("%d",&T);
while(T--){
memset(a,0,sizeof(a)),memset(b,0,sizeof(b));
scanf("%d",&n);
for(i=1;i<=n;i++)scanf("%d",&a[i].val),a[i].pos=i;
sort(a+1,a+1+n,cmp1);
b[1]=a[1],cnt=1;
for(i=2;i<=n;i++)
if(a[i].val!=a[i-1].val)b[++cnt]=a[i];//去重
sort(b+1,b+1+cnt,cmp2);
printf("%d",b[1].val);
for(i=2;i<=cnt;i++)printf(" %d",b[i].val);
printf("\n");
}
return 0;
}
5.bzoj 1432 //1432: [ZJOI2009]Function //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1432
//1432: [ZJOI2009]Function
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1432
//直线两两相交，但没有三条直线交在一起，
//读题多遍，有疑问，分层的依据是什么，没弄明白啊
//题意理解很困难。2019-5-22
//经画图实践即观测不同图形形态，对题中的层作如下定义
//交点为层的边界，同一层，只能在交点的一侧，且两层之间的线段不能重合。
//https://blog.csdn.net/flaze_/article/details/52886497此文思路代码不错
//嗯………………考虑一下，显然可以上下翻转，于是对于求f[k]( k>(n/2))的，等效于求f[n-k]，所以只要保证对于某个n，前(n>>1)层的值最小即可
//考虑一下…………
//当n==1，答案是1
//当n==2，答案是2
//之后每次增加一条线，都不会影响前k层的值，但是因为题目说任意两条直线不平行，于是一定会有交点
//我们从上往下增加直线，当已经放了k-1条直线，在放第ｋ条直线的时候，
//如果直线k和直线k-1的交点在k-1原有交点的右边，坑定会影响前面的答案，不是最优，所以强行让交点在倒数第二个，不会影响f[k-1]，而f[k]显然会因此增加了2（自己画一画）
//n=1输出1
//n>1,有n层，若k<=n/2 输出2*k，反之输出(n-k+1)*2层
//说实话，该题，还是很难想出的。2019-5-22
//样例通过，提交AC。2019-5-22
#include <stdio.h>
int main(){
int n,k;
scanf("%d%d",&n,&k);
if(n==1)printf("1\n");
else printf("%d\n",k<=n-k+1?2*k:(n-k+1)*2);
return 0;
}
6.bzoj 3098 //3098: Hash Killer II //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=3098
//3098: Hash Killer II
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=3098
//如果一个房间里有23个或23个以上的人，那么至少有两个人的生日相同的概率要大于50%。
//https://baike.baidu.com/item/%E7%94%9F%E6%97%A5%E6%82%96%E8%AE%BA/2715290?fr=aladdin
//生日悖论推导如下
#include <stdio.h>
int main(){
int i=1;
double ans=1;//ans所有人生日都不相同的概率
while(1){
ans*=(365-i+1)*1.0/365;
if(ans<0.5)break;
i++;
}
printf("i=%d ans=%lf\n",i,ans);
return 0;
}
//23的得出，可以看作是23=sqrt(365)
//https://www.cnblogs.com/wuminyan/p/5211715.html此文思路不错
//这里取模的数是10^9+7,所以只需要生成sqrt（10^9+7）≈100000的数就会出现冲突
//提交AC。2019-5-22 22:07
//此题关于l的取值，还是有些疑惑，日后要再看看。2019-5-22 22:16
#include <stdio.h>
#include <stdlib.h>
int main(){
int i;
printf("100000 10\n");
for(i=1;i<=100000;i++)printf("%c",rand()%26+'a');
return 0;
}
7.bzoj 1800 //1800: [Ahoi2009]fly 飞行棋 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1800
//1800: [Ahoi2009]fly 飞行棋
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1800
//画图进行研究，发现，只要知道 直径 数量m,C(m,2)即可得到长方形数量
//采用前缀和的方式，存储点间圆弧长度。
//若两点差值为圆周长一半，此两点即为直径。
//样例通过，提交Wrong_Answer，表示希望越大，失望越大。
//重读源代码，发现测试代码竟然没删除，if((sum[j]-sum[i])*2==sum[n])m++,printf("i=%d j=%d\n",i,j);
//天啊，功亏一篑，删除测试代码，提交AC.2019-5-23
//提交之前，还是要测试一遍，对对输出。
//该题独立完成，甚是欣慰。2019-5-23
//思维量比较大。
#include <stdio.h>
int n,sum[25],m=0;
int main(){
int i,j;
scanf("%d",&n);
sum[0]=0;
for(i=1;i<=n;i++)scanf("%d",&sum[i]),sum[i]+=sum[i-1];//前缀和
for(i=0;i<n;i++)//此处写成 for(i=0;i<=n;i++)
for(j=i+1;j<n;j++)//此处写成 for(j=i+1;j<=n;j++)
if((sum[j]-sum[i])*2==sum[n])m++;//统计直径 数量
printf("%d\n",m*(m-1)/2);//任选2条直径。
return 0;
}
8.bzoj 1218 //1218: [HNOI2003]激光炸弹 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1218
//读完题意，第一点，求子矩阵之和
//第二点，正方形占据的位置，可以是目标之间的空隙。
//考虑到数组计算问题，所有坐标值都自加1
//R=1,只能炸1个点
//R=2,只能炸(2-1)^2=1个点
//R=3,只能炸(3-1)^2=4个点
//R=4, 只能炸(4-1)^2=9个点
//故R=1时，使其变成R=2即可。
//样例通过，提交Wrong_Answer，信心慢慢，可惜。
//想了想，提供的数据很可能是长方形，而非正方形，故需补齐成正方形
//修改代码，提交Wrong_Answer
//发现
//for(i=R;i<=S;i++)//此处写成for(i=R;i<=row;i++)
// for(j=R;j<=S;j++)//此处写成 for(j=R;j<=col;j++)
//修改，提交 Wrong_Answer
//在https://www.luogu.org/recordnew/show/19293278提交，30分，测试点1-2,6-10WA
//在洛谷里测试了一遍，发现是边界的问题。
//计算矩阵和的算法，没有问题。
//通过对比代码发先，题目叙述有问题，“若目标位于爆破正方形的边上，该目标将不会被摧毁。”
//还不如这样说R=1,有1个元素会被炸;R=2，有4个元素会被炸;R=3,有9个元素会被炸
//该题最后问题竟然出现在题意的理解问题上。2019-5-23 17:19
//在洛谷里提交通过，继续在bzoj上提交，AC。201-5-23 17:23
//该题的问题，题意不明，可惜，离成功只差一步，这一步，靠自己是难以 跨越的，还请出题人多多费心。
//做完该题，再做 bzoj 2462 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2462会更有收获
//bzoj 2462 题解 详见https://blog.csdn.net/mrcrack/article/details/90228821
#include <stdio.h>
#include <string.h>
int sum[5010][5010];//sum[i][j]表示[1][1]-[i][j]之间的矩阵和
int max(int a,int b){
return a>b?a:b;
}
int main(){
int i,j,n,R,row=0,col=0,x,y,v,ans=0,S=0;//此处写成 int i,j,n,R,row,col,x,y,v;
memset(sum,0,sizeof(sum));
scanf("%d%d",&n,&R);
for(i=1;i<=n;i++)scanf("%d%d%d",&x,&y,&v),x++,y++,sum[x][y]=v,row=max(row,x),col=max(col,y);
S=max(row,col);//漏考虑
for(i=1;i<=S;i++)//此处写成 for(i=1;i<=row;i++)
for(j=1;j<=S;j++)//此处写成 for(j=1;j<=col;j++)
sum[i][j]+=sum[i][j-1];//行计算
for(j=1;j<=S;j++)//此处写成 for(j=1;j<=col;j++)
for(i=1;i<=S;i++)//此处写成 for(i=1;i<=row;i++)
sum[i][j]+=sum[i-1][j];//列计算
for(i=R;i<=S;i++)//此处写成for(i=R;i<=row;i++)
for(j=R;j<=S;j++)//此处写成 for(j=R;j<=col;j++)
ans=max(ans,sum[i][j]-sum[i][j-R]-sum[i-R][j]+sum[i-R][j-R]);
printf("%d\n",ans);
return 0;
}
9.bzoj 4001 //4001: [TJOI2015]概率论 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=4001
//4001: [TJOI2015]概率论
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=4001
//不得不说，树，确实是薄弱环节。
//该题与二叉树的形态有关。
//估计涉及 卡特兰数
//看了 判断 树是否同构的代码，得出上述。
//叶结点，有 2个极端，一是 链式，只有各一个叶结点
//另一个是 满二叉树，若达不到，接近满二叉树。对应最多节点.2019-5-23
//1<=N<=10^9，不用担心，因树中查找，算法的时间复杂度是O(logN)
//回忆了 卡特兰数，还是想了起来，
//f(n)=f(n-1)f(0)+f(n-2)f(1)+f(n-3)f(2)+...+f(1)f(n-2)+f(0)f(n-1)
//根需占用1个节点，f(0)=1,f(1)=2
//期望的分母解决了，就是 卡特兰数
//分子 如何处理？
//大概就到这个程度了，准备开始翻看他人代码。
//翻看全网，发现思考方向正确，确实在分析 叶节点数量 需要学习 他人
//http://blog.miskcoo.com/2015/04/bzoj-4001?replytocom=44067此文思路不错
//卡特兰数 由 递推式 推导 公式，
//详见《数据结构(C语言版)》(严蔚敏 吴伟民）P152-P154 非常详细具体，超越网络中大部分的推导。
//学习上述几页内容的感慨，写书的作者在写书时，早已是熟练者的程度，
//书因篇幅的限制，在推导过程中，往往会略过作者认为的简单部分，
//但对于初学者来说，往往是难以逾越的障碍，很多初学者因此望而却步，
//能弄明白的，往往要花相当的精力及耐力。当然，就此也就相当部分被挡在了门外。2019-5-25 21:00
//找遍全网，从生成函数进行推导，也就此文写得好https://blog.csdn.net/Estia_/article/details/87915439摘抄如下
应该能比较容易的推导出(2n-2)!/(n-1)!/(n-1)! 2019-5-26 22:13
//2019-5-26 21:08
//提交AC。2019-5-26 22:18
#include <stdio.h>
int main(){
long long n;
scanf("%lld",&n);
printf("%.9lf\n",1.0*n*(n+1)/2/(2*n-1));
}
10.bzoj 2760 //2760: [JLOI2011]小A的烦恼 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2760
//2760: [JLOI2011]小A的烦恼
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2760
//","号数+1即为每行元素个数
//一个文件一个文件的处理
//开一个字符数组output[105][205*100],用来存储输出字符串
//每个文件处理时，每行都统计","的个数
//找出同个文件中最大","个数，其它行，不足","，补齐
//样例通过，提交Output_Limit_Exceed
//对比代码，发现网络中全是STL模版编写
//看完代码，觉得没有问题，突然想到，应是input初始化问题
//memset(input,0,sizeof(input)); //漏了此句，查了好久
//修改，提交AC。2019-5-27
#include <stdio.h>
#include <string.h>
char file[105],input[105][205],output[105][205*100];
int dot[105];
int max(int a,int b){
return a>b?a:b;
}
int main(){
int t,m,i,j,len,cnt,row=0;
scanf("%d",&t);
while(t--){
memset(input,0,sizeof(input)); //漏了此句，查了好久
memset(dot,0,sizeof(dot)),cnt=0;
scanf("%d%s",&m,file);
for(i=1;i<=m;i++){
scanf("%s",input[i]);
len=strlen(input[i]);
for(j=0;j<len;j++)
if(input[i][j]==',')
dot[i]++;
cnt=max(cnt,dot[i]);//cnt用来统计当前文件最多","数
}
for(j=1;j<=cnt;j++)strcat(file,",");
for(i=1;i<=100;i++){//此处写成 for(i=1;i<=m;i++)
for(j=dot[i];j<cnt;j++)
strcat(input[i],",");
}
if(t){//不是最后一个文件，补上文件间间隔","
strcat(file,",");
for(i=1;i<=100;i++)//此处写成 for(i=1;i<=m;i++)
strcat(input[i],",");
}
strcat(output[0],file);
for(i=1;i<=100;i++)strcat(output[i],input[i]);//此处写成 for(i=1;i<=m;i++)strcat(output[i],input[i]);
row=max(row,m);
}
for(i=0;i<=row;i++)printf("%s\n",output[i]);
return 0;
}
11.bzoj 1088 //1088: [SCOI2005]扫雷Mine //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1088
//1088: [SCOI2005]扫雷Mine
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1088
//题中少了图，建议可在https://www.luogu.org/problemnew/show/P2327看到缺失的图
//借助图形与样例，很快弄懂题意，不过，却没什么思路。
//有一点感觉，就是int可能溢出，采用long long稳妥
//此文思路代码都写得不错https://cloud.tencent.com/developer/article/1087793摘抄如下
//感叹，直接动手模拟，敢于尝试，这种能力还是比较欠缺。2019-5-27
//样例通过，提交AC。2019-5-27
#include <stdio.h>
#define maxn 10010
int a[maxn],b[maxn];
int main(){
int n,i,flag,cnt=0;
scanf("%d",&n);
a[0]=b[0]=0;//初始化
for(i=1;i<=n;i++)scanf("%d",&b[i]);
for(a[1]=0;a[1]<=1;a[1]++){//a[1]=0无雷，a[1]=1有雷
flag=1;
for(i=2;i<=n;i++){
a[i]=b[i-1]-a[i-1]-a[i-2];//注b[i-1]=a[i-2]+a[i-1]+a[i]这样比较好理解
if(!(a[i]==0||a[i]==1)){//即不是0，也不是1，即非法的情况
flag=0;
break;
}
if(i==n&&!(a[i]==b[i]-a[i-1]))flag=0;//边界还需判定
}
if(flag)cnt++;
}
printf("%d\n",cnt);
return 0;
}
12.bzoj 1034 //1034: [ZJOI2008]泡泡堂BNB //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1034
//1034: [ZJOI2008]泡泡堂BNB
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1034
//田忌赛马
//快排，升序，降序，进行比较
//样例通过，提交Wrong_Answer。2019-5-27 21:05
//此文代码思路都很不错 https://blog.csdn.net/nike0good/article/details/8578176摘抄如下
//田忌赛马加强版。
//尽可能让最弱的赢，最强的赢，都不行则最弱打最强
//测试样例
//输入
//4
//1 2 6 7
//1 2 4 7
//输出
//6 2
//最差情况：
//容易发现2个队伍的总分一定(2-0,1-1,0-2) 即3场比分3:3，总得分6 2*n即是这样得出
//所以ans2=2n-B队最高分
//样例通过，提交AC。2019-5-27 21:37
//重读了一遍田忌赛马的课文，发现课文上升不到算法高度。2019-5-27 21:40
#include <cstdio>
#include <algorithm>
#define maxn 100100
using namespace std;
int a[maxn],b[maxn];
int cal(int *a,int *b,int n){
int i=1,j=1,k=n,m=n,ans=0;
while(i<=k){
if(a[i]>b[j])ans+=2,i++,j++;
else if(a[k]>b[m])ans+=2,k--,m--;
else ans+=a[i]==b[m],i++,m--;
}
return ans;
}
int main(){
int n,i,ans1,ans2;
scanf("%d",&n);
for(i=1;i<=n;i++)scanf("%d",&a[i]);
for(i=1;i<=n;i++)scanf("%d",&b[i]);
sort(a+1,a+1+n);
sort(b+1,b+1+n);
ans1=cal(a,b,n);
ans2=2*n-cal(b,a,n);
printf("%d %d\n",ans1,ans2);
return 0;
}
13.bzoj 1207 //1207: [HNOI2004]打鼹鼠 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1207
//1207: [HNOI2004]打鼹鼠
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1207
//象数字三角形，但又不是
//f[t][x][y]动归，但发现n（n<=1000）, m（m<=10000），空间溢出
//此文思路不错https://www.cnblogs.com/five20/p/9021689.html摘抄如下
//本题吐槽数据太水，直接裸的最长上升子序列（O(n2)
//）也能过。。。
// 先贪心思路，由于给出的鼹鼠出现时间已经单调递增，所以一开始肯定是从第1
//个鼹鼠位置出发，因为就算从后面出发最多能到达的位置1都能到达。
//再考虑到从一个点到下一个点时间就是两点的曼哈顿距离，而且可以停留，
//那么限制条件就是abs(x[i]?x[j])+abs(y[i]?y[j])<=t[i]?t[j]。
// 此时容易想到本题实际在求类似最长上升子序列（只不过多了限制条件）。
// 所以定义状态f[i]
//表示前i个鼹鼠中最多打的鼹鼠个数，则不难想到状态转移方程：
//（前提：在满足限制条件下）f[i]=Max(f[j])+1,j∈[1,i)。
// 但是数据104
//，O(n2)显然超时（事实是数据太水～～直接过）。于是考虑优化，由于有限制条件，
//普通的二分优化就行不通了，此时一个贪心的玄学优化是设置一个mx[i]表示前i个数中最长上升子序列的长度，
//每次更新时同步更新mx，那么在限制条件前可以加一条可行性剪枝if(mx[j]+1<=f[i])break;
//（枚举j时倒序枚举，表示当前j个中的最长上升子序列长度+1都不大于f[i]，就没必要继续往前找了，因为往前长度显然会更小），
//每次更新mx也很简单：mx[i]=max(mx[i?1],f[i])。
// 最后输出目标状态ans
//就OK了。
//样例通过，提交AC。2019-5-28
//总体来讲，建模的能力还需加强。
#include <stdio.h>
#define maxm 10010
int mx[maxm],t[maxm],x[maxm],y[maxm],f[maxm],ans;//mx[i] 1-i之间最大长度
int max(int a,int b){
return a>b?a:b;
}
int abs(int a){
if(a<0)a=-a;
return a;
}
int main(){
int n,m,i,j;
scanf("%d%d",&n,&m);
for(i=1;i<=m;i++)scanf("%d%d%d",&t[i],&x[i],&y[i]);
mx[0]=0,ans=0;
for(i=1;i<=m;i++){//最长上升子序列
f[i]=1;
for(j=i-1;j>=1;j--)//逆向查找
if(mx[j]+1<=f[i])break;//剪枝
else if(f[j]+1>f[i]&&abs(x[i]-x[j])+abs(y[i]-y[j])<=t[i]-t[j])
f[i]=f[j]+1;
mx[i]=max(f[i],mx[i-1]),ans=max(ans,f[i]);//写错位置，将其写在上面else if中
}
printf("%d\n",ans);
return 0;
}
14.bzoj 2748 //2748: [HAOI2012]音量调节 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2748
//2748: [HAOI2012]音量调节
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2748
//改变音量，可加可减。50个点，2^50=1024^5超时
//深搜，剪枝，可以试试。
//过程中，若不在0-maxLevel范畴，则剪枝
//样例解释如下
//5-5+3+7=10
//样例通过，提交Time_Limit_Exceed
//深搜代码如下，虽然无法AC。2019-5-28
//此文思路代码都不错https://www.cnblogs.com/five20/p/9071611.html摘抄如下
//这是一次考试的防爆0的DP题目，5min钟就能刚完。
// 定义状态f[i][j]
//表示到了第j次操作，能否达到i的音量。初始状态f[st][0]=1
// 状态转移方程显而易见：if(f[i][j?1]==1)f[i+a[j]][j]=f[i?a[j]][j]=1
// 然后倒序找满足f[i][n]==1的最大的i。
//有01背包，布尔背包的味道
//样例通过，提交AC。2019-5-28
#include <stdio.h>
#include <string.h>
int c[1010],f[55][1010];
int main(){
int i,j,n,b,m;
scanf("%d%d%d",&n,&b,&m);
for(i=1;i<=n;i++)scanf("%d",&c[i]);
memset(f,0,sizeof(f)),f[0][b]=1;
for(i=1;i<=n;i++)//第一重循环 是 个数
for(j=0;j<=m;j++)
if(f[i-1][j]){
if(j-c[i]>=0)f[i][j-c[i]]=1;
if(j+c[i]<=m)f[i][j+c[i]]=1;
}
for(j=m;j>=0;j--)
if(f[n][j])break;
printf("%d\n",j);//若不在0-m范围内，j=-1
return 0;
}
//2748: [HAOI2012]音量调节
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2748
//改变音量，可加可减。50个点，2^50=1024^5超时
//深搜，剪枝，可以试试。
//过程中，若不在0-maxLevel范畴，则剪枝
//样例解释如下
//5-5+3+7=10
//样例通过，提交Time_Limit_Exceed
//深搜代码如下，虽然无法AC。2019-5-28
#include <stdio.h>
int n,b,m,c[55],ans=-1;
int max(int a,int b){
return a>b?a:b;
}
void dfs(int step,int level){
if(level<0||level>m)return;
if(step==n+1){
ans=max(ans,level);
return;
}
dfs(step+1,level+c[step+1]);
dfs(step+1,level-c[step+1]);
}
int main(){
int i;
scanf("%d%d%d",&n,&b,&m);
for(i=1;i<=n;i++)scanf("%d",&c[i]);
dfs(0,b);
printf("%d\n",ans);
return 0;
}
15.bzoj 4302 //4302: Hdu 5301 Buildings //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=4302
//4302: Hdu 5301 Buildings
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=4302
//对了翻译，发现还是有些地方理解不到位。
//建议看看原题https://vjudge.net/problem/HDU-5301有配图对样例进行说明
//此文题意翻译不错https://blog.csdn.net/yeguxin/article/details/47029057
//此文思路清晰https://blog.csdn.net/Fun_Zero/article/details/47059673摘抄如下
//此文https://blog.csdn.net/controlbear/article/details/52014702代码思路都写得很不错
//算法核心部分，用下图进行理解
//黑色为障碍，红色部分 (n+1)/2即可填充，白色部分 竖铺 max(up,down),横铺 max(left,right),之后在 竖铺 横铺 中选取最小值
//最后，在红色填充，与白色 铺中选取最大值。2019-6-1 8:39
//思路如上图
//样例通过，提交Wrong_Answer
//重读代码，发现 right=m-y+1;//此处写成 right=n-y+1 笔误
//修改，提交Wrong_Answer
//left=y//此处写成left=y+1对比代码，才发现
//修改，提交AC。2019-6-1 9:02
//想到与做到，还是有一定距离。
#include <stdio.h>
void swap(int *a,int *b){
int t;
t=*a,*a=*b,*b=t;
}
int min(int a,int b){
return a<b?a:b;
}
int max(int a,int b){
return a>b?a:b;
}
int main(){
int n,m,x,y,up,down,left,right,ans,h,w,white;
while(scanf("%d%d%d%d",&n,&m,&x,&y)!=EOF){
if(n>m)swap(&n,&m),swap(&x,&y);
ans=(n+1)/2;
if(n==m&&n%2&&x==y&&x==(n+1)/2)ans--;//特判
else{
up=x-1,down=n-x,left=y,right=m-y+1;//此处写成left=y+1对比代码，才发现//此处写成 right=n-y+1 笔误
h=max(up,down),w=min(left,right),white=min(h,w);//白色区域处理
ans=max(white,ans);//白红区域合并
}
printf("%d\n",ans);
}
return 0;
}
16.bzoj 1083 //1083: [SCOI2005]繁忙的都市 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1083
//1083: [SCOI2005]繁忙的都市
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1083
//该题就是一个裸的最小生成树
//准备采用Kruskal算法，因其简单好写。
//样例通过，提交AC。2019-6-1 20:57
#include <cstdio>
#include <algorithm>
using namespace std;
int n,m,f[306];
struct node{
int u,v,c;
}e[306*306/2];
int cmp(node a,node b){
return a.c<b.c;
}
int getf(int u){
return f[u]==u?u:f[u]=getf(f[u]);
}
int merge(int u,int v){
int f1=getf(u),f2=getf(v);
if(f1!=f2){
f[f2]=f1;//左靠
return 1;
}
return 0;
}
int main(){
int i,u,v,c,b=0;
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)f[i]=i;
for(i=1;i<=m;i++)scanf("%d%d%d",&e[i].u,&e[i].v,&e[i].c);
sort(e+1,e+1+m,cmp);
for(i=1;i<=m;i++)
if(merge(e[i].u,e[i].v)){
b++;
if(b==n-1){
printf("%d %d\n",n-1,e[i].c);
break;
}
}
return 0;
}
17.bzoj 2563 //2563: 阿狸和桃子的游戏 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2563
//2563: 阿狸和桃子的游戏
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2563
//最有策略，包括互相干扰吗，尝试了一下，没有对上样例，搜索网络，发现想不到，或者说不会想
//https://www.cnblogs.com/zzyh/p/7640470.html此文思路摘抄如下
//把一条边的边权分给两端的点权，如果一条边的两端被两人选择，那么最后得分相减时会抵消，如果被一个人选择，那么
//这条边的权值都会给他。
//100000/2*10000=5*10^8 int不会溢出
//因权值可能为奇数，故将其扩大1倍，最后处理完后，再将结果除2. 疑虑，除2前 可能为奇数
//样例通过，提交AC。2019-6-2 8:36
#include <cstdio>
#include <algorithm>
using namespace std;
int v[10010];
int cmp(int a,int b){
return a>b;
}
int main(){
int n,m,a,b,c,i,sum1=0,sum2=0;
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)scanf("%d",&v[i]),v[i]*=2;
for(i=1;i<=m;i++)scanf("%d%d%d",&a,&b,&c),v[a]+=c,v[b]+=c;
sort(v+1,v+1+n,cmp);
for(i=1;i<=n;i++)
if(i%2)sum1+=v[i];
else sum2+=v[i];
printf("%d\n",(sum1-sum2)/2);
return 0;
}
18.bzoj 3097 //3097: Hash Killer I //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=3097
//3097: Hash Killer I
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=3097
//第二次遇到该题，这次终于看懂了题干中的代码。觉得写得不错，没什么问题。
//怎么卡，在没有提示的情况下，没有思路。
//打算测试unsigned long long自然溢出，操作系统中的计算器，支持不好
//马上降维，测试unsigned int自然溢出.2019-6-2 15:05
//自然溢出代码如下
#include <stdio.h>
int main(){
unsigned int x;
while(scanf("%u",&x)!=EOF) printf("%u\n",x);
return 0;
}
//win7测试效果如下
4294967295
4294967295
4294967296
0
4294967297
1
4294967298
2
//自然溢出相当于 模 4294967296即2^32 那么 unsigned long long溢出同理。2019-6-2 15:22
//Ubuntu14.04.2 采用 g++ 3097.cpp -o 3097
// ./3097测试效果如下
-1
4294967295
4294967295
4294967295
4294967296
4294967295
4294967297
4294967295
4294967298
4294967295
//对上述结果深表遗憾，在Ubuntu14.04.2 采用 g++ 3097.cpp -o 3097 自然溢出是什么情况？2019-6-2 15:29
//https://blog.csdn.net/regina8023/article/details/43112899此文思路写得不错，摘抄如下
非常神奇的构造题。
首先明白两点：
1.卡hash的关键在于构造两个不同的串对应的hash值相同。
2.爆u64相当于对2^64这个数取模。
如果base是偶数，那么a.........aaa(>64个a)与ba.......aa(a的数量为前面那么串a的数量-1)，这两个串长度相同，hash值相同，显然串是不同的，这样就卡掉了。(通篇看完后，才弄懂此处，ba.......aa中 (‘b’-'a')*base^(前面那么串a的数量-1) 因 base是偶数，故 2^64\(‘b’-'a')*base^(串a的数量-1) mrcrack添加 2019-6-2 22:14 )
如果base是奇数，就比较麻烦了。
看vfk的做法吧：
如果base是奇数的话，现在只考虑a、b两个字母。
a \ b表示a能整除b。（orz 具体数学）
设数学上的函数not(S)表示把字符串S中每个位置的'a'变成'b'，把'b'变成'a'后形成的字符串。比如not("ababaa") = "bababb"
strA . strB代表字符串串联。如"娃" . "哈哈" = "娃哈哈"
|str|表示字符串str的长度。
设字符串序列{orzstr[i]}，orzstr[1] = "a", orzstr[i] = orzstr[i - 1] . not(orzstr[i - 1]) （注，i>=2 mrcrack添加）
那么|orzstr[i]| = |orzstr[i - 1]| * 2。显然这是等比数列，得到：|orzstr[i]| = |orzstr[1]| . 2 ^ (i - 1) = 2 ^ (i - 1)
设hash(str)为str的哈希值。
则：
hash(orzstr[i]) = hash(orzstr[i - 1]) * base ^ |not(orzstr[i - 1])| + hash(not(orzstr[i - 1])) (由题干中 hash算法决定 mrcrack添加)
= hash(orzstr[i - 1]) * base ^ (2 ^ (i - 2)) + hash(not(orzstr[i - 1]))
hash(not(orzstr[i])) = hash(not(orzstr[i - 1])) * base ^ (2 ^ (i - 2)) + hash(orzstr[i - 1]) (对比上式，依次类推 mrcrack添加)
两式相减：
hash(orzstr[i]) - hash(not(orzstr[i]))
= (hash(orzstr[i - 1]) * base ^ (2 ^ (i - 2)) + hash(not(orzstr[i - 1]))) - (hash(not(orzstr[i - 1])) * base ^ (2 ^ (i - 2)) + hash(orzstr[i - 1]))
= (hash(orzstr[i - 1]) - hash(not(orzstr[i - 1]))) * (base ^ (2 ^ (i - 2)) - 1)
这让我们发现，hash(orzstr[i]) - hash(not(orzstr[i]))似乎是个神奇的东西。而我们的目的实际上是要找两个字符串strA, strB使得
hash(strA) % 2^64 = hash(strB) % 2^64
相当与
2^64 \ hash(strA) - hash(strB)
设数列{f[i]}，f[i] = hash(orzstr[i]) - hash(not(orzstr[i]))
这样就有：
f[i] = f[i - 1] * (base ^ (2 ^ (i - 2)) - 1)
还是有点不爽啊……我们再设数列{g[i]}，g[i] = base ^ (2 ^ (i - 1)) - 1
于是能写成：
f[i] = f[i - 1] * g[i - 1]
则f[i] = f[1] * g[1] * g[2] * ... * g[i - 1]
然后发现一个神奇的事情？
base是奇数，则base的任意正整数次方也一定是奇数。所以对于任意的i必有g[i]为偶数，所以2 ^ (i - 1) \ f[i] （即 g[1]=2*k1,g[2]=2*k2,g[3]=2*k3,...,g[i-1]=2*k(i-1),故2^(i-1)\g[1]*g[2]*...*g[i-1] mrcrack添加 2019-6-2 21:54）
问题是不是结束了呢……发现没有……这样的话我们要使2 ^ 64 \ f[i]，至少得让i = 65……然后发现|orzstr[65]|是个天文数字。
发现我们刚才那样分析太坑爹了……
i > 1时有：
g[i] = base ^ (2 ^ (i - 1)) - 1 = (base ^ (2 ^ (i - 2)) - 1) * (base ^ (2 ^ (i - 2)) + 1) = g[i - 1] * 一个偶数
而g[1]显然是偶数吧……
那么4 \ g[2]，8 \ g[3]...
也就是说2 ^ i \ g[i]
所以f[i] 实际上有：
(2 ^ 1) * (2 ^ 2) * (2 ^ 3) * ... * (2 ^ (i - 1)) \ f[i]
2 ^ (i * (i - 1) / 2) \ f[i]
当i取12时，就有66个2了哟！
这就是卡base为奇数时的方法。orzstr[12]和not(orzstr[12])即为所求。
而读入中base既有奇数又有偶数，直接在奇数构造的字符串后面加64个a就可以了。
//感觉该题，属于 出题人 自娱自乐，这种构造怎么想得到。2019-6-2 21:12
//不过，通过该问题学习，对hash的冲突，有了更深刻的了解。
//比较奇怪，该题的输出数据是怎么个形式，因为存在多解。2019-6-3 5:51
//样例通过，提交AC。2019-6-3 8:56
#include <stdio.h>
int str[(1<<12)+65+10];
int main(){
int n,L,i,j,now=1;
printf("%d %d\n",(1<<12)+65,1<<11);//每串 orzstr[12]和not(orzstr[12] 即1<<11,两串拼接,即1<<12,65由来 b+65个a 串中 取b+64个a 和 65个a
str[1]='a';
for(i=1;i<=12;i++){
for(j=1;j<=now;j++)
str[now+j]=str[j]=='a'?'b':'a';
now*=2;
}
for(i=1;i<=(1<<12);i++)printf("%c",str[i]);//打印两串完全相反的串组成的新串。 对应base 为 奇数
for(i=1;i<=65;i++)printf("a");//打印附加的65个a 对应 base 为 偶数
return 0;
}
19.bzoj 1022 //1022: [SHOI2008]小约翰的游戏John //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1022
//1022: [SHOI2008]小约翰的游戏John
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1022
//https://blog.csdn.net/mrcrack/article/details/80562324通过此文学习 Nim游戏 弄明白了 原理
//该题数据，均化成二进制，奇数个1,后手赢，偶数个1,先手赢。
//判断 奇数个1 偶数个1 异或 最合适。
//写着写着，还是觉得上述思路有误，因结果与Nim刚好相反，故思路还是要有大动。2019-6-3 22:16
//分类讨论。此文思路写得不错https://blog.csdn.net/clover_hxy/article/details/53818624
//搜遍网络，速成解析均难以理解，无奈，只好从头开始研习国家集训队的3篇论文。
//https://wenku.baidu.com/view/e119050a326c1eb91a37f111f18583d049640fe2.html?from=search《由感性认识到理性认识——透析一类搏弈游戏的解答过程》2019-6-3 21:38
上文中 若A胜B胜，则有时S胜，有时S负。感到难以理解，A胜，之后B开始时，乙先手，在B中，乙胜，即S负。2019-6-3 22:07
//https://wenku.baidu.com/view/67e4e88583d049649b6658f0.html《解析一类组合游戏》
//https://wenku.baidu.com/view/7cd481e9524de518964b7d1f.html《组合游戏略述——浅谈SG游戏的若干拓展及变形》
//https://blog.csdn.net/bestsort/article/details/88197959SG函数此文介绍得不错 2019-6-6 22:40
//https://blog.csdn.net/strangedbly/article/details/51137432SG函数此文图文并貌。2019-6-7 22:00
摘自https://blog.csdn.net/PoPoQQQ/article/details/40543799
题目大意：反Nim游戏，即取走最后一个的人输
首先状态1：如果所有的堆都是1，那么堆数为偶先手必胜，否则先手必败
然后状态2：如果有两个堆数量相同且不为1，那么后手拥有控场能力，即：
若先手拿走一堆，那么后手可以选择将另一堆留下1个或者全拿走，使这两堆最终只剩1个或0个；
若先手将一堆拿剩一个，那么后手可以选择将另一堆留下一个让先手拿或全拿走，使这两堆最终只剩1个或0个；
若先手将一堆拿走一部分，那么后手可以将另一堆同样拿走一部分，然后同上
状态3：若Xor!=0 那么先手可以先拿走一部分让Xor=0 然后同状态2先手必胜 否则先手必败
※鉴于本人过于沙茶，以上内容仅存在参考价值，无法证明正确性，欢迎指正
于是若所有堆全是1 Xor==0先手必胜 否则后手必胜
若有堆不是1 Xor==0后手必胜 否则先手必胜
摘自https://blog.csdn.net/WerKeyTom_FTD/article/details/48207079
Anti-SG游戏定义
1、决策集合为空的操作者胜。
2、其余规则与SG游戏一致。
SJ定理
对于任意一个Anti-SG游戏，如果定义所有子游戏的SG值为0时游戏结束，先手必胜的条件：
1、游戏的SG值为0且所有子游戏SG值均不超过1。
2、游戏的SG值不为0且至少一个子游戏SG值超过1。
证明
先证明第一个条件：
所有都不超过1，那么显然如果有偶数个1则先手必胜偶数个1即游戏的SG值为0。
再证明第二个条件：
如果有至少两个子游戏SG值超过1，然后其余有偶数个1。那么游戏的SG值就为所有SG值超过1的子游戏的SG值的Nim和。设为X。假设X中1的最高位为第k位，那么一定存在至少一个子游戏的SG值第k位有1，设这个子游戏的SG值为Y。那么我们可以把这个子游戏变为YxorX
。可以知道YxorX<Y，因为会消去第k位上的1，k位之后不变，那么相当于减少了一个2^(k−1)，前k-1位不论如何变，最多增加2^(k−1)−1。根据SG函数的转移可以得知一定能转移到一个决策其SG值为YxorX
。接下来游戏的SG值便变为0，而此时至少存在一个子游戏的SG值超过1。此条件下先手必胜。
如果有至少两个子游戏的SG值超过1，然后其余有奇数个1。我们设所有SG值超过1的子游戏的SG值的Nim和为X。那么游戏的SG值为Xxor1
。因为不为0，那么现在如果X=0，我们去掉一个1即可。于是变为游戏的SG值为0，而此时至少存在两个子游戏的SG值超过1。如果X不为0，那么最高位的1一定不在第一位。于是和1没有关系，像上面一样处理即可。
如果只有一个子游戏的SG值超过1，然后其余有偶数个1。显然我们可以拿将那个超过1的子游戏转移到一个决策使SG值为1，就能造成奇数个1。
如果只有一个子游戏的SG值超过1，然后其余有奇数个1。直接将那个超过1的子游戏转移到一个决策使SG值为0，就能造成奇数个1。
这样证明了先手必胜局面有至少一个后继局面为先手必败局面。
接下来需证明每一个先手必败局面都只能转移到先手必胜局面。
先来证明第二个条件：
所有SG值不超过1，那么有奇数个1先手就必败，无论如何都只能转移到有偶数个1的情况。有奇数个1即游戏的SG值不为0。
第一个条件：
如果有至少两个子游戏的SG值超过1，那么无论如何都会转移到一个先手必胜局面（有至少一个SG值超过1，游戏的SG值不为0），因为无法造成游戏的SG值继续为0。
如果只有一个子游戏的SG值超过1，这是不可能发生的，不可能存在一个局面只有一个子游戏SG值超过1但异或和为0的，一定有至少两个子游戏SG值超过1。
证毕。
//样例通过，提交AC.2019-8-4 17:59
#include <stdio.h>
int main(){
int t,n,x,ans,flag;
scanf("%d",&t);
while(t--){
ans=0,flag=0;
scanf("%d",&n);
while(n--){
scanf("%d",&x);
ans^=x;
if(x>1)flag=1;
}
if(ans==0&&flag==0)printf("John\n");
else if(ans&&flag)printf("John\n");
else printf("Brother\n");
}
return 0;
}
20.bzoj 2465 //2465: [中山市选2009]小球 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2465
//2465: [中山市选2009]小球
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2465
//准备上背包，再次读题，发现没这么复杂
//快排，自大到小，再进行选择
//前缀和，能快速计算，马上加上该功能。
//样例通过，提交Wrong_Answer
//这么简单的题，怎么可能。2019-6-5
//重新读题，发现误解提议，原以为是球装最多球的瓶 及其 分数
//读完题目，发现是用上所有瓶子，装的球数，以及分数。
//样例通过，提交Wrong_Answer。2019-6-5 19:25
//重读代码，发现if(j==m+1)break;//此处写成if(j==n+1)break;昏招
//修改，提交Wrong_Answer
//对比他人代码，发现思路没有问题
//看了数据范围，0 <= c <= 200，没想到c可以等于0,这点没考虑到,在读取数据时，马上将该组数据丢弃
//又看了数据范围，0<=m<=200，发现m可以等于0,马上进行特判。
//提交Wrong_Answer，2019-6-5 21:16
//while(scanf("%d%d",&n,&m)&&(n!=0||m!=0)){//此处写成while(scanf("%d%d",&n,&m)&&!(n==0&&m==0))
//修改，提交Wrong_Answer
//https://blog.csdn.net/Mmh2000/article/details/76599914翻到这里才发现，是if与while的问题
//决定重写代码,样例通过，提交AC。
//该题最大的难点，在于c可以等于0,并且可以有一串的c等于0.
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
int a[210];
struct node{
int c,q;
}b[210];
int cmp(int a,int b){
return a>b;
}
int cmp2(node a,node b){
return a.q>b.q;
}
int main(){
int n,m,i,j,ans1,ans2;
while(scanf("%d%d",&n,&m)&&n&&m){
for(i=1;i<=n;i++)scanf("%d",&a[i]);
sort(a+1,a+1+n,cmp);
for(i=1;i<=m;i++) scanf("%d%d",&b[i].c,&b[i].q);
sort(b+1,b+1+m,cmp2);
ans1=ans2=0,i=1,j=1;
while(i<=n){
while(!b[j].c){//此处是难点，写成if错，写成while对。2019-6-5 22:30
j++;
if(j==m+1)break;
}
if(b[j].q>=a[i])ans1++,ans2+=a[i],b[j].c--;
i++;//漏了此句
}
printf("%d %d\n",ans1,ans2);
}
return 0;
}
//2465: [中山市选2009]小球
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2465
//准备上背包，再次读题，发现没这么复杂
//快排，自大到小，再进行选择
//前缀和，能快速计算，马上加上该功能。
//样例通过，提交Wrong_Answer
//这么简单的题，怎么可能。2019-6-5
//重新读题，发现误解提议，原以为是球装最多球的瓶 及其 分数
//读完题目，发现是用上所有瓶子，装的球数，以及分数。
//样例通过，提交Wrong_Answer。2019-6-5 19:25
//重读代码，发现if(j==m+1)break;//此处写成if(j==n+1)break;昏招
//修改，提交Wrong_Answer
//对比他人代码，发现思路没有问题
//看了数据范围，0 <= c <= 200，没想到c可以等于0,这点没考虑到,在读取数据时，马上将该组数据丢弃
//又看了数据范围，0<=m<=200，发现m可以等于0,马上进行特判。
//提交Wrong_Answer，2019-6-5 21:16
//while(scanf("%d%d",&n,&m)&&(n!=0||m!=0)){//此处写成while(scanf("%d%d",&n,&m)&&!(n==0&&m==0))
//修改，提交Wrong_Answer
//https://blog.csdn.net/Mmh2000/article/details/76599914翻到这里才发现，是if与while的问题
//经过一番修改，自己的代码也AC了。2019-6-5 22:44
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
int a[210];
struct node{
int c,q;
}b[210];
int cmp(int a,int b){
return a>b;
}
int cmp2(node a,node b){
return a.q>b.q;
}
int main(){
int n,m,i,j,ans1,ans2,c,q;
while(scanf("%d%d",&n,&m)&&n&&m){
for(i=1;i<=n;i++)scanf("%d",&a[i]);
sort(a+1,a+1+n,cmp);
j=0;
for(i=1;i<=m;i++){
scanf("%d%d",&c,&q);
if(c>0)j++,b[j].c=c,b[j].q=q;
}
m=j;
if(m==0){//特判
printf("0 0\n");
continue;
}
sort(b+1,b+1+m,cmp2);
ans1=ans2=0,i=1,j=1;
while(i<=n){
if(a[i]<=b[j].q){
ans1++,b[j].c--,ans2+=a[i];
if(b[j].c==0){//此处写成if(b[j].c==0)
j++;
if(j==m+1)break;//此处写成if(j==n+1)break;昏招
}
}
i++;//漏了此句
}
printf("%d %d\n",ans1,ans2);
}
return 0;
}
21.bzoj 2122 //2122: 工作评估 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2122
//2122: 工作评估
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2122
//啥算吗，估计要超时
//计算50000*50000=2.5*10^9超时无疑。
//样例通过，提交Time_Limit_Exceed
//以下为超时代码。2019-6-6
//以下为卡时间，恰好AC代码，发现减少变量，也能缩减一点点时间，分块算法 等百题之后再来吧。2019-6-6
//此文值得一看https://blog.csdn.net/qq_36797743/article/details/74626217
#include <stdio.h>
#define maxn 50010
int D[maxn],L[maxn];
int min(int a,int b){
return a<b?a:b;
}
int max(int a,int b){
return a>b?a:b;
}
int main(){
int n,m,i,left,right,ans,X0,X;
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)scanf("%d",&D[i]);
for(i=1;i<=n;i++)scanf("%d",&L[i]);
while(m--){
scanf("%d%d%d",&left,&right,&X0);
ans=0,X=X0;
for(i=left;i<=right;i++)X=min(X+D[i],L[i]),X=max(X,X0),ans=max(X,ans);//此处写成 for(i=left;i<=right;i++)X[i]=X[i-1]+D[i],X[i]=min(X[i],L[i]),X[i]=max(X[i],X[0]),ans=max(X[i],ans);
printf("%d\n",ans);
}
return 0;
}
//2122: 工作评估
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2122
//啥算吗，估计要超时
//计算50000*50000=2.5*10^9超时无疑。
//样例通过，提交Time_Limit_Exceed
//以下为超时代码。2019-6-6
//
#include <stdio.h>
#define maxn 50010
int D[maxn],L[maxn],X[maxn];
int min(int a,int b){
return a<b?a:b;
}
int max(int a,int b){
return a>b?a:b;
}
int main(){
int n,m,i,j,left,right,ans;
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)scanf("%d",&D[i]);
for(i=1;i<=n;i++)scanf("%d",&L[i]);
while(m--){
ans=0;
scanf("%d%d%d",&left,&right,&X[0]);
X[left-1]=X[0];//初始化,该句比较关键
for(i=left;i<=right;i++)X[i]=X[i-1]+D[i],X[i]=min(X[i],L[i]),X[i]=max(X[i],X[0]),ans=max(X[i],ans);
printf("%d\n",ans);
}
return 0;
}
22.bzoj 1197 //1197: [HNOI2006]花仙子的魔法 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1197
//1197: [HNOI2006]花仙子的魔法
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1197
//对样例表示奇怪，不是000,001,010,011,100,101,110,111,共计8种吗
//输出结果，如何跟n维扯上关系。
//看了https://www.luogu.org/problemnew/solution/P4176中的图，对样例有了点感觉
//尝试画图(如下)，弄明白了样例，可喜可贺。2019-6-8 20:56
//想了想，m次魔法，极限类型为2^m次。
//弄懂题意，但是该题难想。
//https://blog.csdn.net/sunshinezff/article/details/48134717此文思路不错，摘抄如下
//简化一下这个题就是让你求一个n维的类似球一样的东西能把一个n维的空间最多分成几个部分。
//然后我们考虑最后一个球。
//两个n维物体的交一定是一个n-1维物体。第m个球最多能和m-1个球相交。这就相当于用m个n-1维的物体去划分一个n-1
//维空间。显然这个之前已经算出来了。
//设 f[i][j]为用j个i维物体去划分一个i维空间，可以划分成最多几个部分。
//所以dp方程就是 f[i][j]=f[i-1][j]+f[i-1][j-1];
//https://www.cnblogs.com/VisJiao/p/BZOJ1197.html此文边界写得不错，推导也值得一看
//记录ans[k][i]表示k维空间中的i个球最多将空间划分为几部分。
//一维情况非常简单，ans[1][i]=2i。一维上的球就是两个点。
//接下来考虑二维情况。假设现在有x个圆，要加入一个新圆。这个新圆可以与之前的每个圆相交，共有2x个交点，圆弧被分成2x份。
//这2x个圆弧都将其所在部分一分为二，也就增加了2x个部分。ans[2][i]=i2−i+2。
//然后考虑三维情况。现在有x个球，要加入一个新球。这个新球可以与之前的每个球相交，每两个球会交出一个圆，那么新球
//的球面就被x个圆划分成ans[2][x]个曲面。这ans[2][x]个曲面都将其所在部分一分为二，
//也就增加了ans[2][x]个部分。ans[3][x+1]=ans[3][x]+ans[2][x]。
//...于是我们发现了一些规律。一个k维球会与i个k维球交出i个k−1维球，
//会增加ans[k−1][i]个部分，即ans[k][i+1]=ans[k][i]+ans[k−1][i]。有初值ans[k][0]=1，那么我们可以递推得到ans[n][m]。
//思维要求高，比较难想到
//样例通过，提交Wrong_Answer
//一查，for(i=1;i<=m;i++)f[1][i]=2*i;//此处写成for(i=1;i<=m;i++)f[1][i]=2*m;
//提交AC。2019-6-8 23:38
#include <stdio.h>
#define LL long long
LL f[20][105];
int main(){
int n,m,i,j;
scanf("%d%d",&m,&n);//此处写成scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)f[i][0]=1;//0次魔法，对应值为1
for(i=1;i<=m;i++)f[1][i]=2*i;//此处写成for(i=1;i<=m;i++)f[1][i]=2*m;//1维空间，对应花的种类
for(i=2;i<=n;i++)
for(j=1;j<=m;j++)
f[i][j]=f[i][j-1]+f[i-1][j-1];
printf("%lld\n",f[n][m]);
return 0;
}
23.bzoj 1295 //1295: [SCOI2009]最长距离 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1295
//1295: [SCOI2009]最长距离
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1295
//通过样例1弄懂题意,最远距离为sqrt(2)，求最远的两点之间的直线距离，看了数据范围，20分容易得。
//通过样例3弄懂了，取出1个障碍的情况，将右下角的1去除，最远距离为2×sqrt(2)
//样例2,最远距离为sqrt(2^2+3^2)
//题意明白后，怎么算。没所有思路
//此文思路写得不错https://www.cnblogs.com/mjtcn/p/8866936.html摘抄如下
//每次选两个点，spfa，到每个点需要搬多少石头，再枚举两个点，
//判断是否可以在搬得石头的个数小于t的情况下，走到，取最大值。
//此文代码写得不错，值得一读http://hzwer.com/1318.html
//spfa
//h=t=1,q[t].r=rr,q[t].c=cc,t++,vis[rr][cc]=1,dis[rr][cc]=map[rr][cc];//漏了此句dis[rr][cc]=map[rr][cc]
//if(dis[i][j]<=T){//此处写成if(dis[rr][cc]<=T)
//样例通过，提交AC。2019-6-9 21:07
#include <stdio.h>
#include <string.h>
#include <math.h>
#define maxn 33
int map[maxn][maxn],dis[maxn][maxn],n,m,T,vis[maxn][maxn],ans=0;
int next[][2]={{-1,0},{1,0},{0,-1},{0,1}};//上 下 左 右
char str[maxn];
struct node{
int r,c;
}q[maxn*maxn*100];
void getans(int rr,int cc){
int i,j,t;
for(i=rr;i<=n;i++)
for(j=1;j<=m;j++)//此处写成for(j=1;j<=n;j++)
if(dis[i][j]<=T){//此处写成if(dis[rr][cc]<=T)
t=(rr-i)*(rr-i)+(cc-j)*(cc-j);
if(t>ans)ans=t;
}
}
void spfa(int rr,int cc){
int h,t,nowr,nowc,nr,nc,k;
memset(vis,0,sizeof(vis)),memset(dis,127,sizeof(dis));
h=t=1,q[t].r=rr,q[t].c=cc,t++,vis[rr][cc]=1,dis[rr][cc]=map[rr][cc];//漏了此句dis[rr][cc]=map[rr][cc]
while(h<t){
nowr=q[h].r,nowc=q[h].c;
for(k=0;k<4;k++){
nr=nowr+next[k][0],nc=nowc+next[k][1];
if(1<=nr&&nr<=n&&1<=nc&&nc<=m&&dis[nowr][nowc]+map[nr][nc]<dis[nr][nc]){
dis[nr][nc]=dis[nowr][nowc]+map[nr][nc];
if(!vis[nr][nc])vis[nr][nc]=1,q[t].r=nr,q[t].c=nc,t++;
}
}
vis[nowr][nowc]=0;
h++;
}
getans(rr,cc);
}
int main(){
int i,j;
scanf("%d%d%d",&n,&m,&T);
for(i=1;i<=n;i++){
scanf("%s",str+1);
for(j=1;j<=m;j++)map[i][j]=str[j]-'0';
}
for(i=1;i<=n;i++)
for(j=1;j<=m;j++)
spfa(i,j);
printf("%.6lf\n",sqrt(1.0*ans));//此处写成printf(".6lf\n",sqrt(1.0*ans));竟然半天没看出来
return 0;
}
24.bzoj 1821 //1821: [JSOI2010]Group 部落划分 Group //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1821
//1821: [JSOI2010]Group 部落划分 Group
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=1821
//最小生成树，去掉合并的边，第一个纳入部落间的最小边，即为答案
//因计算int没有溢出，故不采用long long
//Kruskal算法
//样例通过，提交Compile_Error，一看#include <cstring>漏了
//样例通过，提交AC。2019-6-10
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#define maxn 1010
using namespace std;
int n,k,cnt=0,f[maxn];
struct Point{
int x,y;
}p[maxn];
struct node{
int a,b,v;
}e[maxn*maxn/2];
int dis(Point a,Point b){
return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}
int cmp(node a,node b){
return a.v<b.v;
}
int getf(int u){
return f[u]==u?u:f[u]=getf(f[u]);
}
int merge(int u,int v){
int f1=getf(u),f2=getf(v);
if(f1!=f2){
f[f2]=f1;
return 1;
}
return 0;
}
int main(){
int i,j,count=0;
scanf("%d%d",&n,&k);
for(i=1;i<=n;i++)scanf("%d%d",&p[i].x,&p[i].y);
for(i=1;i<=n;i++)f[i]=i;
memset(e,0,sizeof(e));
for(i=1;i<=n;i++)
for(j=i+1;j<=n;j++)
cnt++,e[cnt].a=i,e[cnt].b=j,e[cnt].v=dis(p[i],p[j]);
sort(e+1,e+1+cnt,cmp);
for(i=1;i<=cnt;i++)
if(merge(e[i].a,e[i].b)){
count++;
if(count==((n-1)-(k-1)+1)){
printf("%.2lf\n",sqrt(1.0*e[i].v));
break;
}
}
return 0;
}
25.bzoj 2134 //2134: 单选错位 //在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2134
//2134: 单选错位
//在线测评地址https://www.lydsy.com/JudgeOnline/problem.php?id=2134
//第一目标，弄懂样例，弄明白，根据以下程序，弄明白了样例输入数据的作用
#include <stdio.h>
int a[10000010];
int main(){
int n,A,B,C;
scanf("%d%d%d%d%d",&n,&A,&B,&C,a+1);
for (int i=2;i<=n;i++) a[i] = ((long long)a[i-1] * A + B) % 100000001;
for (int i=1;i<=n;i++) a[i] = a[i] % C + 1;
for(int i=1;i<=n;i++)printf("a[%d]=%d\n",i,a[i]);
return 0;
}
//样例中a序列数据如下
a[1]=2
a[2]=3
a[3]=1
//样例弄懂了，题目怎么做，要学新的内容吗。2019-6-10 20:27
//此文思路相当棒https://www.cnblogs.com/five20/p/9383549.html摘抄如下
//此文思路写得不错https://www.luogu.org/problemnew/solution/P1297?page=2摘抄如下
//样例通过，提交AC。2019-6-10 22:43
#include <stdio.h>
int max(int a,int b){
return a>b?a:b;
}
int main(){
int n,A,B,C,a,b,i,x;
double ans=0;
scanf("%d%d%d%d%d",&n,&A,&B,&C,&a);
x=a;
for(i=2;i<=n;i++)b=((long long)a*A+B)%100000001,ans+=1.0/max(a%C+1,b%C+1),a=b;//此处写成b=b%C+1,ans+=1.0/max(a,b),a=b;
ans+=1.0/max(x%C+1,b%C+1);//此处写成ans+=1.0/max(x%C+1,b);
printf("%.3lf\n",ans);
return 0;
}
2019-8-4 18:00 AC该篇内容．
342
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
