官方解题报告:http://blog.sina.com.cn/s/blog_a19ad7a10102uyet.html
ZCC Loves Intersection
ZCC Loves COT
首先考虑一维下的版本。(数列,区间加,区间和)
显然可以使用线段树,但是线段树推广到高维度的难度较大。
针对本题先修改再询问的特点,我们可以在修改时只保存差分后的数组。
在处理询问前求前缀和还原原数列。之后再求一次前缀和就可以做到O(1)回答询问了。
现在试图把它推广到二维。(三角形,子三角形加,子三角形和)
如果朴素地对每一行应用一维的差分法,最后得到的标记会像下面一样:
其中+1标记和-1标记都是连续的一段。
标记也是一种值,所以可以对标记打标记!
于是我们将所有的+1标记按从右上到左下做差分,将所有的-1标记按从左上到右下做差分。
这样每个操作实际上只需要在两个数组中修改四个值,最后利用这两个数组还原原来的标记,利用标记还原原数阵。前缀和的计算也是类似,只不过第二维的前缀和要从两个方向各算一遍。
类似地,只需要4个数组就可以表示三维情况下的所有标记,经过三轮不同方向上的求前缀和就可以得到原四面体,前缀和的计算也需要4个数组。这样单操作/单询问要拆成8个,是O(1)的,复杂度是O(N^3+M+Q)。
ZCC Loves RPG
在程序的每一个位置,都存在若干条形如a[x]>=v或!(a[x]>=v)的限制。它们给每个a[i]确定了一个上界和一个下界。显然,我们应当让a[i]尽量小,因此我们可以令a[i]恰取到它的下界。
那么,考虑一对上下界数列Ui, Di,当前语句可以被执行到的条件就是:
随着我们分析程序的过程,一些位置的下界会发生变化,我们需要随时重新计算这个最大非零子段和,这就是说,我们需要对一个数列支持单点修改和查询最大非零子段和。这是一个线段树的经典应用,这里不再赘述。
下面考虑如何分析程序。
首先需要实现一个词法分析器:它接受字符串作为其输入,每次返回一个记号(常量,符号或标识符等)。这部分可以按照写读入优化的方法来写。
然后就到了对记号流进行分析的步骤,这一步的做法有很多,这里只讲一下std的做法。
首先将game(n, k)读入,获取n, k的值,同时建立线段树。
清空两个栈,这两个栈一个用于保存程序结构(称为栈P),一个用于保存Ui和Di的变化便于之后撤销(称为栈Q)。
读入左花括号,向栈P中压入一个B符号,B符号声明当前语句处于一对未闭合的花括号内,当栈P顶是B符号时,可以向这个花括号内追加语句。
随后重复以下过程直至栈P空。
检查栈P的栈顶:
若为B:从词法分析器取得下一个记号。
若为右花括号:弹出B。(闭合花括号)
否则:把这个记号“塞回去”,向栈P压入S符号。(追加语句)
若为S:弹出S,从词法分析器取得下一个记号。
若为左花括号:向栈P压入B符号。(新建块)
若为cg:读入整个cg语句,查询当前是否可达,更新答案。
若为if:读入if语句的第一句,读入x, v。向栈P先后压入I符号,T符号和S符号。向栈Q先后压入!(a[x]>=v)和a[x]>=v,在线段树上做a[x]>=v的修改。
若为T:弹出T,弹出栈Q栈顶元素,撤销它的修改。从词法分析器取得下一个记号:
若为else:应用栈Q栈顶元素(之前加入了但没有应用的!(a[x]>=v))的修改。向栈P压入S符号。
否则:把记号塞回去。弹出栈Q栈顶元素,弹出栈P栈顶元素(必定是I符号)。
若为I:弹出I,弹出栈Q栈顶元素,撤销它的修改。
B,S,T,I四个符号的意思分别是:块,语句,Then结束,整个If结束。
最后输出答案。
ZCC Loves Minecraft
所求的S即为当前点集的正交凸包。
参见en.wikipedia.org/wiki/Orthogonal_convex_hull的相关介绍。直观地看正交凸包是左上、右上、左下、右下四段折线构成的多边形,每条边都与坐标轴平行。求已知点集的正交凸包,可以先找出最上、最下、最左、最右的点。求右上方的一段折线,可以从最上面的点开始,每次找当前点右边的点中最上面的点,与当前点L形连接,作为新的当前点,不断重复直到选到最右边的点。这一过程实现非常简单,只要先按横坐标排序,然后扫描一遍并维护纵坐标最大值即可。其它三段的求法是类似的,而且可以通过旋转进行转化。那么对于动态加点的问题,可以用4棵平衡树(实现用STL即可)分别维护四段折线。插入点时尝试将其分别插入四段折线的相应位置,若在当前折线外部需要更新折线,并计算面积的增量。以右上折线为例,需要不断判断新点左边的点是否在新点下方,若是则删除。由于每个点最多被插入一次删除一次,总的时间复杂度是O(nlogn)。边界情况需要一定的特判。
ZCC Loves words
这道题我们可以先对单词建出AC自动机。然后在AC自动机上进行计数(dp)。F[i][j]表示当前长度第i位,当前在AC自动机的j号节点。我们考虑主动转移,如果j号节点在接受c这个字符后到了第k个节点,那么F[i+1][k]+=F[i][j]*Get。其中Get表示在第k个节点能乘上的数字。
我们来分析一下Get。,那么如果没有+i这个东西,我们就可以建出矩阵,然后直接快速幂+矩阵乘法。但是有了+i的话,每个矩阵都不同,但是事实上P=179*173*163,这其中的质数非常的小,然后就是乱搞的部分了,i Mod Pi就是O(Pi)的循环,于是我们建出Pi个矩阵,然后对于L/Pi的部分快速幂+矩阵乘法,对于L Mod Pi的部分暴力矩阵乘法即可。对于每一个质数得到的答案,利用中国剩余定理就可算出最后的答案。复杂度大概是O(Pi*tot^3 + log(L)*tot^3)。
如果有更优的做法欢迎指教。
ZCC Loves cards
ZCC Love ranklist
第一问:
当k=1的时候答案显然为n-1。
当k=2t(t∈Z*)时答案显然为0。把比赛分成两两一组,每组和为n+1即可。
当k=2t+1(t∈Z*)时先两两一组到只剩3个。n为偶数时答案显然不能为0,只能做到1。奇数时则可以做到。一种解决办法是:
n为奇数 |
1 |
4 |
7 |
n为偶数 |
1 |
4 |
6 |
2 |
5 |
5 |
2 |
5 |
4 |
||
3 |
6 |
3 |
3 |
6 |
2 |
||
4 |
7 |
1 |
4 |
1 |
5 |
||
5 |
1 |
6 |
5 |
2 |
3 |
||
6 |
2 |
4 |
6 |
3 |
1 |
||
7 |
3 |
2 |
第二问:
先观察进步指数的性质:由于NewRank=OldRank的时候进步指数=0所以只能出现一次所以不应该使用。当NewRank!=OldRank的时候两个进步指数相等当且仅当NewRank和OldRank都相等,两个进步指数为相反数当且仅当NewRank1=OldRank2,OldRank1=NewRank2。
总共能出现n*(n-1)+1种不同的进步指数,而且n>1,所以每种进步指数(除了0)都会出现。并且进步指数为相反数的一对只能出现在同一场考试中。所以n为奇数时显然无解。
当n为偶数时,注意到每场考试实际上是将考生两两组队n-1轮且不重复。于是使用循环赛构造算法即可。
循环赛构造算法:http://www.doc88.com/p-694165485213.html
ZCC Love march
我们用set维护当前的士兵位置,合并的时候将每个士兵暴力合并到该位置,即删掉原位置士兵并使目标位置size+=原位置size。每次移动的时候将原位置size--并在新位置新建一个size为1的士兵。这样最多只会有n+移动数 的士兵,每个士兵只会合并一次,所以时间复杂度为O(nlogn)。实现的时候为了记录每个点的坐标可以用并查集维护所在的块。
ZCC Love traffic lights
首先考虑如果是一般图的话怎么做。可以证明存在一个最优解在某一条边上卡着时间进去或者卡着时间出来。因为如果不这样的话可以将每个点的时间往后推直到卡着时间。所以只要对在每个点上卡时的情况进行判断。
枚举了某个点的到达时间以后那么接下来我们可以使用最短路算法求出到终点的最早时间和从起点出发的最晚时间。具体实现的时候可以对每个点开8个状态记录从哪个方向进入和是否闯过红灯。
ZCC loves Army
Solution:
可以发现上司和下属之间的关系可以构成一棵树,考虑三种操作。
交换下属:为了使交换下属的时候改变的边变少,可以用左儿子右兄弟的方式表示这棵树。那么交换时只需要改变至多4条边。可以用LCT维护。
传送信息:求从x传信息到y所需要的中介人的最少个数。令x为深度较大的一个点,这条路径就是x向上传到和y同一层,再水平传送。可以把编号为1的儿子的权值赋值为1,求x,y简单路径上的权值和。
发送指令:求x可以收到几个士兵发出的指令。表现在左儿子右兄弟的树上即为求该点到根之间的点个数。
Postscript:
为了方便,可以给每个点加一个编号为0的虚孩子,可以避免一些细节问题。
交换下属的时候找孩子可以对每个点开一个平衡树,也可以直接用LCT里的Splay。
ZCC loves Codefires
考察序列中相邻的两题i, j(i在前)。交换它们后,解出它们之前的题目所带来的时间对答案的贡献是不变的,它们对它们后面的题目的贡献也是不变的,其他题目之间对答案的贡献自然也是不变的。唯一的变化就是,原来的EiKj一项变成了EjKi一项。那么,为了使答案变优,需要满足的条件是EjKi≤EiKj。也即Ei/Ki≥Ej/Kj。
那么,最优解序列Ai一定满足,EAi/KAi是递增的。
排序一遍即可。
ZCC Loves Intersection http://acm.hdu.edu.cn/showproblem.php?pid=4873
公式见题解。我没推。
1 import java.math.BigInteger; 2 import java.util.Scanner; 3 4 5 class Main{ 6 public static void main(String arg[]){ 7 Scanner cin=new Scanner(System.in); 8 BigInteger n,fenzi,fenmu,tmp,gcd; 9 int d; 10 while(cin.hasNext()){ 11 n=cin.nextBigInteger(); 12 d=cin.nextInt(); 13 tmp=n.add(BigInteger.valueOf(4)); 14 fenzi=tmp.multiply(tmp); 15 fenzi=fenzi.multiply(BigInteger.valueOf(d)); 16 fenzi=fenzi.multiply(BigInteger.valueOf(d-1)); 17 fenmu=n.pow(d); 18 fenmu=fenmu.multiply(BigInteger.valueOf(18)); 19 gcd=fenmu.gcd(fenzi); 20 fenzi=fenzi.divide(gcd); 21 fenmu=fenmu.divide(gcd); 22 if(!fenzi.equals(fenmu)){ 23 System.out.println(fenzi+"/"+fenmu); 24 } 25 else { 26 System.out.println(1); 27 } 28 } 29 } 30 }
ZCC loves cards http://acm.hdu.edu.cn/showproblem.php?pid=4876
二进制暴力加一个剪枝。勉强过。
1 #include<cstdio> 2 #include<cstring> 3 #include<algorithm> 4 #define mt(a,b) memset(a,b,sizeof(a)) 5 using namespace std; 6 const int M=256; 7 int a[M]; 8 int b[M]; 9 bool vis[M]; 10 int sum[M]; 11 int main(){ 12 int n,k,L; 13 while(~scanf("%d%d%d",&n,&k,&L)){ 14 for(int i=0;i<n;i++){ 15 scanf("%d",&a[i]); 16 } 17 int all=1<<n; 18 int ans=0; 19 for(int i=0;i<all;i++){ 20 int tmp=i; 21 int num=0; 22 while(tmp){ 23 if(tmp&1) num++; 24 tmp>>=1; 25 } 26 if(num==k){ 27 int lb=0; 28 for(int j=0;j<n;j++){ 29 if((i>>j)&1){ 30 b[lb++]=a[j]; 31 } 32 } 33 int big=1<<k; 34 mt(vis,0); 35 for(int j=0;j<big;j++){ 36 int add=0; 37 for(int u=0;u<k;u++){ 38 if((j>>u)&1) add^=b[u]; 39 } 40 vis[add]=true; 41 } 42 int R=0; 43 for(int j=L;j<M;j++){ 44 if(!vis[j]) break; 45 R=j; 46 } 47 if(ans>=R) continue; //cut 48 sort(b,b+lb); 49 do{ 50 for(int j=0;j<k;j++){ 51 b[j+k]=b[j]; 52 } 53 sum[0]=b[0]; 54 int k2=k+k; 55 for(int j=1;j<k2;j++){ 56 sum[j]=sum[j-1]^b[j]; 57 } 58 mt(vis,0); 59 for(int u=0;u<k;u++){ 60 for(int v=u;v<k2&&v-u<k;v++){ 61 int add=sum[v]; 62 if(u-1>=0) add^=sum[u-1]; 63 vis[add]=true; 64 } 65 } 66 int R=0; 67 for(int j=L;j<M;j++){ 68 if(!vis[j]) break; 69 R=j; 70 } 71 ans=max(ans,R); 72 }while(next_permutation(b,b+lb)); 73 } 74 } 75 printf("%d ",ans); 76 } 77 return 0; 78 }
ZCC loves march http://acm.hdu.edu.cn/showproblem.php?pid=4879
数据要离散化,并查集,有点吊。
1 #include<cstdio> 2 #include<iostream> 3 #include<map> 4 #include<set> 5 using namespace std; 6 typedef __int64 LL; 7 typedef pair<LL,LL> pll; 8 const int mod=1000000007; 9 const int M=200010; 10 int f[M],sz[M],MXtot,MYtot,link[M]; 11 LL lans; 12 struct Input { 13 LL x,y; 14 } a[M]; 15 map<pll,LL>pos; 16 set<pll>PX[M],PY[M]; 17 set<pll>::iterator it; 18 map<LL,int>MX,MY; 19 int getroot(int x) { 20 if(x==f[x])return x; 21 return f[x]=getroot(f[x]); 22 } 23 void ins(int id) { 24 pll now=make_pair(a[id].x,a[id].y); 25 if(pos[now]) { 26 sz[pos[now]]++; 27 f[id]=pos[now]; 28 return; 29 } 30 pos[now]=id; 31 if(!MX[a[id].x]) MX[a[id].x]=++MXtot; 32 if(!MY[a[id].y]) MY[a[id].y]=++MYtot; 33 PX[MX[a[id].x]].insert(make_pair(a[id].y,id)); 34 PY[MY[a[id].y]].insert(make_pair(a[id].x,id)); 35 } 36 int main() { 37 int n,q,id; 38 LL m; 39 while(~scanf("%d%I64d",&n,&m)) { 40 for(int i=1; i<=n; i++) { 41 f[i]=i; 42 sz[i]=1; 43 link[i]=i; 44 } 45 for(int i=1; i<=n; i++) { 46 scanf("%I64d%I64d",&a[i].x,&a[i].y); 47 ins(i); 48 } 49 scanf("%d",&q); 50 while(q--) { 51 char op[4]; 52 scanf("%s",op); 53 if(op[0]=='Q') { 54 scanf("%d",&id); 55 id^=lans; 56 id=link[id]; 57 int fa=getroot(id); 58 lans=0; 59 LL X=a[fa].x; 60 LL Y=a[fa].y; 61 LL XX=MX[X]; 62 LL YY=MY[Y]; 63 LL pf,ps; 64 for(it=PX[XX].begin();it!=PX[XX].end();it++){ 65 pf=it->first; //y 66 ps=it->second; //id 67 if(ps!=fa) { 68 LL dy=(pf-Y)%mod; 69 dy=dy*dy%mod; 70 lans=(lans+dy*sz[ps]%mod)%mod; 71 pos.erase(make_pair(X,pf)); 72 PY[MY[pf]].erase(make_pair(X,ps)); 73 sz[fa]+=sz[ps]; 74 f[getroot(ps)]=fa; 75 } 76 } 77 PX[XX].clear(); 78 PX[XX].insert(make_pair(Y,fa)); 79 for(it=PY[YY].begin();it!=PY[YY].end();it++){ 80 pf=it->first; //x 81 ps=it->second; //id 82 if(ps!=fa) { 83 LL dx=(pf-X)%mod; 84 dx=dx*dx%mod; 85 lans=(lans+dx*sz[ps]%mod)%mod; 86 pos.erase(make_pair(pf,Y)); 87 PX[MX[pf]].erase(make_pair(Y,ps)); 88 sz[fa]+=sz[ps]; 89 f[getroot(ps)]=fa; 90 } 91 } 92 PY[YY].clear(); 93 PY[YY].insert(make_pair(X,fa)); 94 printf("%d ",lans); 95 } 96 else { 97 LL y; 98 scanf("%d%I64d",&id,&y); 99 id^=lans; 100 int p=link[id]; 101 sz[getroot(p)]--; 102 n++; 103 link[id]=n; 104 id=getroot(p); 105 a[n]=a[id]; 106 f[n]=n; 107 sz[n]=1; 108 pos.erase(make_pair(a[id].x,a[id].y)); 109 if(op[0]=='U') { 110 a[n].x-=y; 111 } 112 else if(op[0]=='D'){ 113 a[n].x+=y; 114 } 115 else if(op[0]=='L'){ 116 a[n].y-=y; 117 } 118 else if(op[0]=='R'){ 119 a[n].y+=y; 120 } 121 ins(n); 122 } 123 } 124 } 125 return 0; 126 }
ZCC Loves Codefireshttp://acm.hdu.edu.cn/showproblem.php?pid=4882
推出一个优先关系,排序。
1 #include<cstdio> 2 #include<algorithm> 3 using namespace std; 4 typedef __int64 LL; 5 const int M=100010; 6 struct G{ 7 int t,k; 8 friend bool operator <(G a,G b){ 9 return a.t*b.k<a.k*b.t; 10 } 11 }g[M]; 12 int main(){ 13 int n; 14 while(~scanf("%d",&n)){ 15 for(int i=0;i<n;i++){ 16 scanf("%d",&g[i].t); 17 } 18 for(int i=0;i<n;i++){ 19 scanf("%d",&g[i].k); 20 } 21 sort(g,g+n); 22 LL ans=0,tmp=0; 23 for(int i=0;i<n;i++){ 24 ans+=(tmp+g[i].t)*g[i].k; 25 tmp+=g[i].t; 26 } 27 printf("%I64d ",ans); 28 } 29 return 0; 30 }
end