2020-2021 “Orz Panda” Cup Programming Contest
A. Accordion Artist And Orz Pandas
题意
给定(n)个相邻的矩形，矩形的宽和高给定，问最大的子矩形面积
分析
单调栈板题，维护出每个矩形左边第一个比他低的点和右边第一个比它低的点，两点距离就是宽，高就是当前矩形的高，维护宽X高的最大值即可
代码
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll rd(){
ll x;
scanf("%lld",&x);
return x;
}
int main(){
int n = rd();
vector<ll> h(n + 1),w(n + 1);
vector<ll> sum(n + 1);
vector<int> lpos(n + 1),rpos(n + 1);
for(int i = 1;i <= n;i++){
h[i] = rd();
w[i] = rd();
sum[i] = sum[i - 1] + w[i];
}
stack<int> st;
for(int i = 1;i <= n;i++){
while(!st.empty() && h[st.top()] > h[i]) {
rpos[st.top()] = i;
st.pop();
}
st.push(i);
}
while(!st.empty()) {
rpos[st.top()] = n + 1;
st.pop();
}
for(int i = n;i >= 1;i--){
while(!st.empty() && h[st.top()] > h[i]) {
lpos[st.top()] = i;
st.pop();
}
st.push(i);
}
while(!st.empty()){
lpos[st.top()] = 0;
st.pop();
}
ll ans = 0;
for(int i = 1;i <= n;i++){
ll W = sum[rpos[i] - 1] - sum[lpos[i]];
ans = max(ans,W * h[i]);
}
cout << ans;
}
C. Closestools of Orz Pandas
题意
模拟题
给定(n)个房间，定义i,j号房间距离(|i - j|)
操作1：选择一个空房间进入，要求该空房间得到的和其他非空房间生成的距离排序后字典序最小
操作2：将第x次操作的房间变为空房间
分析
贪心策略：
把空房间看成线段，每次找到最长的线段，取最长线段的中点。
通常维护线段可以用set实现。再用另一个线段维护非空房间，可以快速得到左右点位置。
代码
#include<bits/stdc++.h>
#define re register
using namespace std;
typedef long long ll;
int rd(){
int x = 0;
char ch = getchar();
while(ch < '0' || ch > '9') {
ch = getchar();
}
while(ch >= '0' && ch <= '9') {
x = x * 10 + ch - '0';
ch = getchar();
}
return x;
}
struct Seg{
int l,r; //[l,r]
Seg(int _l = 0,int _r = 0){
l = _l;
r = _r;
}
int d() const {
return (r - l) / 2 + 1;
}
bool operator < (const Seg& b) const{
int d1 = d();
int d2 = b.d();
return (d1 > d2) || (d1 == d2 && l < b.l);
}
};
int main(){
int n,m;
while(~scanf("%d%d",&n,&m)){
set<Seg> s;
s.insert(Seg(1,n));
set<int> st;
vector<int> v(m + 1);
for(re int i = 1;i <= m;i++){
int op = rd();
if(op == 1) {
auto best = s.begin();
int bestd = best -> d();
if(!st.empty() && *st.begin() > bestd) {
best = s.lower_bound(Seg(1,*st.begin() - 1));
bestd = (*st.begin() - 1) / 2;
}
if(!st.empty() && n - *st.rbegin() > bestd) {
best = s.lower_bound(Seg(*st.rbegin() + 1,n));
bestd = (n - *st.rbegin()) / 2;
}
if(st.empty())
best = s.begin();
int pos = best -> l + bestd - 1;
if(best -> l == 1)
pos = 1;
else if(best -> r == n)
pos = n;
printf("%d
",pos);
v[i] = pos;
st.insert(pos);
Seg L = Seg(best -> l,pos - 1);
Seg R = Seg(pos + 1,best -> r);
s.erase(best);
if(L.l <= L.r)
s.insert(L);
if(R.l <= R.r)
s.insert(R);
}
else {
int x = rd();
int pos = v[x];
auto it = st.lower_bound(pos);
auto itt = it;
int L = (it != st.begin() ? *(--it) + 1 : 1);
itt++;
int R = (itt != st.end() ? *itt - 1 : n);
if(L <= pos - 1)
s.erase(Seg(L,pos - 1));
if(pos + 1 <= R)
s.erase(Seg(pos + 1,R));
s.insert(Seg(L,R));
st.erase(pos);
}
}
}
}
D. Data Structure Master and Orz Pandas
题意
大概意思就是求期望(如有不对的地方请指正)
对一颗树中的某点，找到根到这个点的路径把这条路径上的边都变成实边（LCT中，某个点的儿子中只能连一条实边），类似LCT的access操作
求稳态下拉一个点，边从虚边变为实边的期望次数
分析
考虑答案组成，枚举每个点，那么这个点到路径上所有点被修改的期望次数的和就是答案（期望的可加性）
即
[ans = frac{1}{n}sum_u sum_{v是u到root的路径上的点} E(v和fa的边不是实边)
]
考虑(E),显然只和(v)和(fa_v)有关，简单推导（直觉）就可以得到(E = 1 - frac{siz(v)}{siz(fa_v) - 1})
于是变换求和顺序即可得
[ans = frac{1}{n} sum_u siz(u) (1 - frac{siz(u)}{siz(fa_u) - 1})
]
dfs即可
代码
#include<bits/stdc++.h>
#define re register
using namespace std;
typedef long long ll;
int rd(){
int x = 0;
char ch = getchar();
while(ch < '0' || ch > '9') {
ch = getchar();
}
while(ch >= '0' && ch <= '9') {
x = x * 10 + ch - '0';
ch = getchar();
}
return x;
}
const int maxn = 1e5 + 5;
const int MOD = 998244353;
ll ksm(ll a,ll b = MOD - 2,ll m = MOD){
ll base = a;
ll ans = 1;
while(b){
if(b & 1) ans *= base,ans %= MOD;
base *= base;
base %= MOD;
b >>= 1;
}
return ans;
}
int Fa[maxn];
int siz[maxn];
vector<int> e[maxn];
int ans;
void dfs(int u,int fa){
Fa[u] = fa;
siz[u] = 1;
for(auto v:e[u]) {
dfs(v,u);
siz[u] += siz[v];
}
}
int main(){
int n = rd();
for(int i = 2;i <= n;i++){
int x = rd();
e[x].push_back(i);
}
dfs(1,0);
for(int u = 2;u <= n;u++){
ans += (ll)siz[u] * (1 - (ll)siz[u] * ksm(siz[Fa[u]] - 1 + MOD) % MOD + MOD) % MOD,ans %= MOD;
}
ans = (ll) ans * ksm(n) % MOD;
cout << ans;
}
E. Encryption of Orz Pandas
题意
经过转化可以把问题变为：
给定长度(n)的01序列，对序列做k次前缀和，求k次前缀和后的序列
(n leq 1e5)
分析
做k次前缀和可以用生成函数看成乘上((1 + x+ x^2...+x^n)^k)
即若有序列((a_0,a_1...a_n))，那么k次前缀和后的(a_i = [x^i](a_0 + a_1 ...+a_n) imes (1+x+x^2...x^n)^k)
做法可以用题解写的泰勒展开后转化为组合数再FFT，复杂度(O(nlogn))
不过我看不懂，比较朴素的想法是对(FFT)做快速幂复杂度(O(nlognlogk)) 考虑到时限5s可以冲一冲
我觉得FFT调试起来有点麻烦，事实上当时没调出来，赛后发现用NTT即可，01序列做1次不会爆模数
代码
#include<bits/stdc++.h>
#define re register
using namespace std;
typedef long long ll;
ll rd(){
ll x = 0;
char ch = getchar();
while(ch < '0' || ch > '9') {
ch = getchar();
}
while(ch >= '0' && ch <= '9') {
x = x * 10 + ch - '0';
ch = getchar();
}
return x;
}
const int maxn = 2e5 + 5;
class NTTClass{
public:
static const int MAXL=22;
static const int MAXN=1<<MAXL;
static const int root=3;
static const int MOD=998244353;
int rev[MAXN];
int fast_pow(int a,int b){
int ans=1;
while(b){
if(b&1)ans=1ll*ans*a%MOD;
a=1ll*a*a%MOD;
b>>=1;
}
return ans;
}
void transform(int n,int *t,int typ){
for(int i=0;i<n;i++)
if(i<rev[i])swap(t[i],t[rev[i]]);
for(int step=1;step<n;step<<=1){
int gn=fast_pow(root,(MOD-1)/(step<<1));
for(int i=0;i<n;i+=(step<<1)){
int g=1;
for(int j=0;j<step;j++,g=1ll*g*gn%MOD){
int x=t[i+j],y=1ll*g*t[i+j+step]%MOD;
t[i+j]=(x+y)%MOD;
t[i+j+step]=(x-y+MOD)%MOD;
}
}
}
if(typ==1)return;
for(int i=1;i<n/2;i++)swap(t[i],t[n-i]);
int inv=fast_pow(n,MOD-2);
for(int i=0;i<n;i++)t[i]=1ll*t[i]*inv%MOD;
}
void ntt(int p,int *A,int *B,int *C){
transform(p,A,1);
transform(p,B,1);
for(int i=0;i<p;i++)C[i]=1ll*A[i]*B[i]%MOD;
transform(p,C,-1);
}
void mul(int *A,int *B,int *C,int n,int m) {
int p=1,l=0;
while(p<=n+m)p<<=1,l++;
//printf("n = %d, m = %d
",n,m);
for (int i=n+1;i<p;i++) A[i] = 0;
for (int i=m+1;i<p;i++) B[i] = 0;
//for (int i=0;i<p;i++) printf("%d ",A[i]);puts("");
//for (int i=0;i<p;i++) printf("%d ",B[i]);puts("");
for(int i=0;i<p;i++) rev[i]=(rev[i>>1]>>1)|((i&1)<<(l-1));
ntt(p,A,B,C);
//puts("C:");for (int i=0;i<p;i++) printf("%d ",C[i]);puts("");
}
}NTT;
int _[maxn << 1],f[maxn << 1];
int ans[maxn << 1],p[maxn << 1],ttt[maxn << 1];
int n,m,mm;
inline void ksm(int *a,ll b){
ans[0] = 1;
for(re int i = 1;i < n;i++)
ans[i] = 0;
while(b){
if(b & 1){
for(re int i = 0;i < n;i++)
p[i] = a[i];
NTT.mul(ans,a,_,n,n);
for(re int i = 0;i < n;i++)
a[i] = p[i],ans[i] = (_[i] & 1 ? 1 : 0);
}
for(re int i = 0;i < n;i++)
p[i] = a[i];
NTT.mul(a,p,_,n,n);
for(re int i = 0;i < n;i++)
a[i] = (_[i] & 1 ? 1 : 0);
//for(re int i = 0;i <= mm;i++) {
// int tmp = (int)(ans[i].x / n + 0.49);
// if(tmp & 1) ans[i].x = 1;
// else ans[i].x = 0;
//}
b >>= 1;
}
}
int main(){
n = rd();
mm = n;
m = n;
ll k = rd();
int cnt = 0;
int M = mm;
while(M) M /= 2,cnt++;
if((1ll << cnt) != mm) cnt++;
k = (k - 1) % (1 << cnt);
k++;
vector<ll> v(mm);
vector<ll> res(mm);
for(int i = 0;i < mm;i++)
v[i] = rd();
for(int i = 0;i < n;i++)
f[i] = 1;
ksm(f,k);
//for(int i = 0;i <= mm;i++)
// printf("%d ",ans[i]);
for(int i = 0;i <= 17;i++){
for(re int j = 0;j < mm;j++){
int tttt = ((v[j] >> i) & 1);
ttt[j] = ans[j];
f[j] = tttt;
}
for(int j = mm;j < n;j++){
ttt[j] = f[j] = 0;
}
/*
if(i == 2) for(int j = 0;j < n;j++){
cout << ttt[j].x << ' ' << f[j].x << '
';
}*/
NTT.mul(f,ttt,_,n,n);
for(re int j = 0;j < mm;j++){
int ok = _[j] & 1 ? 1 : 0;
//if(ok) cout << "i = " << i << "j = " << j << '
';
if(ok) res[j] |= (1ll << i);
}
}
for(int i = 0;i < mm;i++)
if(i) putchar(' '),cout << res[i];
else cout << res[i];
}
H. Hamming Code and Orz Pandas
模拟题，因为我没有参与此题，所以无法提供解释
代码
#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll rd(){
ll x = 0;
int f = 1;
char ch = getchar();
while(ch < '0' || ch > '9') {
if(ch == '-') f = -1;
ch = getchar();
}
while(ch >= '0' && ch <= '9') {
x = x * 10 + ch - '0';
ch = getchar();
}
return x * f;
}
int k;
string d;
vector<int>a;
vector<int>have[20];
void init()
{
int num=1;
a.push_back(0);
for(int i=1;i<=17;i++){
a.push_back(num);
num=num*2;
}
int p=1;
for(int i=3;i<=1024*64*2;i++){
p=i;
for(int j=16;j>=1;j--){
if(i==a[j]){
p=0;
break;
}
}
if(p==0)continue;
int hh=0;
while(p){
for(int j=17;j>=1;j--){
hh=0;
while(p>=a[j]){
p-=a[j];
hh=1;
}
if(hh)have[j].push_back(i);
if(p==0)break;
}
}
}
}
int main()
{
init();
while(~scanf("%d",&k)){
cin>>d;
int num=0;
for(int i=0;i<d.size();i++){
num=num^d[i];
}
if(num==0){
int wa=0;
for(int i=1;i<=k;i++){
int pos=0;
//cout<<i<<":"<<endl;
for(auto x:have[i]){
if(x>=a[k+1])break;
//cout<<x<<endl;
pos=pos^d[x];
}
if(pos!=d[a[i]]){
wa=1;
break;
}
}
if(wa){
printf("broken
");
}else{
printf("good
");
}
}else{
vector<int>ans;
for(int i=1;i<=k;i++){
int pos=0;
for(auto x:have[i]){
if(x>=a[k+1])break;
pos=pos^d[x];
}
if(pos!=d[a[i]]){
ans.push_back(a[i]);
}
}
//cout<<ans.size()<<endl;
if(ans.size()==1){
printf("d(%d) is changed
",ans[0]);
}else{
int sum=0;
for(auto x:ans){
sum+=x;
}
printf("d(%d) is changed
",sum);
}
}
}
}
/*
4
1011101110111011
4
1011101110111010
4
1011101110111110
*/
I. Irregular Shape of Orz Pandas
题意
顺序给出点，求多边形面积
[3 leq n leq 1e3\
-1e9 leq x,y leq 1e9
]
分析
一开始用double算，发现爆了，因为给的是整点，所以直接用__int128，事实上longlong即可，最后除以二可以人为加小数 这题花了很长时间，实在是不应该
代码
#include<bits/stdc++.h>
using namespace std;
typedef __int128 ll;
ll rd(){
ll x = 0;
int f = 1;
char ch = getchar();
while(ch < '0' || ch > '9') {
if(ch == '-') f = -1;
ch = getchar();
}
while(ch >= '0' && ch <= '9') {
x = x * 10 + ch - '0';
ch = getchar();
}
return x * f;
}
void Put(ll x){
if(x == 0) return;
if(x < 10) {
putchar(x + '0');
return;
}
Put(x / 10);
putchar(x % 10 + '0');
}
struct Point{
ll x;
ll y;
Point(){}
Point(ll _x,ll _y):x(_x),y(_y){}
};
ll operator ^ (const Point &a,const Point &b){
return a.x * b.y - a.y * b.x;
}
Point operator - (const Point &a,const Point &b){
return Point(a.x - b.x,a.y - b.y);
}
struct Pol{
int n;
Point p[1005];
void input(int _n){
n = _n;
for(int i = 0;i < n;i++){
p[i].x = rd();
p[i].y = rd();
}
}
ll getarea(){
ll sum = 0;
for(int i = 0;i < n;i++){
sum += ((p[i]) ^ (p[(i + 1) % n]));
}
return sum;
}
};
int main(){
int n;
while(~scanf("%d",&n)){
Pol poly;
poly.input(n);
ll tmp = poly.getarea();
if(tmp < 0) tmp *= -1;
if(tmp % 2 == 1) {
Put(tmp / 2);
putchar('.');
putchar('5');
putchar('0');
}
else {
Put(tmp / 2);
putchar('.');
putchar('0');
putchar('0');
}
puts("");
}
}