东莞市2014年市队选拔赛 分组 构造二分图最小覆盖方案

xiaoxiao2021-02-28  102

题意

分析

显然判断前缀和后缀是否相等可以用hash。我在比赛的时候打了个贪心,每次都选比较大的组,水了59分。

其实可以把每个前缀和后缀都当成一个点,建立二分图,把每个字符串当成一条边,那么问题就转变成了求一个二分图最小覆盖的方案。

方法如下: 首先把所有匹配边和匹配点求出来,然后对于每个连通块讨论:若该连通块没有未匹配点,则把该连通块某一侧的点全部选掉即可。否则就从任意一个未匹配点开始,搜索出一棵交替树(未匹配边->匹配边->未匹配边…),那么交替树上深度为偶数的点就都可以选掉。显然每个点对应一条匹配边,如果某条匹配边两端都是奇数点,则会出现奇环,与二分图矛盾。 但考虑到一棵交替树未必能找到一个连通块的所有点,那么每找到一棵交替树,就把树上的点删掉,再继续重复就好了。这样的话就可以保证每一步都是正确的。

代码

#include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<algorithm> #include<map> #include<vector> #include<queue> using namespace std; typedef long long LL; const int N=20005; const int MOD=1000000007; const int inf=0x3f3f3f3f; int n,m,cnt,last[N],dis[N],s,t,cur[N*2],ans,sz,f1,f[N]; struct edge{int to,c,next,use;}e[N]; queue<int> que; vector<int> vec[N]; map<int,int> suf,pre; bool ty[N*2],cho[N],vis[N],use[N]; char ch[1005]; void addedge(int u,int v,int c) { e[++cnt].to=v;e[cnt].c=c;e[cnt].next=last[u];last[u]=cnt; e[++cnt].to=u;e[cnt].c=0;e[cnt].next=last[v];last[v]=cnt; } bool bfs() { for (int i=s;i<=t;i++) dis[i]=0; while (!que.empty()) que.pop(); dis[s]=1;que.push(s); while (!que.empty()) { int u=que.front();que.pop(); for (int i=last[u];i;i=e[i].next) if (e[i].c&&!dis[e[i].to]) { dis[e[i].to]=dis[u]+1; if (e[i].to==t) return 1; que.push(e[i].to); } } return 0; } int dfs(int x,int maxf) { if (x==t||!maxf) return maxf; int ret=0; for (int &i=cur[x];i;i=e[i].next) if (e[i].c&&dis[e[i].to]==dis[x]+1) { int f=dfs(e[i].to,min(e[i].c,maxf-ret)); e[i].c-=f; e[i^1].c+=f; ret+=f; if (maxf==ret) break; } return ret; } void dinic() { while (bfs()) { for (int i=s;i<=t;i++) cur[i]=last[i]; ans+=dfs(s,inf); } } void print(int x) { int size=0; for (vector<int>::iterator it=vec[x].begin();it!=vec[x].end();it++) if (!use[*it]) size++; if (!size) return; printf("%d ",size); for (vector<int>::iterator it=vec[x].begin();it!=vec[x].end();it++) if (!use[*it]) use[*it]=1,printf("%d ",*it); putchar('\n'); } void solve(int x,int dis) { if (dis) print(x); vis[x]=1; for (int i=last[x];i;i=e[i].next) if (!vis[e[i].to]&&e[i].use==dis) solve(e[i].to,dis^1); } void build() { vis[s]=vis[t]=1; for (int i=2;i<=cnt;i+=2) if (!ty[e[i^1].to]&&ty[e[i].to]&&!e[i].c) cho[e[i^1].to]=1,cho[e[i].to]=1,e[i].use=e[i^1].use=1; for (int i=1;i<=sz;i++) if (!vis[i]&&!cho[i]) solve(i,0); for (int i=1;i<=sz;i++) if (!vis[i]&&cho[i]&&!ty[i]) print(i); } int main() { freopen("group.in","r",stdin);freopen("group.out","w",stdout); scanf("%d%d",&n,&m); cnt=1; for (int i=1;i<=n;i++) { scanf("%s",ch); int len=strlen(ch),x=0,y=0; for (int j=0;j<m;j++) { x=((LL)x*27%MOD+ch[j]-'A'+1)%MOD; y=((LL)y*27%MOD+ch[len-1-j]-'A'+1)%MOD; } if (!pre[x]) pre[x]=++sz,ty[sz]=0; if (!suf[y]) suf[y]=++sz,ty[sz]=1; int p=pre[x],q=suf[y]; vec[q].push_back(i);vec[p].push_back(i); addedge(p,q,1); } s=0;t=sz+1; for (int i=1;i<=sz;i++) if (!ty[i]) addedge(s,i,1); else addedge(i,t,1); dinic(); printf("%d\n",ans); build(); return 0; }
转载请注明原文地址: https://www.6miu.com/read-43105.html

最新回复(0)