topcoder srm 435 div1

problem1 link

遍历未被删除的叶子结点即可。

problem2 link

首先,将所有的蛋白质原子编号,设为$[0,m-1]$,每个原子可能对应多个长度为3的$ACGT$。设$n$为DNA串的长度。用$g[i][j]$表示从$i$开始匹配数字$j$后的最小位置,$0\leq i \leq n-1,0\leq j \leq m-1$。

然后就是动态规划.设$dp[i][j]$表示位置$i$之前最后匹配的数字为$j$的方案数。对所有的$0\leq t < m$有$dp[g[i][t]][t]+=dp[i][j] $

problem3 link

首先将所有节点分成若干个$group$,每个$group$里的节点可以互相到达。然后一个一个$group$进行处理。

对于某个$group$,建一个新的二分图$G$。对于其中的一个点$x$,若存在不在该$group$中的一点$y$使得存在边$y$到$x$那么在新图$G$中连一条$x$到$y$的边。那么该$group$中需要作为新的$division$的根结点的个数为$|group|-maxMatch(G)$。

 

code for problem1

import java.util.*;
import java.math.*;
import static java.lang.Math.*;

public class CellRemoval {

	List<List<Integer>> lists=null;

	public int cellsLeft(int[] parent, int deletedCell) {
		final int n=parent.length;
		lists=new ArrayList<>();
		for(int i=0;i<n;++i) {
			lists.add(new ArrayList<>());
		}
		int root=-1;
		for(int i=0;i<n;++i) {
			if(parent[i]==-1) {
				root=i;
			}
			else {
				lists.get(parent[i]).add(i);
			}
		}
		return dfs(root,deletedCell);
	}

	int dfs(int u,int d) {
		if(u==d) {
			return 0;
		}
		if(lists.get(u).size()==0) {
			return 1;
		}
		int result=0;
		for(int i=0;i<lists.get(u).size();++i) {
			result+=dfs(lists.get(u).get(i),d);
		}
		return result;
	}
}

  


code for problem2

import java.util.ArrayList;
import java.util.HashMap;
import java.util.List;
import java.util.Map;

/**
 * Created by mcl on 2017/9/27.
 */
public class DNADeletion {

	final static int MOD=1000000007;

	public int differentProteins(String[] DNASequence, String[] codonTable) {
		StringBuilder sb=new StringBuilder();
		for(int i=0;i<DNASequence.length;++i) {
			sb.append(DNASequence[i]);
		}
		final String S=sb.toString();
		final int n=S.length();
		int[][] nxt=new int[n+1][4];
		nxt[n][0]=nxt[n][1]=nxt[n][2]=nxt[n][3]=-1;
		for(int i=n-1;i>=0;--i) {
			for(int j=0;j<4;++j) {
				nxt[i][j]=nxt[i+1][j];
			}
			nxt[i][getIndex(S.charAt(i))]=i;
		}
		List<List<String>> lists=init(codonTable);
		final int m=lists.size();
		int[][] g=new int[n][m];
		for(int i=0;i<n;++i) {
			for(int j=0;j<m;++j) {
				g[i][j]=-1;
				for(int k=0;k<lists.get(j).size();++k) {
					String s=lists.get(j).get(k);
					int cur=i;
					for(int t=0;t<s.length();++t) {
						final int id=getIndex(s.charAt(t));
						if(nxt[cur][id]==-1) {
							cur=-1;
							break;
						}
						cur=nxt[cur][id]+1;
					}
					if(cur!=-1&&(g[i][j]==-1||g[i][j]>cur)) {
						g[i][j]=cur;
					}
				}
			}
		}
		int[][] dp=new int[n+1][m];
		for(int i=0;i<m;++i) {
			if(g[0][i]!=-1) {
				dp[g[0][i]][i]=1;
			}
		}
		for(int i=0;i<n;++i) {
			for(int j=0;j<m;++j) {
				if(dp[i][j]==0) {
					continue;
				}
				for(int k=0;k<m;++k) {
					if(g[i][k]!=-1) {
						dp[g[i][k]][k]=(dp[g[i][k]][k]+dp[i][j])%MOD;
					}
				}
			}
		}
		int result=0;
		for(int i=0;i<=n;++i) {
			for(int j=0;j<m;++j) {
				result=(result+dp[i][j])%MOD;
			}
		}
		return result;
	}

	List<List<String>> init(String[] codonTable) {
		Map<String,Integer> map=new HashMap<>();
		int id=0;
		List<List<String>> lists=new ArrayList<>();
		for(int i=0;i<codonTable.length;++i) {
			String s=codonTable[i];
			String key=s.substring(0,3);
			String value=s.substring(4);
			if(!map.containsKey(value)) {
				map.put(value,id++);
				lists.add(new ArrayList<>());
			}
			lists.get(map.get(value)).add(key);
		}
		return lists;
	}

	int getIndex(char c) {
		if(c=='A') {
			return 0;
		}
		else if(c=='C') {
			return 1;
		}
		else if(c=='T') {
			return 2;
		}
		return 3;
	}
}

  


code for problem3

import java.util.*;
import java.math.*;
import static java.lang.Math.*;

public class CompanyRestructuring {
	
	public int fewestDivisions(String[] hasManaged) {
		final int n=hasManaged.length;
		boolean[][] g=new boolean[n][n];
		for(int i=0;i<n;++i) {
			for(int j=0;j<n;++j) {
				if(hasManaged[i].charAt(j)=='Y') {
					g[i][j]=true;
				}
			}
		}
		for(int k=0;k<n;++k) {
			for(int i=0;i<n;++i) {
				for(int j=0;j<n;++j) {
					if(g[i][k]&&g[k][j]) {
						g[i][j]=true;
					}
				}
			}
		}
		boolean[] visited=new boolean[n];
		int result=0;
		for(int i=0;i<n;++i) {
			if(visited[i]) {
				continue;
			}
			List<Integer> left=new ArrayList<>();
			List<Integer> right=new ArrayList<>();
			for(int j=0;j<n;++j) {
				if(i==j||g[i][j]&&g[j][i]) {
					left.add(j);
					visited[j]=true;
				}
				else {
					right.add(j);
				}
			}
			if(right.size()==0) {
				result+=left.size();
				continue;
			}
			final int lNum=left.size();
			final int rNum=right.size();
			boolean[][] graph=new boolean[lNum][rNum];
			for(int ll=0;ll<lNum;++ll) {
				for(int rr=0;rr<rNum;++rr) {
					if(hasManaged[right.get(rr)].charAt(left.get(ll))=='Y') {
						graph[ll][rr]=true;
					}
				}
			}
			result+=lNum-maxMatch(graph);
		}
		return result;
	}

	int maxMatch(boolean[][] g) {
		final int n=g.length;
		final int m=g[0].length;
		int[] match=new int[m];
		Arrays.fill(match,-1);
		int result=0;
		for(int i=0;i<n;++i) {
			if(find(i,m,g,new boolean[m],match)) {
				++result;
			}
		}
		return result;
	}
	boolean find(int u,int m,boolean[][] g,boolean[] visited,int[] match) {
		for(int i=0;i<m;++i) {
			if(!visited[i]&&g[u][i]) {
				visited[i]=true;
				int t=match[i];
				match[i]=u;
				if(t==-1||find(t,m,g,visited,match)) {
					return true;
				}
				match[i]=t;
			}
		}
		return false;
	}
}

  

posted @ 2017-09-28 19:42  朝拜明天19891101  阅读(445)  评论(0编辑  收藏  举报