2019.01.27【NOIP提高组】模拟 B 组

JZOJ 1273 袁绍的刁难

题目

用3的非负整数幂组成一个数，使这个数第$k$小，询问这个数

---

分析

那么首先这些幂是互不受影响的，而且递增排列，这可以想到位运算也是这样操作的，那么$k$的二进制就是表示答案是否选择的二进制状态，那么把$k$进行二进制，然后再用3的非负整数幂表达

---

代码

#include <cstdio>
#include <cctype>
#define rr register
using namespace std;
long long p3[41];
inline signed iut(){
	rr int ans=0; rr char c=getchar();
	while (!isdigit(c)) c=getchar();
	while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();
	return ans;
}
inline void print(long long ans){
	if (ans>9) print(ans/10);
	putchar(ans%10+48);
}
signed main(){
	freopen("recruitment.in","r",stdin);
	freopen("recruitment.out","w",stdout); 
	p3[0]=1;
	for (rr int i=1;i<32;++i) p3[i]=p3[i-1]*3;
	for (rr int t=iut();t;--t){
		rr long long ans=0;
		for (rr int n=iut(),len=0;n;n>>=1,++len)
			if (n&1) ans+=p3[len];
		print(ans); putchar(10);
	}
	return 0; 
} 

---

JZOJ 5230 队列统计

题目

现在有n个人要排成一列，编号为1->n 。但由于一些不明原因的关系，人与人之间可能存在一些矛盾关系，具体有m条矛盾关系(u,v),表示编号为u的人想要排在编号为v的人前面。要使得队伍和谐，最多不能违背k条矛盾关系（即不能有超过k条矛盾关系(u,v),满足最后v排在了u前面）。问有多少合法的排列。答案对10^9+7取模。

---

分析

那么这道题是一个状压dp的题目，设$f[s][k]$表示选定的集合为$s$，矛盾关系$k$条的排列数量，那么${\sum }_{i=0}^{k-矛盾数}f[s][与s产生的矛盾+i]+=f[上一个位置][i]$，但是问题是这样时间复杂度是$O(2^n{n}^{2}k)$，铁定会超时，问题是这样还是会超时,那么要把1提取出来，优化为$O(能过)$。

---

代码

#include <cstdio>
#include <cctype>
#define rr register
using namespace std;
const int mod=1000000007;
unsigned short n,m,k,cnt[1050001];
int all,f[1050001][21],cont[21],ans;
inline signed iut(){
	rr int ans=0; rr char c=getchar();
	while (!isdigit(c)) c=getchar();
	while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();
	return ans;
}
inline signed p(int x){if (x>=mod) x-=mod; return x;}
signed main(){
	freopen("count.in","r",stdin);
	freopen("count.out","w",stdout);
	all=1<<(n=iut()); m=iut(); k=iut(); f[0][0]=1;
	for (rr int i=1;i<=m;++i) cont[iut()-1]|=1<<(iut()-1);//在它后面的，用一个集合表示
	for (rr int i=1;i<all;++i) cnt[i]=cnt[i-(-i&i)]+1;//提前提取出1的个数，用lowbit实现
	for (rr int i=1;i<all;++i){
		for (rr int x=i,w=-x&x;x;x-=w,w=-x&x)
		for (rr int *p1=f[i]+cnt[cont[cnt[w-1]]&i],*p2=f[i^w];p1<=f[i]+k;++p1,++p2) *p1=p(*p1+*p2);
	}
	for (rr int i=0;i<=k;++i) ans=p(f[all-1][i]+ans);
	printf("%d",ans);
	return 0;
}

---

JZOJ 5231 序列问题

题目

∑ l = 1 n ∑ r = l n m a x { a [ l ∼ r ] } × m i n { a [ l ∼ r ] } \sum_{l=1}^n\sum_{r=l}^nmax\{a[l\sim r]\}\times min\{a[l\sim r]\} l=1∑n​r=l∑n​max{a[l∼r]}×min{a[l∼r]}

---

分析(我太菜了)

---

代码

#include <cstdio>
#include <cctype>
#define rr register
#define min(a,b) ((a)<(b)?(a):(b))
#define max(a,b) ((a)>(b)?(a):(b))
using namespace std;
const int mod=1000000007,N=500010;
int ans,a[N],maxx[N],minx[N],sh[N],ss[N],sl[N],n;
inline signed iut(){
	rr int ans=0; rr char c=getchar();
	while (!isdigit(c)) c=getchar();
	while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();
	return ans;
}
inline void print(int ans){
	if (ans>9) print(ans/10);
	putchar(ans%10+48);
}
inline void merg(int l,int r){
	if (l==r) ans=(1ll*a[l]*a[r]+ans)%mod;
	if (l>=r) return; rr int mid=(l+r)>>1,l1=mid,l2=mid;
	maxx[mid]=minx[mid]=a[mid]; ss[l-1]=sl[l-1]=sh[l-1]=0;
	for (rr int i=mid-1;i>=l;--i)//向两边扩展
	    maxx[i]=max(maxx[i+1],a[i]),minx[i]=min(minx[i+1],a[i]);
	for (rr int i=mid+1;i<=r;++i)
	    maxx[i]=max(maxx[i-1],a[i]),minx[i]=min(minx[i-1],a[i]);
	for (rr int i=l;i<=mid;++i)//记录前缀和
	    ss[i]=(1ll*maxx[i]*minx[i]+ss[i-1])%mod,
	    sl[i]=(sl[i-1]+minx[i])%mod,sh[i]=(sh[i-1]+maxx[i])%mod;
	for (rr int i=mid;i<=r;++i){
		while (l1>=l&&maxx[l1]<maxx[i]) --l1;
		while (l2>=l&&minx[l2]>minx[i]) --l2;
		ans=(ans+ss[min(l1,l2)])%mod;//直接用maxx[i]*minx[i]的前缀和实现
		ans=(1ll*(mid-max(l1,l2))*maxx[i]%mod*minx[i]%mod+ans)%mod;//确定的区间乘上它的长度
		if (l1<=l2) ans=(1ll*maxx[i]*(sl[l2]-sl[l1])+ans)%mod;//不知道最小值，用前缀和实现
		else ans=(1ll*minx[i]*(sh[l1]-sh[l2])+ans)%mod;//不知道最大值，用前缀和实现
	}
	merg(l,mid-1),merg(mid+1,r);
}
signed main(){
	freopen("seq.in","r",stdin);
	freopen("seq.out","w",stdout);
	n=iut();
	for (rr int i=1;i<=n;++i) a[i]=iut();
	merg(1,n); printf("%d",ans);
}

---

后续

沉迷奥赛，日渐变菜，力不从心，心花怒放，放荡不羁