2017 ACM/ICPC Asia Regional Qingdao Online

A Apple

不会Java和高精度

B Bomberman

unsolved

C The Dominator of Strings

取最长的串跑AC自动机，跑过的点需要标记。

#include <bits/stdc++.h>
#define Next next
using namespace std;
const int MAXN=1e5+5;
struct Trie
{
    int vis[MAXN];
    int next[MAXN][26],fail[MAXN],end[MAXN];
    int root,L;
    int newnode()
    {
        for(int i = 0;i < 26;i++)
            next[L][i] = -1;
        vis[L] = 0;
        end[L++] = 0;
        return L-1;
    }
    void init()
    {
        L = 0;
        root = newnode();
    }
    void insert(char buf[],int len)
    {
        int now = root;
        for(int i = 0;i < len;i++)
        {
            if(next[now][buf[i]-'a'] == -1)
                next[now][buf[i]-'a'] = newnode();
            now = next[now][buf[i]-'a'];
        }
        end[now]++;
    }
    void build()
    {
        queue<int>Q;
        fail[root] = root;
        for(int i = 0;i < 26;i++)
            if(next[root][i] == -1)
                next[root][i] = root;
            else
            {
                fail[next[root][i]] = root;
                Q.push(next[root][i]);
            }
        while( !Q.empty() )
        {
            int now = Q.front();
            Q.pop();
            for(int i = 0;i < 26;i++)
                if(next[now][i] == -1)
                    next[now][i] = next[fail[now]][i];
                else
                {
                    fail[next[now][i]]=next[fail[now]][i];
                    Q.push(next[now][i]);
                }
        }
    }
    int query(char buf[],int len)
    {
        int now = root;
        int res = 0;
        for(int i = 0;i < len;i++)
        {
            now = next[now][buf[i]-'a'];
            int temp = now;
            while( temp != root )
            {
                if(vis[temp]) break;
                vis[temp] = 1;
                res += end[temp];
                end[temp] = 0;
                temp = fail[temp];
            }
        }
        return res;
    }
}aho;

int n;
char tmp[MAXN];
char str[MAXN];
int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        int Len=0;
        aho.init();
        scanf("%d",&n);
        for(int i=0;i<n;i++){
            scanf("%s",str);
            int len=strlen(str);
            aho.insert(str,len);
            if(Len<=len){
                Len=len;
                strcpy(tmp,str);
            }
        }
        aho.build();
        int ans=aho.query(tmp,Len);
        if(ans==n){
            printf("%s
",tmp);
        }
        else printf("No
");
    }

    return 0;
}

D The Intersection

unsolved

E King's Visit

unsolved

F Pythagoras

本原勾股三角形，所有边互质。

1.有m,n互质，并且异号，则（n*n+m*m, 2*n*m, n*n-m*m）组成一组本原勾股三角形。

2.Stern-Brocot Tree中能够构造出所有互质的数对。

跑一下Stern-Brocot Tree，然后计数，但是2.3e8有点慢，勉强卡过去

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

typedef long long LL;
const LL inf=1e9;

LL n;
LL a[(1<<18)];
LL num[(1<<18)];
LL cnt;

inline void solve(int a,int b,int c,int d){
    int L=(a+c),R=(b+d);
    if(1LL*L*L+1LL*R*R>inf) return;
    if((R-L)%2) num[((max(2*R*L,R*R-L*L))%(1<<17))]++;
    solve(a,b,L,R); solve(L,R,c,d);
}

int main(){
    solve(0,1,1,1);
    LL t;
    scanf("%lld",&t);
    while(t--){
        LL ans=0;
        scanf("%d",&n);
        for(int i=0;i<(1<<n);i++)
            scanf("%d",&a[i]);
        for(int i=0;i<(1<<17);i++)
            ans+=num[i]*a[i%(1<<n)];
        printf("%lld
",ans);
    }

    return 0;
}

G Zuma

区间dp，先处理出连续相同的块。

1.普通合并。

2.两侧相同的情况，合并两侧。

3.两侧和中间合并，且需要保证两侧和不超过3个并且中间的个数只能是1。

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

const int inf=0x3f3f3f3f;

int n;
char s[205];

int tot;
int a[205];
int num[205];
int dp[205][205];

int main(){
    int t;
    scanf("%d",&t);
    int T=t;
    while(t--){
        tot=0;
        memset(num,0,sizeof(num));
        memset(dp,inf,sizeof(dp));
        scanf("%s",s+1);
        n=strlen(s+1);
        int cnt=1;
        for(int i=2;i<=n;i++){
            if(s[i]!=s[i-1]){
                cnt%=3;
                if(cnt==0) continue;
                tot++;
                a[tot]=s[i-1]-'0';
                num[tot]=cnt;
                cnt=1;
            }
            else cnt++;
            if(i==n){
                if(cnt>=3) continue;
                tot++;
                a[tot]=s[i]-'0';
                num[tot]=cnt;
            }
        }
        for(int i=1;i<=tot;i++){
            if(num[i]==1) dp[i][i]=2;
            else dp[i][i]=1;
        }
        for(int len=2;len<=tot;len++){
            for(int i=1;i<=tot-len+1;i++){
                int j=i+len-1;
                for(int k=i;k<j;k++){
                    dp[i][j]=min(dp[i][j],dp[i][k]+dp[k+1][j]);
                }
                if(a[i]!=a[j]) continue;
                dp[i][j]=min(dp[i][j],dp[i+1][j-1]+(num[i]+num[j]<3));
                if(num[i]+num[j]<4){
                    for(int k=i+1;k<j;k++){
                        if(num[k]==1&&a[k]==a[i]){
                            dp[i][j]=min(dp[i][j],dp[i+1][k-1]+dp[k+1][j-1]);
                        }
                    }
                }
            }
        }
        printf("Case #%d: %d
",T-t,dp[1][tot]);
    }

    return 0;
}

H Chinese Zodiac

直接模拟

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <map>
using namespace std;
map<string,int> mp;
int main(){
    mp["rat"]=11;
    mp["ox"]=10;
    mp["tiger"]=9;
    mp["rabbit"]=8;
    mp["dragon"]=7;
    mp["snake"]=6;
    mp["horse"]=5;
    mp["sheep"]=4;
    mp["monkey"]=3;
    mp["rooster"]=2;
    mp["dog"]=1;
    mp["pig"]=0;
    int t;
    cin>>t;
    while(t--){
        string s1,s2;
        cin>>s1>>s2;
        int ans=mp[s1]+12-mp[s2];
        if(ans>12) ans-=12;
        cout<<ans<<endl;
    }


    return 0;
}

I Smallest Minimum Cut

边权*1000+1，结果为最大流对1000取模。

#include <bits/stdc++.h>
using namespace std;
const int inf=0x3f3f3f3f;
const int MAXN=205;
int level[MAXN];
int iter[MAXN];
int head[MAXN],tot;
struct edge{
    int to,cap,Next;
} e[2005];
void init(){
    memset(head,-1,sizeof(head));
    tot=0;
}
void add(int from,int to,int cap){
    e[tot].Next=head[from];
    e[tot].to=to;
    e[tot].cap=cap;
    head[from]=tot;
    tot++;
}
void addedge(int from,int to,int cap){
    add(from,to,cap);
    add(to,from,0);
}
void bfs(int s){
    memset(level,-1,sizeof(level));
    queue<int> q;
    level[s]=0;
    q.push(s);
    while(!q.empty()){
        int v=q.front(); q.pop();
        for(int i=head[v];~i;i=e[i].Next){
            edge &ed=e[i];
            if(ed.cap>0&&level[ed.to]<0){
                level[ed.to]=level[v]+1;
                q.push(ed.to);
            }
        }
    }
}
int dfs(int v,int t,int f){
    if(v==t) return f;
    for(int &i=iter[v];~i;i=e[i].Next){
        edge &ed=e[i];
        if(ed.cap>0&&level[v]<level[ed.to]){
            int d=dfs(ed.to,t,min(f,ed.cap));
            if(d>0){
                ed.cap-=d;
                e[i^1].cap+=d;
                return d;
            }
        }
    }
    return 0;
}
int max_flow(int s,int t){
    int flow=0;
    while(1){
        bfs(s);
        if(level[t]<0) return flow;
        memcpy(iter,head,sizeof(iter));
        int f;
        while((f=dfs(s,t,inf))>0){
            flow+=f;
        }
    }
}

int n,m;

int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        init();
        int s,t;
        scanf("%d%d",&n,&m);
        scanf("%d%d",&s,&t);
        for(int i=1;i<=m;i++){
            int u,v,w;
            scanf("%d%d%d",&u,&v,&w);
            addedge(u-1,v-1,w*1000+1);
        }
        int ans=max_flow(s-1,t-1);
        printf("%d
",ans%1000);
    }

    return 0;
}

J Brute Force Sorting

链表维护当前剩余的点，队列维护可能被炸掉的点。

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

const int MAXN=1e5+5;

int n;
int a[MAXN];
int in[MAXN];
int vis[MAXN];
int pre[MAXN];
int Next[MAXN];

void solve(int x){
    pre[Next[x]]=pre[x];
    Next[pre[x]]=Next[x];
    vis[x]=1;
}

int main(){
    int t;
    scanf("%d",&t);
    while(t--){
        queue<int> q;
        scanf("%d",&n);
        for(int i=1;i<=n;i++){
            scanf("%d",&a[i]);
            Next[i]=i+1;
            pre[i]=i-1;
            vis[i]=0;
            in[i]=1;
            q.push(n+1-i);
        }
        while(!q.empty()){
            int u=q.front();
            q.pop(); in[u]=0;
            if(vis[u]) continue;
            int pos=u;
            if(pre[pos]>=1&&a[pos]<a[pre[pos]]){
                if(Next[pos]<=n&&!in[Next[pos]]){
                    q.push(Next[pos]);
                    in[Next[pos]]=1;
                }
                while(pre[pos]>=1&&a[pos]<a[pre[pos]]){
                    solve(pos);
                    pos=pre[pos];
                }
                if(pos>=1) solve(pos);
            }
        }
        int num=0;
        for(int i=1;i<=n;i++)
            if(!vis[i]) num++;
        printf("%d
",num);
        for(int i=1;i<=n;i++)
            if(!vis[i]) printf("%d ",a[i]);
        printf("
");
    }

    return 0;
}

K A Cubic number and A Cubic Number

立方差公式为(a-b)(a^2+ab+b^2)，所以a,b需要相邻，预处理出所有的值，然后直接查找即可。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
vector<LL> v;
int main(){
    int t;
    cin>>t;
    for(LL i=1;i<=577350;i++){
        LL a=i,b=i+1;
        v.push_back(a*a+a*b+b*b);
    }
    sort(v.begin(),v.end());
    while(t--){
        LL x;
        cin>>x;
        LL id=lower_bound(v.begin(),v.end(),x)-v.begin();
        if(v[id]==x) cout<<"YES"<<endl;
        else cout<<"NO"<<endl;
    }


    return 0;
}