XVI Open Cup named after E.V. Pankratiev. GP of Siberia

A. Passage

枚举两个点，看看删掉之后剩下的图是否是二分图。

#include <bits/stdc++.h>
using namespace std ;

const int MAXN = 205 ;

vector < int > G[MAXN] ;
int vis[MAXN] , col[MAXN] ;
int n ;

int dfs ( int u ) {
	for ( int i = 0 ; i < G[u].size () ; ++ i ) {
		int v = G[u][i] ;
		if ( vis[v] == 0 ) continue ;
		if ( !col[v] ) {
			col[v] = 3 - col[u] ;
			if ( !dfs ( v ) ) return 0 ;
		}
		if ( col[u] + col[v] != 3 ) return 0 ;
	}
	return 1 ;
}

int check () {
	for ( int i = 1 ; i <= n ; ++ i ) {
		col[i] = 0 ;
	}
	for ( int i = 1 ; i <= n ; ++ i ) {
		if ( col[i] || !vis[i] ) continue ;
		col[i] = 1 ;
		if ( !dfs ( i ) ) return 0 ;
	}
	return 1 ;
}

void solve () {
	for ( int i = 1 ; i <= n ; ++ i ) {
		int x , y ;
		scanf ( "%d" , &x ) ;
		for ( int j = 0 ; j < x ; ++ j ) {
			scanf ( "%d" , &y ) ;
			G[i].push_back ( y ) ;
			G[y].push_back ( i ) ;
		}
		vis[i] = 1 ;
	}
	if ( n <= 3 ) {
		printf ( "Hurrah!\n" ) ;
		return ;
	}
	for ( int i = 1 ; i <= n ; ++ i ) {
		for ( int j = i + 1 ; j <= n ; ++ j ) {
			vis[i] = vis[j] = 0 ;
			if ( check () ) {
				printf ( "Hurrah!\n" ) ;
				return ;
			}
			vis[i] = vis[j] = 1 ;
		}
	}
	printf ( "Fired.\n" ) ;
}

int main () {
	freopen ( "input.txt" , "r" , stdin ) ;
	freopen ( "output.txt" , "w" , stdout ) ;
	while ( ~scanf ( "%d" , &n ) ) solve () ;
	return 0 ;
}

　　

B. Files list

按题意模拟。

#include <bits/stdc++.h>
using namespace std ;

const int MAXN = 10005 ;

map < string , int > mp ;
map < int , string > mp2 ;
char s[MAXN] , p[MAXN] ;
int n ;

void solve () {
	mp.clear () ;
	mp2.clear () ;
	int cnt = 0 ;
	for ( int i = 0 ; i < n ; ++ i ) {
		scanf ( "%s" , s ) ;
		int t = 0 , f = 0 ;
		for ( int j = 0 ; s[j] ; ++ j ) {
			if ( f ) p[t ++] = s[j] ;
			else if ( s[j] == '.' ) f = 1 ;
		}
		p[t] = 0 ;
		if ( mp.count ( p ) ) mp[p] ++ ;
		else {
			mp[p] ++ ;
			mp2[++ cnt] = p ;
		}
	}
	for ( int i = 1 ; i <= cnt ; ++ i ) {
		cout << mp2[i] << ": " << mp[mp2[i]] << endl ;
	}
}

int main () {
	freopen ( "input.txt" , "r" , stdin ) ;
	freopen ( "output.txt" , "w" , stdout ) ;
	while ( ~scanf ( "%d" , &n ) ) solve () ;
	return 0 ;
}

　　

C. Graph optimization

将所有1类限制的边加入，分块bitset判定限制2是否都满足即可。

时间复杂度$O(\frac{nm}{64})$。

#include<algorithm>
#include<cstdio>
#include<bitset>
#include<set>
#include<ctime>
using namespace std;
typedef bitset<4096>B;
typedef unsigned long long ll;
const int N=300010;

int n,m,i,j,k,x,y,z,ans[N][2],quailty;
int g[N],G[N],v[N],nxt[N],ed;

int vis[N],q[N],h,t,cnt,f[N],d[N];
B dp[100010];//dp[x] : this can reach x

set<int>SET[100010];

struct E{int x,y,z;}e[N];

const int BUF=5000000;
char Buf[BUF],*buf=Buf;
inline void read(int&a){for(a=0;*buf<48;buf++);while(*buf>47)a=a*10+*buf++-48;}

inline bool cmp(const E&a,const E&b){return a.z<b.z;}
inline void add(int x,int y){
  v[++ed]=y;nxt[ed]=g[x];g[x]=ed;
}
inline void add2(int x,int y){
  v[++ed]=y;nxt[ed]=G[x];G[x]=ed;
}
void dfs1(int x){
  vis[x]=1;
  for(int i=g[x];i;i=nxt[i])if(!vis[v[i]])dfs1(v[i]);
  q[++t]=x;
}
void dfs2(int x,int y){
  vis[x]=0;f[x]=y;
  for(int i=G[x];i;i=nxt[i])if(vis[v[i]])dfs2(v[i],y);
}
inline bool solve(int L,int R){
  if ( clock () > 2.95 * CLOCKS_PER_SEC ) return 1 ;
  int z=e[L].z;
  int i;
  for(i=1;i<=cnt;i++){
    dp[i].reset();
    //if ( clock () > 2.95 * CLOCKS_PER_SEC ) return 1 ;
  }
  for(i=max(z<<12,1);(i>>12)==z;i++){
    dp[f[i]][i&4095]=1;
  }
  for(i=1;i<=cnt;i++){
    int x=q[i];
    for(int j=g[x];j;j=nxt[j]){
      dp[v[j]]|=dp[x];
      //if ( clock () > 2.95 * CLOCKS_PER_SEC ) return 1 ;
    }
  }
  for(i=L;i<=R;i++){
    int x=e[i].x,y=e[i].y;
    //printf("-> %d %d %d %llu\n",x,y,z,dp[f[y]]);
    if(dp[f[y]][x&4095])return 0;
  }
  return 1;
}
int main (){
  freopen ( "input.txt" , "r" , stdin ) ;
  freopen ( "output.txt" , "w" , stdout ) ;
  fread(Buf,1,BUF,stdin);
  read(n),read(m);
  for(i=1;i<=m;i++){
    //x=i,y=i+1;
    read(x),read(y);
    ans[i][0]=x,ans[i][1]=y;
    add(x,y);
    add2(y,x);
  }

  for(i=1;i<=n;i++)if(!vis[i])dfs1(i);
  for(i=n;i;i--)if(vis[q[i]])dfs2(q[i],++cnt);

  for(ed=0,i=1;i<=cnt;i++)g[i]=0;
  for(i=1;i<=m;i++){
    x=ans[i][0],y=ans[i][1];
    if(f[x]==f[y])continue;
    if(SET[f[x]].find(f[y])!=SET[f[x]].end())continue;
    SET[f[x]].insert(f[y]);
    add(f[x],f[y]);
    d[f[y]]++;
  }
  for(h=1,t=0,i=1;i<=cnt;i++)if(!d[i])q[++t]=i;
  while(h<=t)for(i=g[q[h++]];i;i=nxt[i])if(!(--d[v[i]]))q[++t]=v[i];
  //now q is topo sequence of SCC

  read(quailty);
  for(i=1;i<=quailty;i++){
    //e[i].x=i+1;
    //e[i].y=i;
    read(e[i].x),read(e[i].y);
    e[i].z=e[i].x>>12;
  }
  sort(e+1,e+quailty+1,cmp);
  for(i=1;i<=quailty;i=j){
    for(j=i;j<=quailty&&e[i].z==e[j].z;j++);
    if(!solve(i,j-1))return puts("NO"),0;
  }
  puts("YES");
  printf("%d\n",m);
  for(i=1;i<=m;i++)printf("%d %d\n",ans[i][0],ans[i][1]);
  return 0 ;
}

　　

D. Housing payments

设$f[i]$表示前$i$个月付清，且第$i$个月进行了交易的最小代价，那么因为债务指数级增长，所以可用决策只有$O(\log n)$项。

#include <bits/stdc++.h>
using namespace std ;

typedef pair < int , int > pii ;

const int MAXN = 100005 ;
const int INF = 0x3f3f3f3f ;

int s[MAXN] , x[MAXN] , p[MAXN] ;
double dp[MAXN] ;
int n ;

void solve () {
	for ( int i = 1 ; i <= n ; ++ i ) {
		scanf ( "%d%d%d" , &s[i] , &x[i] , &p[i] ) ;
		dp[i] = 1e18 ;
	}
	dp[0] = 0 ;
	for ( int i = 0 ; i <= n ; ++ i ) {
		double sum1 = 0 ;
		for ( int j = i + 1 ; j <= i + 100 && j <= n ; ++ j ) {
			dp[j] = min ( dp[j] , dp[i] + sum1 + s[j] + x[j] ) ;
			sum1 = ( sum1 + s[j] ) * ( 1 + 0.01 * p[j] ) ;
		}
	}
	printf ( "%.10f\n" , dp[n] ) ;
}

int main () {
	freopen ( "input.txt" , "r" , stdin ) ;
	freopen ( "output.txt" , "w" , stdout ) ;
	while ( ~scanf ( "%d" , &n ) ) solve () ;
	return 0 ;
}

　　

E. Arithmetic expressions

$f[i][j]$表示长度为$i$的值为$j$的表达式个数，然后枚举后面接上什么转移即可。

#include <bits/stdc++.h>
using namespace std ;
const int mod=1e9+7;
int n,m,p;
int dp[55][222];
int len[333];
void up(int &x,int y){x+=y;if(x>=mod)x-=mod;}
int main () {
	freopen ( "input.txt" , "r" , stdin ) ;
	freopen ( "output.txt" , "w" , stdout ) ;
	while(scanf("%d%d%d",&n,&m,&p)!=EOF){
		memset(dp,0,sizeof dp);
		for(int i=0;i<m;i++){
			int x=i;
			if(!x){len[i]=1;continue;}
			int l=0;
			while(x)l++,x/=10;
			len[i]=l;
		}
		for(int i=1;i<=n;i++){
			for(int j=0;j<m;j++){
				if(len[j]==i)up(dp[i][j],1);
				if(i>=3)up(dp[i][j],dp[i-2][j]);
			}
			//if(i==2)up(dp[i][0],1);
			for(int j=0;j<m;j++){
				for(int k=1;k+1<i;k++){
					int nj,ni=i-k-1;
					for(int nm=0;nm<m;nm++){
						for(int ty=0;ty<2;ty++){
							if(!ty)nj=(j+nm)%m;
							else nj=(j-nm+m)%m;
							if(ni>2)up(dp[i][nj],1LL*dp[k][j]*dp[ni-2][nm]%mod);
							if(len[nm]==ni)up(dp[i][nj],dp[k][j]);
						}
					}
				}
			}
		}
		printf("%d\n",dp[n][p]);
	}
        return 0;
}

　　

F. Sputnik

留坑。

G. Voting

设$f[i][j]$表示$i$的子树里选举情况为$j$的最小代价，转移则用另一个$dp[i][a][b][c]$表示考虑了前$i$个儿子，三个人选票各自为$a,b,c$的最小代价来转移。

#include <bits/stdc++.h>
using namespace std ;
const int mod=1e9+7;
const int Maxn=10050,M=10,Inf=1e9;
typedef pair<int,int>pi;
int n,m,K;
vector<int>G[Maxn];
int f[Maxn][4];
int g[23][23][23][23];
int ori[Maxn];
vector<pi>ans;
vector<pi>res[Maxn][4];
pi pre[23][23][23][23];
void init(int cs){
	for(int i=0;i<21;i++){
		for(int j=0;j<21;j++){
			for(int k=0;k<21;k++)
				g[cs][i][j][k]=Inf;
		}
	}
}
struct State{
	int x,y,z;
	State(){}
	State(int x,int y,int z):x(x),y(y),z(z){}
}st[Maxn][4];
int get(int a,int b,int c){
	vector<pi>v;
	v.push_back(pi(a,1));
	v.push_back(pi(b,2));
	v.push_back(pi(c,3));
	sort(v.begin(),v.end());
	reverse(v.begin(),v.end());
	if(v[0].first==v[1].first)return 0;
	return v[0].second;
}
void dfs(int u){
	if(u>m){
		for(int i=0;i<4;i++){
			if(i>K)f[u][i]=Inf;
			else {
				res[u][i].push_back(pi(u,i));
				f[u][i]=i==ori[u]?0:1;
			}
		}
		return ;
	}
	for(int i=0;i<G[u].size();i++){
		int v=G[u][i];
		dfs(v);
	}
	int n=G[u].size();
	init(0);
	g[0][0][0][0]=0;
	for(int i=0;i<n;i++){
		int v=G[u][i];
		init(i+1);
		int it[4];
		for(it[0]=0;it[0]<=20;it[0]++)
		for(it[1]=0;it[1]<=20;it[1]++)
		for(it[2]=0;it[2]<=20;it[2]++){
			if(g[i][it[0]][it[1]][it[2]]==Inf)continue;
			int w=g[i][it[0]][it[1]][it[2]];
			for(int ty=0;ty<4;ty++){
				if(f[v][ty]==Inf)continue;
				if(ty)it[ty-1]++;
				int nw=w+f[v][ty];
				int &t=g[i+1][it[0]][it[1]][it[2]];
				if(nw<t){
					t=nw;
					pre[i+1][it[0]][it[1]][it[2]]=pi(v,ty);
				}
				if(ty)it[ty-1]--;
			}
		}
	}
	for(int i=0;i<4;i++)f[u][i]=Inf;
	for(int i=0;i<=20;i++){
		for(int j=0;j<=20;j++){
			for(int k=0;k<=20;k++){
				if(g[n][i][j][k]==Inf)continue;
				int ty=get(i,j,k);
				if(f[u][ty]>g[n][i][j][k]){
					f[u][ty]=g[n][i][j][k];
					st[u][ty]=State(i,j,k);
				}
			}
		}
	}
	for(int i=0;i<4;i++){
		if(f[u][i]==Inf)continue;
		res[u][i].clear();
		State tmp=st[u][i];
		for(int cur=n;cur>=1;cur--){
			res[u][i].push_back(pre[cur][tmp.x][tmp.y][tmp.z]);
			int ty=pre[cur][tmp.x][tmp.y][tmp.z].second;
			if(ty==1)tmp.x--;
			if(ty==2)tmp.y--;
			if(ty==3)tmp.z--;
		}
		reverse(res[u][i].begin(),res[u][i].end());
	}
	//printf("u=%d\n",u);
	//for(int i=0;i<4;i++)printf("%d ",f[u][i]);puts("");
}
void dfs2(int u,int ty){
	//printf("u=%d ty=%d\n",u,ty);
	if(u>m){
		if(ori[u]!=ty)ans.push_back(pi(u,ty));
		return;
	}
	//printf("sz=%d\n",(int)res[u][1].size());
	for(int i=0;i<G[u].size();i++){
		//printf("v=%d\n",G[u][i]);
		dfs2(G[u][i],res[u][ty][i].second);
	}
}
int main () {
	freopen ( "input.txt" , "r" , stdin ) ;
	freopen ( "output.txt" , "w" , stdout ) ;
	scanf("%d%d%d",&n,&m,&K);
	for(int i=1;i<=n;i++)scanf("%d",&ori[i+m]);
	for(int i=1;i<=m;i++){
		int k;scanf("%d",&k);
		for(int j=0;j<k;j++){
			int x;scanf("%d",&x);
			if(x<0)x=-x;
			else x=m+x;
			G[i].push_back(x);
		}
	}
	dfs(1);
	//puts("ok1");
	dfs2(1,1);
	printf("%d\n",(int)ans.size());
	for(int i=0;i<ans.size();i++)printf("%d %d\n",ans[i].first-m,ans[i].second);
        return 0;
}

　　

H. Novice urbanist

枚举每个点和每个区间，那么它能贡献的距离是一段区间，差分前缀和即可。

#include <bits/stdc++.h>
using namespace std ;

int n,m,l,r,cnt,i,j;
int a[11111],f[4444444];
struct P{int l,r;P(){}P(int _l,int _r){l=_l,r=_r;}}b[1111],c[1111];
inline bool cmp(const P&a,const P&b){return a.l<b.l;}

int main () {
  freopen ( "input.txt" , "r" , stdin ) ;
  freopen ( "output.txt" , "w" , stdout ) ;
  scanf("%d%d",&n,&m);
  for(i=1;i<=n;i++)scanf("%d",&a[i]);
  for(i=1;i<=m;i++)scanf("%d%d",&b[i].l,&b[i].r);
  sort(b+1,b+m+1,cmp);
  l=1e9;
  r=-l;
  for(i=1;i<=m;i++){
    if(b[i].l>r+1){
      if(l<=r)c[++cnt]=P(l,r);
      l=b[i].l;
    }
    r=max(r,b[i].r);
  }
  if(l<=r)c[++cnt]=P(l,r);
  for(i=1;i<=n;i++){
    for(j=1;j<=cnt;j++){
      f[c[j].l-a[i]+2000000]++;
      f[c[j].r-a[i]+2000000+1]--;
    }
  }
  for(i=1;i<4000000;i++)f[i]+=f[i-1];
  int ans=-1,ans2;
  for(i=0;i<2000000;i++){
    int t=max(f[2000000+i],f[2000000-i]);
    if(t>ans)ans=t,ans2=i;
  }
  printf("%d %d",ans2,ans);
  return 0 ;
}

　　

I. Rangefinder

留坑。

J. Hive

求出方向向量$(x,y)$，顺时针旋转$60$°后是$(x+y,-x)$。

#include <bits/stdc++.h>
using namespace std ;
const int Maxn=100020;
int n,x,y,a,b;
int main () {
	freopen ( "input.txt" , "r" , stdin ) ;
	freopen ( "output.txt" , "w" , stdout ) ;
	scanf("%d%d%d",&n,&x,&y);
	while(n--){
	  scanf("%d%d",&a,&b);
	  a-=x,b-=y;
	  printf("%d %d\n",x+a+b,y-a);
        }
        return 0;
}

　　

K. Side effects

每次加入一条边之后暴力染色，然后将经过的边删掉，均摊复杂度$O(n+m)$。

#include <bits/stdc++.h>
using namespace std ;
const int Maxn=100020;
int n,m,q,ans;
set<int>G[Maxn];
int col[Maxn];
void dfs(int u){
	//printf("u=%d\n",u);
	col[u]=1;
	for(set<int>::iterator it=G[u].begin();it!=G[u].end();){
		//printf("iv=%d\n",*it);
		if(!col[*it]){
			dfs(*it);
			ans++;
		}
		G[u].erase(it++);
	}
}
void solve(){
	ans=0;
	for(int i=1;i<=m;i++){
		int x;scanf("%d",&x);
		col[x]=1;
		ans++;
	}
	while(q--){
		//puts("ok");
		int u,v;scanf("%d%d",&u,&v);
		if(u!=v)G[v].insert(u);
		if(col[v]==1)dfs(v);
		printf("%d\n",ans);
	}
}
int main () {
	freopen ( "input.txt" , "r" , stdin ) ;
	freopen ( "output.txt" , "w" , stdout ) ;
	while ( ~scanf ( "%d%d%d" , &n,&m,&q ) ) solve () ;
	return 0 ;
}

　　

L. Cypher

留坑。

---

总结：

• 在分块bitset时，应尽量设大bitset的大小，减少bitset的操作次数，这样常数反而小。