2017 CCPC Qinhuangdao Site
A. Balloon Robot
假设机器人$0$时刻位于$0$号位置，那么每个气球所需的时间为$(s_a-b)\bmod m$。
将所有气球按这个时间排序，枚举每个气球的时间作为偏移量，得出最优解即可。
时间复杂度$O(p\log p)$。
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=1000010;
int T,n,m,p,i,s[N];ll f[N],init,ans;
int main(){
scanf("%d",&T);
while(T--){
scanf("%d%d%d",&n,&m,&p);
for(i=1;i<=n;i++)scanf("%d",&s[i]);
init=0;
for(i=1;i<=p;i++){
int a,b;
scanf("%d%d",&a,&b);
f[i]=((s[a]-b)%m+m)%m;
init+=f[i];
}
ans=init;
sort(f+1,f+p+1);
for(i=1;i<=p;i++)ans=min(ans,init-1LL*p*f[i]+1LL*(i-1)*m);
printf("%lld\n",ans);
}
}
B. Expected Waiting Time
设$f_n$表示长度为$n$的合法括号序列个数，则奇数为$0$，偶数为卡特兰数。
分母一定是$f_{2n}$，而对于分子，可以考虑每个位置作为左右括号的贡献，假设是第$i$个位置作为左括号，那么枚举与其配对的右括号的距离$j$，则中间的方案数为$f_{j-1}$，剩下部分的方案数为$f_{2n-j-1}$，故总方案数为$\sum_{j=1}^{2n-i}f_{j-1}f_{2n-j-1}$。
注意到$i$只影响上式的求和上界，故直接计算出每个的值即可。
时间复杂度$O(n)$。
#include<cstdio>
typedef long long ll;
const int N=5000000;
int n,P,b,A,B,T,i,a[N],sum,ans;
int h[N],inv[N];
inline ll po(ll a,ll b){
ll t=1;
for(;b;b>>=1,a=a*a%P)if(b&1)t=t*a%P;
return t;
}
inline int C(int n){
if(n&1)return 0;
return h[n>>1];
}
int main(){
scanf("%d",&T);
while(T--){
scanf("%d%d%d%d%d",&n,&P,&b,&A,&B);
h[0]=1;
h[1]=1;
inv[0]=inv[1]=1;
for(i=2;i<=n+1;i++)inv[i]=1LL*(P-inv[P%i])*(P/i)%P;
for(i=2;i<=n;i++)h[i]=1LL*h[i-1]*(i*4-2)%P*inv[i+1]%P;
n*=2;
for(i=1;i<=n;i++){
b=(1LL*b*A+B)%P;
a[i]=(1LL*a[i-1]+1LL*b+1)%P;
}
sum=0;
ans=0;
for(i=1;i<n;i++){
//n-x = i
//x=n-i
sum=(1LL*C(i-1)*C(n-i-1)+sum)%P;
ans=(1LL*(P-a[n-i])*sum+ans)%P;
ans=(1LL*a[i+1]*sum+ans)%P;
}
printf("%d\n",1LL*ans*po(C(n),P-2)%P);
}
}
C. Crusaders Quest
爆搜所有可行方案即可。
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<string>
#include<ctype.h>
#include<math.h>
#include<set>
#include<map>
#include<vector>
#include<queue>
#include<bitset>
#include<algorithm>
#include<time.h>
using namespace std;
void fre() { freopen("c://test//input.in", "r", stdin); freopen("c://test//output.out", "w", stdout); }
#define MS(x, y) memset(x, y, sizeof(x))
#define ls o<<1
#define rs o<<1|1
typedef long long LL;
typedef unsigned long long UL;
typedef unsigned int UI;
template <class T1, class T2>inline void gmax(T1 &a, T2 b) { if (b > a)a = b; }
template <class T1, class T2>inline void gmin(T1 &a, T2 b) { if (b < a)a = b; }
const int N = 0, M = 0, Z = 1e9 + 7, inf = 0x3f3f3f3f;
template <class T1, class T2>inline void gadd(T1 &a, T2 b) { a = (a + b) % Z; }
int casenum, casei;
char w[4] = "gao";
int p[4];
int ANS;
void dfs(int o, string s, int val)
{
if(o == 2)
{
gmax(ANS, val + 1);
return;
}
int len = s.length();
string nxt = "";
for(int j = 0; j < len; ++j)if(s[j] != w[p[o]])nxt = nxt + s[j];
bool flag = 0;
for(int i = 0; i < len; ++i)if(s[i] == w[p[o]])
{
flag = (s[i] == s[i + 1] && s[i] == s[i + 2]);
break;
}
dfs(o + 1, nxt, val + flag);
}
int main()
{
scanf("%d", &casenum);
for (casei = 1; casei <= casenum; ++casei)
{
string s;
cin >> s;
for(int i = 0; i < 3; ++i)p[i] = i;
ANS = 0;
do
{
dfs(0, s, 0);
}while(next_permutation(p, p + 3));
printf("%d\n", ANS);
}
return 0;
}
/*
【trick&&吐槽】
【题意】
【分析】
【时间复杂度&&优化】
*/
D. Graph Generator
从后往前还原整个过程。
若图只有一个点，那么方案显然。
若图由多个连通块构成，那么每个连通块之间是独立的，分开处理即可。
否则图只剩下一个大小至少为$2$的连通块，那么最后一步操作必然选择一个到其它每个点都有边的点，然后将图继续分裂成若干连通块。
注意到对于$n$个点的连通块，边数会减少$n-1$，故迭代层数不超过$O(\sqrt{n})$。
时间复杂度$O(n\sqrt{n})$。
#include<cstdio>
#include<vector>
#include<algorithm>
using namespace std;
typedef pair<int,int>P;
const int N=100010,M=200010;
int Case,n,m,i,qa[N],qb[M],tmp[M];
int fin[N];//should reverse
vector<int>ans[N];
int stage;
bool err;
struct E{
int x,y;
}e[M];
int deg[N],g[N],v[M<<1],nxt[M<<1],ed;
int vis[N];
int cnt;
int q[N],cq;
int pos[M];
inline void add(int x,int y){
v[++ed]=y;
nxt[ed]=g[x];
g[x]=ed;
deg[x]++;
}
void dfs(int x){
if(vis[x])return;
vis[x]=cnt;
q[++cq]=x;
for(int i=g[x];i;i=nxt[i])dfs(v[i]);
}
void solve(int nl,int nr,int ml,int mr,vector<int>&old){
if(err)return;
//for(int i=nl;i<=nr;i++)printf("%d ",qa[i]);puts("");
if(nl==nr){
fin[++stage]=qa[nl];
old.push_back(qa[nl]);
return;
}
//find connect comp
int i,j;
cnt=ed=0;
for(i=nl;i<=nr;i++){
int x=qa[i];
deg[x]=g[x]=vis[x]=0;
}
for(i=ml;i<=mr;i++)add(e[qb[i]].x,e[qb[i]].y),add(e[qb[i]].y,e[qb[i]].x);
cq=0;
for(i=nl;i<=nr;i++){
int x=qa[i];
if(!vis[x]){
old.push_back(x);
cnt++;
dfs(x);
}
}
if(cnt==1){
int cv=nr-nl;
for(i=nl;i<=nr;i++)if(deg[qa[i]]==cv)break;
if(i>nr){
err=1;
return;
}
int x=qa[i];
swap(qa[nl],qa[i]);
int L=ml,R=ml-1;
for(i=ml;i<=mr;i++){
int u=e[qb[i]].x,v=e[qb[i]].y;
if(u==x||v==x)continue;
tmp[++R]=qb[i];
}
for(i=L;i<=R;i++)qb[i]=tmp[i];
fin[++stage]=x;
solve(nl+1,nr,L,R,ans[stage]);
}else{
vector<P>nson,mson;
int l=nl;
for(i=1;i<=cq;i=j){
int r=l-1;
for(j=i;j<=cq&&vis[q[i]]==vis[q[j]];j++){
qa[++r]=q[j];
}
nson.push_back(P(l,r));
l=r+1;
}
for(i=1;i<=cnt;i++)pos[i]=0;
pos[0]=ml-1;
for(i=ml;i<=mr;i++)pos[vis[e[qb[i]].x]]++;
for(i=1;i<=cnt;i++)pos[i]+=pos[i-1];
for(i=1;i<=cnt;i++)mson.push_back(P(pos[i-1]+1,pos[i]));
for(i=ml;i<=mr;i++){
int x=vis[e[qb[i]].x];
tmp[pos[x]--]=qb[i];
}
for(i=ml;i<=mr;i++)qb[i]=tmp[i];
int _cnt=cnt;
for(i=0;i<_cnt;i++){
ans[0].clear();
solve(nson[i].first,nson[i].second,mson[i].first,mson[i].second,ans[0]);
}
}
}
int main(){
scanf("%d",&Case);
while(Case--){
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)qa[i]=i;
for(i=1;i<=m;i++)qb[i]=i;
for(i=1;i<=m;i++)scanf("%d%d",&e[i].x,&e[i].y);
stage=err=0;
for(i=0;i<=n;i++)ans[i].clear();
solve(1,n,1,m,ans[0]);
if(err)puts("No");else{
puts("Yes");
for(i=n;i;i--){
printf("%d %d",fin[i],ans[i].size());
for(int j=0;j<ans[i].size();j++)printf(" %d",ans[i][j]);
puts("");
}
}
}
}
E. String of CCPC
设$f[i][j][k]$表示考虑前$i$个字符，目前购买了$j$次，与“CCPC”KMP的指针为$k$的最大净收益。
注意到最优解中$j$为个位数，令其不超过$9$即可。
时间复杂度$O(n)$。
#include<cstdio>
const int N=200010,M=10;
int T,n,m,i,j,k,t,ans,b[N],nxt[N],g[M][2],w[M][2];
int f[N][M][4];
char a[N];
inline void up(int&a,int b){a<b?(a=b):0;}
int main(){
m=4;
b[1]=0;
b[2]=0;
b[3]=1;
b[4]=0;
for(i=2;i<=m;nxt[i++]=j){
while(j&&b[j+1]!=b[i])j=nxt[j];
if(b[j+1]==b[i])j++;
}
for(i=0;i<m;i++){
for(j=0;j<2;j++){
int k=i,o=0;
while(k&&b[k+1]!=j)k=nxt[k];
if(b[k+1]==j)k++;
if(k==m)k=nxt[k],o++;
g[i][j]=k;
w[i][j]=o;
}
}
scanf("%d",&T);
while(T--){
scanf("%d%s",&n,a+1);
for(i=1;i<=n;i++)a[i]=a[i]=='P';
for(i=0;i<=n;i++)for(j=0;j<M;j++)for(k=0;k<m;k++)f[i][j][k]=-10000000;
f[0][0][0]=0;
for(i=0;i<=n;i++)for(j=0;j<M;j++)for(k=0;k<m;k++){
if(j+1<M){
for(t=0;t<2;t++)up(f[i][j+1][g[k][t]],f[i][j][k]+w[k][t]-j);
}
if(i<n)up(f[i+1][j][g[k][a[i+1]]],f[i][j][k]+w[k][a[i+1]]);
}
ans=0;
for(j=0;j<M;j++)for(k=0;k<m;k++)up(ans,f[n][j][k]);
printf("%d\n",ans);
}
}
F. Getting Lost
留坑。
G. Numbers
在二进制下从高位到低位贪心考虑。
若答案这一位可以为$0$，那么全部填$0$，否则尽量填$1$。
import java.util.*;
import javax.swing.text.TabableView;
import java.io.*;
import java.math.*;
public class Main
{
static final int N = (int)1e5 + 10;
static Scanner cin = new Scanner(System.in);
static BigInteger v0 = BigInteger.valueOf(0);
static BigInteger v1 = BigInteger.valueOf(1);
static BigInteger v2 = BigInteger.valueOf(2);
static BigInteger b[] = new BigInteger [4000];
public static void main(String args[])
{
b[0] = v1;
for(int i = 1; i < 4000; ++i)b[i] = b[i - 1].multiply(v2);
int casenum = cin.nextInt();
for(int casei = 1; casei <= casenum; ++casei)
{
BigInteger n = cin.nextBigInteger();
BigInteger m = cin.nextBigInteger();
BigInteger top = n.add(m).subtract(v1).divide(m);
int w = 0;
while(b[w].compareTo(top) <= 0)++w;
BigInteger ans = v0;
for(int k = w; k >= 0; --k)
{
if(n.compareTo( b[k].subtract(v1).multiply(m) ) <= 0)
{
}
else
{
BigInteger g = n.divide(b[k]);
if(g.compareTo(m) > 0)g = m;
n = n.subtract(g.multiply(b[k]));
ans = ans.add(b[k]);
}
}
System.out.println(ans);
}
}
}
H. Prime Set
按奇偶建立二分图匹配模型，先忽略$1$求出最大匹配，再考虑$1$继续求最大匹配，每个匹配都将贡献$2$。
对于剩下的数，再检查能否与已选的数配对，每次配对贡献$1$。
注意特殊处理$1$内部的匹配。
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<string>
#include<ctype.h>
#include<math.h>
#include<set>
#include<map>
#include<vector>
#include<queue>
#include<bitset>
#include<algorithm>
#include<time.h>
using namespace std;
void fre() { freopen("c://test//input.in", "r", stdin); freopen("c://test//output.out", "w", stdout); }
#define MS(x, y) memset(x, y, sizeof(x))
#define ls o<<1
#define rs o<<1|1
typedef long long LL;
typedef unsigned long long UL;
typedef unsigned int UI;
template <class T1, class T2>inline void gmax(T1 &a, T2 b) { if (b > a)a = b; }
template <class T1, class T2>inline void gmin(T1 &a, T2 b) { if (b < a)a = b; }
const int N = 3030, M = N * N, Z = 1e9 + 7, inf = 0x3f3f3f3f;
//edge_num = ???
template <class T1, class T2>inline void gadd(T1 &a, T2 b) { a = (a + b) % Z; }
int casenum, casei;
bool is_prime[(int)2e6 + 10];
int nn, K;
int o[N], e[N];
int ST, ED, pone;
int first[N], id;
int w[M], cap[M], nxt[M];
void ins(int x, int y, int cap_)
{
w[++id] = y;
cap[id] = cap_;
nxt[id] = first[x];
first[x] = id;
w[++ id] = x;
cap[id] = 0;
nxt[id] = first[y];
first[y] = id;
}
int d[N];
bool bfs()
{
//MS(d, -1);
for(int i = 0; i <= ED; ++i)d[i] = -1;
queue<int> q; q.push(ST); d[ST] = 0;
while(! q.empty()){
int x = q.front(); q.pop();
for(int z = first[x]; z; z = nxt[z]) if(cap[z]){
int y = w[z];
if(d[y] == -1){
d[y] = d[x] + 1;
q.push(y);
if(y == ED) return 1;
}
}
}
return 0;
}
int dfs(int x, int all)
{
if(x == ED) return all;
int use = 0;
for(int z = first[x]; z; z = nxt[z] ) if(cap[z]){
int y = w[z];
if(d[y] == d[x] + 1){
int tmp = dfs(y, min(cap[z], all - use));
cap[z] -= tmp;
cap[z ^ 1] += tmp;
use += tmp;
if(use == all) break;
}
}
if(use == 0) d[x] = -1;
return use;
}
int dinic()
{
int ret = 0;
while(bfs()) ret += dfs(ST, inf);
return ret;
}
void prime_init()
{
int top = 2e6;
MS(is_prime, 1); is_prime[0] = is_prime[1] = 0;
for(int i = 2; i <= top; ++i)if(is_prime[i])
{
for(int j = i + i; j <= top; j += i)
{
is_prime[j] = 0;
}
}
}
int main()
{
prime_init();
scanf("%d", &casenum);
for (casei = 1; casei <= casenum; ++casei)
{
scanf("%d%d", &nn, &K);
int odd = 0;
int even = 0;
int one = 0;
for(int i = 1; i <= nn; ++i)
{
int x; scanf("%d", &x);
if(x == 1)++one;
else if(x & 1)o[++odd] = x;
else e[++even] = x;
}o[0] = 1;
pone = 0; ST = odd + even + 1; ED = ST + 1;
for(int i = 0; i <= ED; ++i)first[i] = 0; id = 1;
for(int i = 0; i <= odd; ++i)
{
for(int j = 1; j <= even; ++j)if(is_prime[o[i] + e[j]])
{
ins(i, odd + j, 1);
}
}
for(int i = 1; i <= odd; ++i)ins(ST, i, 1);
for(int i = 1; i <= even; ++i)ins(odd + i, ED, 1);
int ORI_ONE = one;
int now = dinic();
//
//printf("two-two pair = %d\n", now);
//
while(one)
{
ins(ST, pone, 1);
if(dinic())
{
++now;
--one;
}
else break;
}
now += one / 2;
one %= 2;
if(K <= now)
{
printf("%d\n", K * 2);
continue;
}
//K > now
K -= now;
int sum = 0;
for(int i = (1 - one); i <= odd; ++i)if(i == 0 || cap[first[i]] == 0)////
{
bool flag = 0;
for(int j = 1; j <= even; ++j)if(is_prime[o[i] + e[j]])
{
//
//printf("single odd for even = %d %d\n", o[i], e[j]);
//
sum += 1;
flag = 1;
break;
}
if(!flag && i == 0 && ORI_ONE > 1)sum += 1;
}
//
//for(int i = 1; i <= even; ++i)printf("cap[%d] = %d\n", e[i], cap[first[odd + i]]);
//
for(int i = 1; i <= even; ++i)if(cap[first[odd + i]] == 1)////
{
int st = ORI_ONE ? 0 : 1;
for(int j = st; j <= odd; ++j)if(is_prime[e[i] + o[j]])
{
//
//printf("single even for one = %d %d\n", e[i], o[j]);
//
sum += 1;
break;
}
}
int ans = now * 2 + min(sum, K);
printf("%d\n", ans);
}
return 0;
}
/*
【trick&&吐槽】
【题意】
【分析】
【时间复杂度&&优化】
6 2
1 1 1 1 2 2
6 3
1 1 1 1 2 2
6 2
1 1 10 10 10 10
7 2
1 1 1 1 1 2 2
7 3
1 1 1 1 1 2 2
7 4
1 1 1 1 1 2 2
*/
I. Triangulation
留坑。
J. Tree Equation
留坑。
K. Diversity and Variance
留坑。
L. One-Dimensional Maze
答案为$\min([2,m]中R的个数,[m,n-1]中L的个数)$。
#include<cstdio>
#include<algorithm>
using namespace std;
const int N=1000010;
int T,n,m,i,A,B;char a[N];
int main(){
scanf("%d",&T);
while(T--){
scanf("%d%d%s",&n,&m,a+1);
A=B=0;
for(i=2;i<=m;i++)if(a[i]=='R')A++;
for(i=m;i<n;i++)if(a[i]=='L')B++;
printf("%d\n",min(A,B));
}
}
M. Safest Buildings
下一轮安全区的圆心可行范围可以用圆表示，圆交求出概率即可。
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<string>
#include<ctype.h>
#include<math.h>
#include<set>
#include<map>
#include<vector>
#include<queue>
#include<bitset>
#include<algorithm>
#include<time.h>
using namespace std;
void fre() { freopen("c://test//input.in", "r", stdin); freopen("c://test//output.out", "w", stdout); }
#define MS(x, y) memset(x, y, sizeof(x))
#define ls o<<1
#define rs o<<1|1
typedef long long LL;
typedef unsigned long long UL;
typedef unsigned int UI;
template <class T1, class T2>inline void gmax(T1 &a, T2 b) { if (b > a)a = b; }
template <class T1, class T2>inline void gmin(T1 &a, T2 b) { if (b < a)a = b; }
const int N = 0, M = 0, Z = 1e9 + 7, inf = 0x3f3f3f3f;
template <class T1, class T2>inline void gadd(T1 &a, T2 b) { a = (a + b) % Z; }
int casenum, casei;
const double eps = 1e-8, PI = acos(-1.0);
struct circle
{
long double x, y, r;
}p1, p2;
long double sqr(long double x)
{
return x * x;
}
long double dis(circle a, circle b)
{
return sqrt(sqr(a.x - b.x) + sqr(a.y - b.y));
}
long double solve(circle a, circle b)
{
long double d = dis(a, b);
if(d >= a.r + b.r) return 0;
if(d <= fabs(a.r - b.r)){
long double r = min(a.r, b.r);
return PI * r * r;
}
long double ang1 = acos((a.r * a.r + d * d - b.r * b.r) / (2. * a.r * d)) * 2;
long double ang2 = acos((b.r * b.r + d * d - a.r * a.r) / (2. * b.r * d)) * 2;
double tmp1 = 0.5 * ang1 * a.r * a.r - 0.5 * a.r * a.r * sin(ang1);
double tmp2 = 0.5 * ang2 * b.r * b.r - 0.5 * b.r * b.r * sin(ang2);
long double ret = tmp1 + tmp2;
return ret;
}
circle a[110];
int b[110];
double area[110];
int sgn(double x)
{
if(fabs(x) < eps) return 0;
return x > 0 ? 1 : -1;
}
int main()
{
double r;
int n;
scanf("%d", &casenum);
for (casei = 1; casei <= casenum; ++casei)
{
double R;
scanf("%d%lf%lf", &n, &R, &r);
a[0].x = 0, a[0].y = 0; a[0].r = R- r;
double ans = 0;
for(int i = 1; i <= n; i ++){
double x, y;
scanf("%lf%lf", &x, &y);
a[i].x = x; a[i].y = y;
a[i].r = r;
area[i] = solve(a[0], a[i]);
gmax(ans, area[i]);
}
int num = 0;
for(int i = 1; i <= n; i ++){
if(sgn(area[i] - ans) == 0){
b[++ num] = i;
}
}
printf("%d\n", num);
for(int i = 1; i < num; i ++) printf("%d ", b[i]); printf("%d\n", b[num]);
}
return 0;
}
/*
【trick&&吐槽】
【题意】
【分析】
【时间复杂度&&优化】
*/
