北京Day 12

今日得分：100+0+60+0

T1

题目大意：把一个正整数n分解k次，每次得到的收益为分解得到的两个数的值的乘积，求最大收益。n,k<=1e9

正解：打表找规律

对于一种分解后的序列，任意两个数的乘积， 都会被计算到最终的收益中。于是我们可以得到这样一个式子：

F n = x ∗ (x − 1) /2 设最终的收益为 P 则： P = F x −Σ F(li) 最小化后面的部分，可以让 P 更大，因此尽量平分就好。复杂度 O(1)。

T1AC代码（找规律）

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<vector>
#include<cmath>
#include<bitset>
using namespace std;
inline int re_ad()
{
	int x=0,f=1;char ch=getchar();
	while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9')x=x*10+ch-48,ch=getchar();
	return x*f;
}
inline long long ma(long long x,long long y){return x>y?x:y;}
long long n,k;
unsigned long long ans=0;
int main()
{
	register int i;
	register long long x,y,z;;
	n=re_ad();k=re_ad();
	x=k+1;y=n/x;
	z=x*(x-1)>>1;
	ans=z*y;x=x*k;
	--y;
	ans+=(1+y)*y/2*x;
	++y;
	x=y*k;z=x+(n%(k+1))-1;
	ans+=(x+z)*(z-x+1)/2;
	cout<<ans<<endl;
}

T2

题目大意：给 n 组操作，每组操作形式为 x y p。 当 p 为 1 时，如果第 x 变量和第 y 个变量可以相等，则输出 YES，并限制他们 相等；否则输出 NO，并忽略此次操作。 当 p 为 0 时，如果第 x 变量和第 y 个变量可以不相等，则输出 YES，并限制他 们不相等；否则输出 NO，并忽略此次操作。n<=1e5 x,y<=1e8。

题解：用并查集维护相等关系，用 set 维护不等关系。相等的放在一个并查集里面，这个集合开个 set，其中存放和它不等的并查集序号。 如果限制不相等，则在 set 中添加元素； 如果限制相等，通过启发式合并将两个并查集以及他们对应的 set 合并

T2AC代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<vector>
#include<cmath>
#include<set>
using namespace std;
inline int re_ad()
{
	int x=0;int f=1;char ch=getchar();
	while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9')x=x*10+ch-48,ch=getchar();
	return x*f;
}
int n,cnt=0;
int a[200010],fa[200010];
struct node{int x,y;int p;}ask[100010];
set<int> S[200010];set<int>::iterator it;
int find(int x){return x==fa[x]?fa[x]:fa[x]=find(fa[x]);}
int main()
{
	register int i,j,x,y;
	n=re_ad();
	for(i=1;i<=n;i++)
	{
	ask[i].x=re_ad();ask[i].y=re_ad();ask[i].p=re_ad();
	a[++cnt]=ask[i].x;a[++cnt]=ask[i].y;
	}
	sort(a+1,a+cnt+1);
	int q=unique(a+1,a+cnt+1)-a;
	for(i=1;i<=n;i++)
	{ask[i].x=lower_bound(a+1,a+q,ask[i].x)-a;ask[i].y=lower_bound(a+1,a+q,ask[i].y)-a;}
	--q;
	for(i=1;i<=q;++i)fa[i]=i;
	for(i=1;i<=n;i++)
	{
	x=ask[i].x;y=ask[i].y;
	x=find(x);y=find(y);
	if(ask[i].p)
	{
	it=S[x].lower_bound(y);
	if(it!=S[x].end()&&(*it)==y){printf("NO\n");continue;}
	printf("YES\n");if(x==y){continue;}
	if(S[x].size()<S[y].size())swap(x,y);
	fa[y]=x;
	for(it=S[y].begin();it!=S[y].end();++it)
	{S[x].insert(*it);S[*it].erase(y);S[*it].insert(x);}
	}
	else
	{
	if(x==y){printf("NO\n");continue;}
	printf("YES\n");S[x].insert(y);S[y].insert(x);
	}
	}
	
}

T3

题目大意：定义一个长度为奇数的区间的值为其所包含的的元素的中位数。 现给出 n 个数的集合 a，求将所有长度为奇数的区间的值排序后，第 K 大的值为 多少。n<=1e5。

题解：二分答案 t，统计中位数大于等于 t 的区间有多少个。 设 a[i]为前 i 个数中有 a[i]个数>=t，若奇数区间[l,r]的中位数>=t，则 (a[r]-a[l-1])*2>r-l+1,即(a[r]*2-r)>(a[l-1]*2-l+1)。 设 b[i]=a[i]*2-i，统计每个 b[i]有多少个 b[j](j和i奇偶性不同)。

T3AC代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<vector>
#include<cmath>
#include<bitset>
using namespace std;
inline int re_ad()
{
	int x=0,f=1;char ch=getchar();
	while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9')x=x*10+ch-48,ch=getchar();
	return x*f;
}
int n,k,m;
int a[100010],b[100010],num[100010],c[100010],d[100010];
int t[100010];
priority_queue<int> dgd,xgd;
inline void change(int x){while(x<=m)++t[x],x+=x&(-x);}
inline long long ask(int x){long long ret=0;while(x)ret+=t[x],x-=x&(-x);return ret;}
inline long long check(int x)
{
	register int i;
	long long ret=0;c[0]=0;
	for(i=1;i<=n;i++)c[i]=c[i-1]+(a[i]>=x);
	for(i=0;i<=n;i++)c[i]=2*c[i]-i,d[i+1]=c[i];
	sort(d+1,d+n+2);
	int q=unique(d+1,d+n+2)-d;
	for(i=0;i<=n;i++)c[i]=lower_bound(d+1,d+q,c[i])-d;
	m=q-1;
	memset(t,0,sizeof(t));
	for(i=0;i<=n;i++)
	{if(i&1)change(c[i]);else ret+=ask(c[i]-1);}
	memset(t,0,sizeof(t));
	for(i=0;i<=n;i++)
	{if(!(i&1))change(c[i]);else ret+=ask(c[i]-1);}
	return ret;
}
int main()
{
	register int i,j,x,l,r,mid,ans;
	n=re_ad();k=re_ad();
	for(i=1;i<=n;i++)a[i]=b[i]=re_ad();
	sort(b+1,b+n+1);
	int q=unique(b+1,b+n+1)-b;
	for(i=1;i<=n;i++)
	{a[i]=lower_bound(b+1,b+q,a[i])-b;}
	--q;
	l=1;r=q;
	while(l<r)
	{
	mid=(l+r)>>1;++mid;
	x=check(mid);
	if(x==k){cout<<b[mid]<<endl;return 0;}
	if(x>k){l=mid;}else r=mid-1;
	}
	cout<<b[l]<<endl;
}

T4

题目大意：一个立方体，是由 1*1*1 的小立方体拼成的。三条边的长度分别为 A, B, C。将立方体沿着3个顶点所在的平面切开，问将会切到多少个 1*1*1 的小立方体？答案mod 1e9+7。

题解：题解证明太可怕，式子玄学，自己看代码。

T4AC代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<queue>
#include<vector>
#include<cmath>
#include<bitset>
using namespace std;
inline long long re_ad()
{
	long long x=0;int f=1;char ch=getchar();
	while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9')x=x*10+ch-48,ch=getchar();
	return x*f;
}
const int mo=1000000007;
int T;
long long A,B,C,ans,mul;
long long gcd(long long x,long long y){return y?gcd(y,x%y):x;}
int main()
{
	T=re_ad();
	while(T--)
	{
	A=re_ad(),B=re_ad();C=re_ad();mul=gcd(gcd(A,B),C)%mo;
	A/=mul;B/=mul;C/=mul;
	ans=((gcd(A,B)%mo)*(gcd(A,C)%mo)%mo)*(gcd(B,C)%mo)%mo;
	mul=mul*mul%mo;
	A%=mo;B%=mo;C%=mo;
	ans=-ans+(A*B%mo)+(A*C%mo)+(B*C%mo);
	ans%=mo;ans=ans*mul%mo;
	if(ans&1)ans+=mo;
	ans>>=1;
	cout<<ans<<endl;
	}
}