CQOI2016 day1 总结

最新推荐文章于 2020-04-03 11:05:04 发布
最新推荐文章于 2020-04-03 11:05:04 发布
阅读量837 收藏
点赞数
文章标签： 最小割 数位dp 乱搞 分治 kd树
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/lcrtest/article/details/51136020
版权

怎么发一次格式抽一次

T1
打表发现答案小于n
考虑求出S,T的最小割，并且提取出在S这一边和在T这一边的点集
考虑S'属于S这一边的点集，T'属于T这一边的点集，那么S'和T'的最小割可以用其他两个点代替
借助bzoj2229黄学长的一句话：
注意这样一个事实：如果(X,Y)是某个s1-t1最小割，(Z,W)是某个s2-t2最小割，那么X∩Z、X∩W、Y∩Z、Y∩W这四项不可能均非空。也就是说，最小割不可能相互跨立。
这个蕴含了，最多一共有N-1个不同的s-t最小割。只需把这些割找出来即可。
寻找的方法：首先，在V中任意找两个点a,b，求最大流，把V划分为割X-Y，之后对X、Y分别递归地进行划分。这样就能得到N-1个割了。
那么我们可以分治求最小割，这样求得的最小割一定能够代表所有的最小割


T2
第一种：考场上YY出来的，因为太懒所以没写
首先求凸包，求出凸包上每一个点所能到达的最远点，push到堆中
对全局维护一个堆，每次弹出最大值，然后讲两个点的次远点在kd树中找出来
找出来的复杂度是sqrt(n)的，外加二份答案，共k次，总复杂度k*sqrt(n)*log（没有代码）
第二种
还是求出凸包和每一个点所能到达的最远点，每次拿出堆中的时候暴力扫n个点，强行插入次远的点
所以下面这个是nklog的。也能过


T3
傻逼数位dp
dp[当前位数][是否有8][是否有4][上一位][连续几位][已经有没有3位连续的][是否达到上限]
dfs一遍记忆化转移就好 
顺便说一句，函数记得return啊！！！关键是小数据还是对的调崩我了


代码

T1

//Copyright(c)2016 liuchenrui
#include<bits/stdc++.h>
#define short unsigned short
#define inf 1000000000
using namespace std;
template<typename T>
inline void splay(T &v){
	v=0;char c=0;T p=1;
	while(c<'0' || c>'9'){if(c=='-')p=-1;c=getchar();}
	while(c>='0' && c<='9'){v=(v<<3)+(v<<1)+c-'0';c=getchar();}
	v*=p;
}
struct Edge{
	short to,next;int flow;
}edge[100010];
short first[851],size;
short deep[851],dl[851];
void addedge(short x,short y,int z){
	size++;
	edge[size].to=y;
	edge[size].next=first[x];
	first[x]=size;
	edge[size].flow=z;
}
void add(short x,short y,int z){
	addedge(x,y,z),addedge(y,x,0);
}
short aim;
int dfs(short now,int flow){
	if(now==aim)return flow;
	int F=0;
	for(short u=first[now];u&&flow;u=edge[u].next){
		if(deep[edge[u].to]==deep[now]+1&&edge[u].flow){
			int tmp=dfs(edge[u].to,min(flow,edge[u].flow));
			F+=tmp;edge[u].flow-=tmp;edge[u^1].flow+=tmp;flow-=tmp;
		}
	}
	if(!F)deep[now]=-5;
	return F;
}
bool bfs(short S,short T){
	memset(deep,0,sizeof(deep));
	dl[1]=S;deep[S]=1;
	short head=0,tail=1;
	while(head!=tail){
		head++;
		for(int u=first[dl[head]];u;u=edge[u].next){
			if(!deep[edge[u].to]&&edge[u].flow){
				dl[++tail]=edge[u].to;
				deep[edge[u].to]=deep[dl[head]]+1;
			}
		}
	}
	return deep[T];
}
int maxflow(short S,short T){
	int ret=0;aim=T;
	while(bfs(S,T)){
		ret+=dfs(S,inf);
	}
	return ret;
}
short fr[8600],to[8600];
int fl[8600];
int ans[1000010];
short cnt,n,m,p[851];
void rebuild(){
	memset(first,0,sizeof first);
	size=1;
	for(int k=1;k<=m;k++){
		add(fr[k],to[k],fl[k]);
		add(to[k],fr[k],fl[k]);
	}
}
bool vis[860];
void dfs2(short now){
	vis[now]=1;
	for(int u=first[now];u;u=edge[u].next){
		if(edge[u].flow&&!vis[edge[u].to]){
			dfs2(edge[u].to);
		}
	}
}
void cdq_calc(short l,short r){
	if(l>=r)return;
	rebuild();
	ans[++cnt]=maxflow(p[l],p[r]);
	memset(vis,0,sizeof vis);
	dfs2(p[l]);int t=l-1;
	for(short i=l;i<=r;i++){
		if(vis[p[i]]){
			swap(p[i],p[++t]);
		}
	}
	cdq_calc(l,t),cdq_calc(t+1,r);
}
int main(){
	freopen("cuts.in","r",stdin);
	freopen("cuts.out","w",stdout);
	splay(n),splay(m);
	for(short i=1;i<=m;i++){
		splay(fr[i]),splay(to[i]),splay(fl[i]);
	}
	for(short i=1;i<=n;i++){
		p[i]=i;
	}
	cdq_calc(1,n);
	sort(ans+1,ans+cnt+1);
	cnt=unique(ans+1,ans+cnt+1)-ans-1;
	printf("%d\n",cnt);
	//cerr<<clock()<<endl;
}
T2是czh的代码，自己的WA没调出来 

#include <iostream>
#include <cstring>
#include <cstdlib>
#include <string>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <ctime>
#include <map>
#define inf (((unsigned long long)1 << (unsigned long long)63) - (unsigned long long)1)
#define Int long long
using namespace std;
struct point {Int x;Int y;int id;
};point p[100010],q[100010],stack[200010];
Int val[800010],minx[850],maxx[850];
int E,belong[100010],cnt,n,k;
int head,tail,father[100010],t,num[100010];
int tot[850];
map < pair<int,int>,int> mp;
Int calc(point x,point y,point z) {
	point A,B;
	A.x = x.x - z.x;
	A.y = x.y - z.y;
	B.x = y.x - z.x;
	B.y = y.y - z.y;
	return A.x * B.y - A.y * B.x;
}
Int dis(point x,point y) {
	Int ret = (x.x - y.x) * (x.x - y.x) + (x.y - y.y) * (x.y - y.y);
	return ret;
}
bool comp(const point &x,const point &y) {
	Int v = calc(x,y,p[1]);
	if(v == 0) return dis(x,p[1]) < dis(y,p[1]);
	return v > 0;
}
void build(int Now,int l,int r) {
	val[Now] = -1;
	if(l == r) return;
	int Mid = (l + r) >> 1;
	build(Now << 1,l,Mid);
	build(Now << 1 | 1,Mid + 1,r);
}
void change(int Now,int l,int r,int x,Int y) {
	if(l == r) {val[Now] = y;return ;}
	int Mid = (l + r) >> 1;
	if(x <= Mid) change(Now << 1,l,Mid,x,y);
	else change(Now << 1 | 1,Mid + 1,r,x,y);
	if(val[Now << 1] > val[Now << 1 | 1])
		val[Now] = val[Now << 1];
	else val[Now] = val[Now << 1 | 1];
}
int Ask(int Now,int l,int r) {
	if(l == r) return r;
	int Mid = (l + r) >> 1;
	if(val[Now << 1] == val[Now])
		return Ask(Now << 1,l,Mid);
	else return Ask(Now << 1 | 1,Mid + 1,r);
}
int find(int Now,int l,int r) {
	if(l == r) return r;
	int Mid = (l + r) >> 1;
	if(minx[Now] == minx[Now << 1])
		return find(Now << 1,l,Mid);
	else return find(Now << 1 | 1,Mid + 1,r);
}
void insert(int Now,int l,int r,int x,Int y) {
	if(l == r) {minx[Now] = maxx[Now] = y;tot[Now] = 1;return ;}
	int Mid = (l + r) >> 1;
	if(x <= Mid) insert(Now << 1,l,Mid,x,y);
	else insert(Now << 1 | 1,Mid + 1,r,x,y);
	tot[Now] = tot[Now << 1] + tot[Now << 1 | 1];
	minx[Now] = min(minx[Now << 1],minx[Now << 1 | 1]);
	maxx[Now] = max(maxx[Now << 1],maxx[Now << 1 | 1]);
}
int bf(int x) {
	int Z = num[x] + 1;
	for(int i = 1;i <= 810;i ++)
		minx[i] = inf,maxx[i] = -1,tot[i] = 0;
	for(int i = 1;i <= n;i ++)
	{
		Int G = dis(q[i],q[x]);
		if(tot[1] < Z)
			insert(1,1,100,tot[1] + 1,G);
		else if(minx[1] < G) 
			insert(1,1,100,find(1,1,100),G);
	}
	for(int i = 1;i <= n;i ++)
		if(dis(q[i],q[x]) == minx[1] && mp[make_pair(x,i)] == 0)
			return i;			
}
int getint() {
	char c = 'd';
	int ret = 0;
	while(c < '0' || c > '9') c = getchar();
	while(c >= '0' && c <= '9') ret = ret * 10 + c - '0',c = getchar();
	return ret;
}
int main() {
	freopen("farthest.in","r",stdin);
	freopen("farthest.out","w",stdout);
	n = getint();
	k = getint();
	for(int i = 1;i <= n;i ++)
		p[i].x = getint(),p[i].y = getint(),q[i].x = p[i].x,q[i].y = p[i].y,p[i].id = i;
	t = 1;
	for(int i = 2;i <= n;i ++)
		if(p[i].y < p[t].y || (p[i].y == p[t].y && p[i].x < p[t].x))
			t = i;
	swap(p[1],p[t]);
	sort(p + 2,p + n + 1,comp);
	head = 1;tail = 0;
	for(int i = 1;i <= n;i ++) 
	{
		while(head < tail && calc(p[i],stack[tail],stack[tail - 1]) > 0)
			tail --;
		stack[ ++ tail] = p[i];
	}
	for(int i = 1;i <= tail;i ++)
		stack[i + tail] = stack[i];
	int M = 1;
	build(1,1,n);
	for(int i = 1;i <= tail;i ++)
	{
		while(M < 2 * tail && dis(stack[M + 1],stack[i]) > dis(stack[M],stack[i]))
			M = M + 1;
		change(1,1,n,stack[i].id,dis(stack[i],stack[M]));father[stack[i].id] = stack[M].id;
	}
	for(int i = 1;i < k;i ++) 
	{
		int t = Ask(1,1,n),w = father[t];
		mp[make_pair(w,t)] = mp[make_pair(t,w)] = 1;
		num[t] ++;
		num[w] ++;
		father[t] = bf(t);
		father[w] = bf(w);
		change(1,1,n,t,dis(q[t],q[father[t]]));
		change(1,1,n,w,dis(q[w],q[father[w]]));
	}
	cout<<val[1];
}

T3 

//Copyright(c)2016 liuchenrui
#include<bits/stdc++.h>
#define LL long long
using namespace std;
inline void splay(int &v){
	v=0;char c=0;int p=1;
	while(c<'0' || c>'9'){if(c=='-')p=-1;c=getchar();}
	while(c>='0' && c<='9'){v=(v<<3)+(v<<1)+c-'0';c=getchar();}
	v*=p;
}
int s[15];
LL dp[13][2][2][10][12][2][10];//当前位数，是否有8,是否有4,上一位，连续几位，已经有没有3位连续的，是否达到上限 
LL dfs(int a,int b,int c,int d,int e,int f,int g){
	if(a>11){
		if(b&&c)return 0;
		if(!f)return 0;
		return 1;
	}
	if(dp[a][b][c][d][e][f][g]!=-1)return dp[a][b][c][d][e][f][g];
	LL ret=0;int t;
	if(g==0){
		for(int i=0;i<=9;i++){
			if(i==d)t=e+1;else t=1;
			ret+=dfs(a+1,(b|(i==8)),(c|(i==4)),i,t,(t>=3)|(f),0);
		}
	}
	else{
		for(int i=0;i<s[a];i++){
			if(i==d)t=e+1;else t=1;
			ret+=dfs(a+1,(b|(i==8)),(c|(i==4)),i,t,(t>=3)|(f),0);
		}
		if(s[a]==d)t=e+1;else t=1;
		ret+=dfs(a+1,(b|(s[a]==8)),(c|(s[a]==4)),s[a],t,(t>=3)|(f),1);
	}
	return dp[a][b][c][d][e][f][g]=ret;
}
LL calc(LL a){
	if(a<10000000000LL)return 0;
	for(int i=11;i>=1;i--){
		s[i]=a%10;a/=10;
	}
	memset(dp,-1,sizeof dp);
	LL ret=0;
	for(int i=1;i<s[1];i++){
		ret+=dfs(2,i==8,i==4,i,1,0,0);
	}
	ret+=dfs(2,s[1]==8,s[1]==4,s[1],1,0,1);
	return ret;//nmb 记得return啊卧槽！！！！！！！ 
}
int main(){
	freopen("number.in","r",stdin);
	freopen("number.out","w",stdout);
	LL x,y;
	cin>>x>>y;
	LL a=calc(x-1);
	LL b=calc(y);
	cout<<b-a<<endl;
}



Mstdream
关注
[CQOI2016]不同的最小割
少主大人殿下的博客
12-13 469
题目 【问题描述】 　　学过图论的同学都知道最小割的概念：对于一个图，某个对图中结点的划分将图中所有结点分成两个部分，如果结点s,t不在同一个部分中，则称这个划分是关于s,t的割。对于带权图来说，将所有顶点处在不同部分的边的权值相加所得到的值定义为这个割的容量，而s,t的最小割指的是在关于s,t的割中容量最小的割。 　　而对冲刺NOI竞赛的选手而言，求带权图中两点的最小割已经不是什么难...
CQOI2016省选游记
liuyunhui246的补给站
04-12 1041
4.9～4.10。重庆邮电大学。每天5小时，3道题。今年完成更新的CQTSC。     第二次省选。依稀记得第一次省选，重庆最后一次5T5H心态赛，每道题都做不来的绝望和无助。网络流是啥？可持久化数据结构？莫比乌斯反演？连通性状压？二项式定理？统统SMG……然后没有等成绩直接回来，记得那天在讲 置换群 ，现在那几道题还没做完……
CQOI2016 总结
weixin_30755709的博客
04-10 113
2016.4.9 day 1 第一题：按照题目描述暴力，N^2枚举点对再跑网络流，最后50.据说正解是分治最小割，还是一道分治最小割的裸题，然而不知道 =_= 心塞 第二题：看到题目就被吓到了，最笨的N^2*log(N)的暴力肯定过不了，但是看到k<=100的时候心情又好了，觉得肯定有什么奇特的东西可以乱搞搞。最后还是暴力的，不过优化了一些，得了50分。用双链表维护...
cqoi2016
某蒟蒻的博客
05-23 477
cqoi2016不同的最小割 做法同zjoi2011最小割，题解略。#include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<string> #include<cmath> #include<algorithm> #define inf 1050000000 #define N 855 #de
CQOI2016
jzhang1的专栏
04-12 792
T1:     大意：一个850个点，8500条边的图，问所有点对之间的最小割有多少种不同的权值。     分治最小割。(然而做这道题之前完全不知道是什么..%YY出来的lcr)。     一开始所有点在一个集合中，随便找两个点求一次最小割，然后会把点集分成两半，继续直到所有点都变成一个点为止 (反正我们这儿唯一A的人就是这么说的...但是奇慢无比，应该是姿势不太对...)。     至
[CQOI2016]游记
liuguangzhe
04-11 991
CQOI2016
ceph day 2016
09-11
CEPH DAYS 2016 APAC ROADSHOW @BEIJING ceph 对象存储论坛
hibernate框架开发2016版视频 day01-day04
07-24
hibernate框架开发2016版视频 day01 hibernate框架开发2016版视频 day02 hibernate框架开发2016版视频 day03 hibernate框架开发2016版视频 day04
CQOI2016游记
鲜衣怒马少年郎
04-12 751
来一篇CQYZ视角的嘤嘤 分析什么的留到最后 两周之前开始停的课，成功水掉一次月考加一次运动会加晚自习无数，之前都是晚自习过来水一水，然而并没有水出什么东西 HB权衡利弊之后果断放弃了计算几何（FLAG），带领大家专攻高（级）数（据结构）和网络流建模（FLAG*2），大家平时在各班班主任的（唠叨）谆谆教诲下过着又不专心准备竞赛又不专心准备高考的生活。 大家一致认为自己复赛发挥不错（FLAG
CQOI2016滚粗记
蒟蒻的奋斗史
04-11 2153
day0 上去试机发现鼠标是坏的，旁边是巴蜀的大牛和教师机… 随手敲了一个树剖和Pollard_rho，顺路用批处理对拍了一下，然后想去交一发意外的发现竟然木有网….晚上复习了一下板子…然后写得有点兴奋了就睡不着辣…day1 上歌乐山的时候还在下雨… 然后我7点就悠悠的到了门口…然后就在外面晃荡了40多分钟… 8点过就进了考场，发现解压密码一点都不好玩…看到T1的时候心里一阵窃喜，嘿嘿我做
[cqoi2016]伪光滑数 解题报告
TA
04-18 1272
这题有点意思。。 考虑对于i个质因子，最大的质因子至多为j能生成的数。我们需要每次在其中取最大值，显然它可以用可持久化左偏树来维护。有leftist(i,j)=leftist(i−1,j)∗j[ji≤n]+leftist(i,j−1)leftist(i,j)=leftist(i-1,j)*j[j^i\le n]+leftist(i,j-1) 然后我们再用一个堆来维护所有可持久化左偏树的根的最小值
4520: [Cqoi2016]K远点对
CRZbulabula的博客
09-26 667
4520: [Cqoi2016]K远点对 Time Limit: 30 Sec  Memory Limit: 512 MB Submit: 594  Solved: 314 [Submit][Status][Discuss] Description 已知平面内 N 个点的坐标，求欧氏距离下的第 K 远点对。 Input 输入文件第一行为用空格隔开的两个整数 N， K。
[CQOI2016]不同的最小割 分治最小割
Shelter
03-23 407
[CQOI2016] 求一个无向图中所有点对之间的本质不同的最小割的个数 一张图最多有n-1种本质不同的最小割，考虑如何快速地找到所有不同的最小割 只要让每次划分的S，T集合存在区别，就对应了不同的最小割 我们考虑分治来解决这问题，首先对于我们当前分治的点集，我们任选 两个点 s,t，然后我们计算这两个点之间的最小割，这个时候整张图会 被我们划分为 S,T 两个点集 这个时候当前对于当前...
【CQOI2016】密钥破解
Facico的博客
04-20 1857
DescriptionSolution要求什么题目看起来挺复杂的，但是理一理思路，只有四个要求：1、能快速提取一个大整数的因数；2、会求逆元；3、会打快速幂；4、会打快速乘。 首先后面三个都是数论十分基础的操作。 但是第一个是怎么回事？怎么分解？有一个算法叫做Pollard’s Rho算法 想说，这个算法真的非常的神奇。 一个random算法，复杂度O（玄学），实际上普遍概率O(n14)O
[CQOI2016]伪光滑数
C202044zxy的博客
04-03 219
一、题目 点此看题 二、解法 首先处理出128128128以内的质数，kkk的范围比较小，可以一个一个找。 我们维护一个四元组(val,mn,lst,sec)(val,mn,lst,sec)(val,mn,lst,sec)，分别表示当前数的值，最小的质因子，最小的质因子次数，次小的质因子次数。一开始我们把每个质数结尾的最大数放进堆中，容易发现此时质因数都相同，我们需要保证结尾的质因数一定是当前钦定...
NOI2016 DAY2：区间选择最小花费问题
NOI2016（中国国家青少年信息学奥林匹克竞赛）的第二天竞赛中，有一个传统型题目叫做“区间”（Interval），其主要涉及了数轴上的区间选择问题。给定的数轴上有 n 个闭区间 [l1, r1], [l2, r2], ..., [ln, rn]，参赛...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值