码迷，mamicode.com
首页 > 其他好文 > 详细

Codeforces Round #431 Div. 1

时间：2019-03-09 23:38:19      阅读：171      评论：0      收藏：0      [点我收藏+]

标签：erase   set   for   data   相同   复杂   sort   a.out   ace   

　　A：显然每种字符的代价互不相关，dp并打表可得合并i个字符的最小代价是i*(i-1)/2。然后直接贪心分配每个字符即可。因为每次分配都将剩余代价降到了根号级别所以字符数量是足够的。

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
#define N 100010
char getc(){char c=getchar();while ((c<‘A‘||c>‘Z‘)&&(c<‘a‘||c>‘z‘)&&(c<‘0‘||c>‘9‘)) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
int read()
{
	int x=0,f=1;char c=getchar();
	while (c<‘0‘||c>‘9‘) {if (c==‘-‘) f=-1;c=getchar();}
	while (c>=‘0‘&&c<=‘9‘) x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x*f;
}
int n;
signed main()
{
#ifndef ONLINE_JUDGE
	freopen("a.in","r",stdin);
	freopen("a.out","w",stdout);
#endif
	n=read()*2;if (n==0) {cout<<‘a‘;return 0;}
	for (int i=0;i<26;i++)
	{
		for (int j=n;j>=2;j--)
		if (1ll*j*(j-1)<=n)
		{
			for (int k=1;k<=j;k++) putchar(‘a‘+i);
			n-=1ll*j*(j-1);
			break;
		}
		if (n==0) break;
	}
	return 0;
	//NOTICE LONG LONG!!!!!
}

　　B：先找出每个会撞在一起的集合，然后可以发现相当于是在网格图中向某方向走遇到一个格点就拐个弯，讨论一下即可。

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
#define N 100010
char getc(){char c=getchar();while ((c<‘A‘||c>‘Z‘)&&(c<‘a‘||c>‘z‘)&&(c<‘0‘||c>‘9‘)) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
int read()
{
	int x=0,f=1;char c=getchar();
	while (c<‘0‘||c>‘9‘) {if (c==‘-‘) f=-1;c=getchar();}
	while (c>=‘0‘&&c<=‘9‘) x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x*f;
}
int n,w,h,ansx[N],ansy[N];
struct data
{
	int op,x,t,i;
	bool operator <(const data&a) const
	{
		return x-t<a.x-a.t;
	}
}a[N];
bool cmp(const data&a,const data&b)
{
	return a.op<b.op||a.op==b.op&&a.x<b.x;
}
signed main()
{
#ifndef ONLINE_JUDGE
	freopen("a.in","r",stdin);
	freopen("a.out","w",stdout);
#endif
	n=read(),w=read(),h=read();
	for (int i=1;i<=n;i++) a[i].op=read(),a[i].x=read(),a[i].t=read(),a[i].i=i;
	sort(a+1,a+n+1);
	for (int i=1;i<=n;i++)
	{
		int t=i;
		while (t<n&&a[t+1].x-a[t+1].t==a[i].x-a[i].t) t++;
		sort(a+i,a+t+1,cmp);
		int x=0,y=0;
		for (int j=i;j<=t;j++) if (a[j].op==1) x++;else y++;
		for (int j=i;j<i+x;j++)
		{
			int W=i+x-j-1,H=y,X=min(W,H);
			if (W>=H) ansx[a[j].i]=a[j+X].x,ansy[a[j].i]=h;
			else ansx[a[j].i]=w,ansy[a[j].i]=a[i+x+X].x;
		}
		for (int j=i+x;j<=t;j++)
		{
			int W=t-j,H=x,X=min(W,H);
			if (W>=H) ansx[a[j].i]=w,ansy[a[j].i]=a[j+X].x;
			else ansx[a[j].i]=a[i+X].x,ansy[a[j].i]=h;
		}
		i=t;
	}
	for (int i=1;i<=n;i++) printf("%d %d\n",ansx[i],ansy[i]);
	return 0;
	//NOTICE LONG LONG!!!!!
}

　　C：显然相当于求区间相邻两相同数的位置差之和。于是可以看成一个二维查询，即r前缀中所有前驱>=l的位置的权值和。可以树状数组套treap动态维护。

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<set> 
using namespace std;
#define ll long long
#define N 100010
#define lson tree[k].ch[0]
#define rson tree[k].ch[1]
char getc(){char c=getchar();while ((c<‘A‘||c>‘Z‘)&&(c<‘a‘||c>‘z‘)&&(c<‘0‘||c>‘9‘)) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
int read()
{
	int x=0,f=1;char c=getchar();
	while (c<‘0‘||c>‘9‘) {if (c==‘-‘) f=-1;c=getchar();}
	while (c>=‘0‘&&c<=‘9‘) x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x*f;
}
int n,m,a[N],root[N],cnt;
set<int> pos[N];
struct data{int ch[2],p,x,v;ll s;
}tree[N<<5];
int pre(int p,int x)
{
	auto it=pos[x].find(p);it--;return *it;
}
void up(int k){tree[k].s=tree[lson].s+tree[rson].s+tree[k].v;}
void move(int &k,int p)
{
	int t=tree[k].ch[p];
	tree[k].ch[p]=tree[t].ch[!p],tree[t].ch[!p]=k,up(k),up(t),k=t;
}
void ins(int &k,int x,int v)
{
	if (k==0) {k=++cnt;tree[k].p=rand();tree[k].x=x;tree[k].s=tree[k].v=v;return;}
	tree[k].s+=v;
	if (tree[k].x<x) {ins(rson,x,v);if (tree[rson].p>tree[k].p) move(k,1);}
	else {ins(lson,x,v);if (tree[lson].p>tree[k].p) move(k,0);}
}
void del(int &k,int x,int v)
{
	if (tree[k].x==x) tree[k].v-=v;
	else if (tree[k].x<x) del(rson,x,v);
	else del(lson,x,v);
	up(k);
}
ll query(int k,int x)
{
	if (k==0) return 0;
	if (tree[k].x<x) return query(rson,x);
	else return tree[rson].s+tree[k].v+query(lson,x);
}
void Insert(int k,int x){int y=k-x;while (k<=n) ins(root[k],x,y),k+=k&-k;}
void Delete(int k,int x){int y=k-x;while (k<=n) del(root[k],x,y),k+=k&-k;}
ll Query(int k,int x){ll s=0;while (k) s+=query(root[k],x),k^=k&-k;return s;}
signed main()
{
#ifndef ONLINE_JUDGE
	freopen("a.in","r",stdin);
	freopen("a.out","w",stdout);
#endif
	srand(20020509);
	n=read(),m=read();
	for (int i=1;i<=n;i++) a[i]=read();
	for (int i=1;i<=n;i++) pos[i].insert(0);
	for (int i=1;i<=n;i++)
	{
		pos[a[i]].insert(i);
		Insert(i,pre(i,a[i]));
	}
	for (int i=1;i<=m;i++)
	{
		int op=read();
		if (op==1)
		{
			int p=read(),x=read(),y=pre(p,a[p]);
			if (a[p]==x) continue;
			auto it=pos[a[p]].find(p);it++;
			if (it!=pos[a[p]].end())
			{
				Delete(*it,p);
				Insert(*it,y);
			}
			pos[a[p]].erase(p);
			a[p]=x;
			pos[x].insert(p);
			int z=pre(p,x);
			it=pos[x].find(p);it++;
			if (it!=pos[x].end())
			{
				Delete(*it,z);
				Insert(*it,p);
			}
			Delete(p,y);Insert(p,z);
		}
		if (op==2)
		{
			int l=read(),r=read();
			printf("%I64d\n",Query(r,l));
		}
	}
	return 0;
	//NOTICE LONG LONG!!!!!
}

　　D：技术图片

　　搬上sol里的图，然后就变得很简单了。考虑初始边所衍生出的每个部分，这样就变成了一个子问题。大体类似于本质不同的有根树的计数，只是稍微复杂一点。

　　设f[i][j]为操作i次后最小割为j的方案数。考虑先求出一个部分的方案数。即设g[i][j]为操作i次最小割为j（不考虑初始边）且钦定初始边只被操作一次的方案数，转移显然，这一部分是O(n4)的。

　　然后考虑通过g推回f，由于相同的部分之间是无序的，类似上面的问题做一些操作即可。即设h[i][j][x][y]为操作i次后最小割为j且当前只考虑到(x,y)及其之前（如x‘<x||x‘==x&&y‘<y）的部分的方案数。转移时枚举(x,y)选择多少个，进行一个插板即可完成无序统计。复杂度是极小常数O(n5logn)。

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
#define P 1000000007
#define N 52
char getc(){char c=getchar();while ((c<‘A‘||c>‘Z‘)&&(c<‘a‘||c>‘z‘)&&(c<‘0‘||c>‘9‘)) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
int read()
{
	int x=0,f=1;char c=getchar();
	while (c<‘0‘||c>‘9‘) {if (c==‘-‘) f=-1;c=getchar();}
	while (c>=‘0‘&&c<=‘9‘) x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x*f;
}
int n,m,f[N][N],g[N][N],h[N][N][N][N],inv[N],ans;
void inc(int &x,int y){x+=y;if (x>=P) x-=P;}
int C(int n,int m)
{
	if (m>n) return 0;
	int s=1;
	for (int i=n-m+1;i<=n;i++) s=1ll*s*i%P;
	return 1ll*s*inv[m]%P;
}
int F(int n,int m){return C(n+m-1,m);}
signed main()
{
#ifndef ONLINE_JUDGE
	freopen("a.in","r",stdin);
	freopen("a.out","w",stdout);
#endif
	n=read(),m=read();int tmp=m;m=n+1;
	inv[0]=inv[1]=1;for (int i=2;i<=51;i++) inv[i]=P-1ll*(P/i)*inv[P%i]%P;
	for (int i=2;i<=50;i++) inv[i]=1ll*inv[i]*inv[i-1]%P;
	f[0][1]=1;g[0][0]=1;
	for (int x=0;x<=n;x++)
		for (int y=0;y<=m;y++)
		h[0][1][x][y]=1;
	for (int i=1;i<=n;i++)
	{
		for (int j=1;j<=m;j++) 
		{
			for (int x=0;x<i;x++)
			{
				for (int y=j+1;y<=m;y++)
				inc(g[i][j],1ll*f[x][y]*f[i-x-1][j]%P),
				inc(g[i][j],1ll*f[x][j]*f[i-x-1][y]%P);
				inc(g[i][j],1ll*f[x][j]*f[i-x-1][j]%P);
			}
		}
		for (int j=1;j<=m;j++)
		{
			for (int x=1;x<=n;x++)
				for (int y=1;y<=m;y++)
				{
					int ux=x,uy=y-1;if (uy==0) ux--,uy=m;
					h[i][j][x][y]=h[i][j][ux][uy];
					for (int k=1;k*x<=i&&k*y<j;k++)
					inc(h[i][j][x][y],1ll*h[i-x*k][j-y*k][ux][uy]*F(g[x][y],k)%P);
				}
			f[i][j]=h[i][j][n][m];
		}
	}
	m=tmp;
	cout<<f[n][m];
	return 0;
	//NOTICE LONG LONG!!!!!
}

　　

Codeforces Round #431 Div. 1

标签：erase   set   for   data   相同   复杂   sort   a.out   ace   

原文地址：https://www.cnblogs.com/Gloid/p/10503568.html
踩
(0)
赞
(0)
   
举报
评论 一句话评论（0）
登录后才能评论！
分享档案
更多>
2021年07月29日 (22)
2021年07月28日 (40)
2021年07月27日 (32)
2021年07月26日 (79)
2021年07月23日 (29)
2021年07月22日 (30)
2021年07月21日 (42)
2021年07月20日 (16)
2021年07月19日 (90)
2021年07月16日 (35)
周排行
mamicode.com排行更多图片
更多
友情链接
兰亭集智   国之画   百度统计   站长统计   阿里云   chrome插件   新版天听网
关于我们 - 联系我们 - 留言反馈
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码！