- 浏览: 380150 次
- 性别:
- 来自: 北京
文章分类
- 全部博客 (229)
- java编程 (4)
- java实用程序 (2)
- 算法设计 (34)
- 数据库 (8)
- ACM模板 (12)
- 技术术语 (1)
- java_web (3)
- php (22)
- eclipse (3)
- linux (25)
- linux命令使用心得 (3)
- web服务器 (8)
- IT知识 (2)
- 前端技术 (17)
- 开源软件 (5)
- vim (3)
- linux多线程 (9)
- web开发经验 (3)
- lua (5)
- linux编程 (3)
- smarty (1)
- mysql (4)
- Hive (2)
- 数据挖掘 (9)
- python (2)
- 生活 (1)
- C++ (2)
- 计算机 (1)
- objective-c (11)
- css (2)
- 游戏 (1)
- Mac (1)
最新评论
-
lr544463316:
我的怎么不行呀.....
Mysql Access denied for user ''@'localhost' to database 的一种解决方法 -
babaoqi:
使用时需要注意group_concat函数返回值的最大长度=g ...
mysql中的group_concat函数 -
代码能力弱成渣:
可以帮我看下我的代码么?我自己写的sam,也有ac过题的,但是 ...
求两个字符串的最长公共连续子序列(SAM实现) -
atgoingguoat:
有1000个?不过还是收藏下。
jquery常用的插件1000收集(转载)
/* * AC自动机,先对资源串和病毒串构成的字符串集合建立AC自动机,然后在trie树上做BFS求出在安全图上每个资源串 * 到其他资源的最短路径,最后做一遍状态压缩dp即可 */ #include <cstdio> #include <cstring> #include <map> using namespace std; const int N=10*1005+50005; //节点个数的最大值 const int S=2; //不同的字符 个数 const int NUM=10; const int INF=0X7FFFFFFF; struct node{ node *sons[S], *fail; int d; //d为-1表示既不是病毒也不是资源的串,为-2表示病毒串,大于等于0表示资源串,并表示为资源的id }nodes[N], *root; int cnt; //cnt是树中节点的 个数 node *que[N]; int n, m; node *ns[NUM]; int dis[NUM][NUM], lens[NUM]; int mark[N], d[N]; int dp[NUM][1<<NUM]; void clear(){ cnt=0; root=NULL; } node* newNode(){ node* ans=&nodes[cnt++]; memset(ans->sons, 0, S*sizeof(node*)); ans->fail=NULL; ans->d=-1; return ans; } int hash(char ch){ //字符的哈希函数,根据不同的需要而定 return ch-'0'; } int id(node*& t){ //计算该节点的id return t-nodes; } //j表示当前关键字的标号 node* insert(node*& root, char* str, int d){ node* t=root; int i, k; for(i=0; str[i]; i++){ if(t->sons[k=hash(str[i])]==NULL){ t->sons[k] = newNode(); } t=t->sons[k]; } t->d=d; return t; } void getFail(node*& root){ int l, r, i; node *t; l=r=0; root->fail=root; //这样可以保证每个节点的fail指针都是非空的 for(que[r++]=root; l!=r; ){ t=que[l++]; for(i=0; i<S; i++){ if(t->sons[i]){ que[r++]=t->sons[i]; if(t==root){ t->sons[i]->fail=t; }else{ t->sons[i]->fail=t->fail->sons[i]; } }else{ //增设虚拟节点 if(t==root) t->sons[i]=t; else t->sons[i]=t->fail->sons[i]; } } } } bool input(){ scanf("%d%d", &n, &m); if(n==0 && m==0) return false; char str[1005]; int i; clear(); root=newNode(); for(i=0; i<n; i++){ scanf("%s", str); lens[i]=strlen(str); ns[i]=insert(root, str, i); } for(i=0; i<m; i++){ scanf("%s", str); insert(root, str, -2); } return true; } void BFS(int s){ int l, r, j; node *u, *v; l=r=0; for(que[r++]=&nodes[s], mark[s]=s, d[s]=0; l!=r; ){ u=que[l++]; for(j=0; j<S; j++){ v=nodes[id(u)].sons[j]; if(v->d==-2) continue; if(mark[id(v)]!=s){ mark[id(v)]=s; d[id(v)]=d[id(u)]+1; que[r++]=v; } } } } void solve(){ getFail(root); int i, j, k, tmp; for(i=0; i<n; i++){ for(j=0; j<n; j++){ dis[i][j]=(i==j ? 0: INF); } } for(i=0; i<cnt; i++){ mark[i]=-1; } for(i=0; i<n; i++){ BFS(id(ns[i])); for(j=0; j<n; j++){ if(mark[id(ns[j])]!=id(ns[i])) continue; if(dis[i][j]>d[id(ns[j])]){ dis[i][j]=d[id(ns[j])]; } } } int sum=1<<n; for(i=0; i<n; i++){ for(j=0; j<sum; j++){ dp[i][j]=INF; } dp[i][1<<i]=lens[i]; } for(j=1; j<sum; j++){ for(i=0; i<n; i++){ if(!(j&(1<<i)) || dp[i][j]==INF) continue; for(k=0; k<n; k++){ if(j&(1<<k) || dis[i][k]==INF) continue; tmp=dp[i][j]+dis[i][k]; if(dp[k][j|(1<<k)] > tmp){ dp[k][j|(1<<k)] = tmp; } } } } int ans=INF; for(i=0; i<n; i++){ if(ans>dp[i][sum-1]){ ans=dp[i][sum-1]; } } printf("%d\n", ans); } int main(){ //freopen("in.txt", "r", stdin); while(input()) solve(); return 0; }
发表评论
-
升序数组中求一个key出现的次数
2013-01-09 23:08 1068算法思路: 在排好序的数组,相同的数字是排列在一起的,所以只需 ... -
判断单链表是否有环
2013-01-08 19:07 888算法思路: 指针p1和p2的起始值均为链表的表头,指针p1每次 ... -
hdu3684
2011-11-15 20:11 894/* 刚开始打了个记录上下左右四个点的,一直tle。 ... -
hdu3686
2011-11-14 20:43 1006/* 无向图边的双连通分量,在同一个连通分量里的边之间 ... -
poj3968
2011-11-14 04:45 1379source: http://poj.org/problem ... -
uva2819
2011-11-13 02:20 863source: http://livearchive.onli ... -
manacher算法
2011-11-11 00:06 2412const int LEN=110005; const ... -
hdu4118
2011-11-09 21:53 1163枚举每条边最多被经过的次数即可 #include ... -
hdu4115
2011-11-09 16:27 998source: http://acm.hdu.edu.cn/ ... -
uva(Transitive Closure)
2011-11-08 14:45 884source: http://livearchive.onli ... -
zoj3500
2011-11-07 17:41 926求两个球的体积交或者并 #include <cs ... -
zoj3545
2011-11-04 18:18 842/* AC自动机 相当暴力的 解法: mark[i ... -
zoj3228
2011-11-04 16:12 928/* * AC自动机,每个节点 添加一个d表示节点代 ... -
poj3691(DNA Repair)
2011-11-04 13:18 1454/* AC自动机,增设虚拟节点,求长度为n的字符串中包 ... -
hdu2825
2011-11-04 11:53 961/* AC自动机,增设虚拟节点,求长度为n的字符 ... -
hdu4095
2011-11-03 13:19 991/* 第一步,构建BST,用第一个数作为bst的 ... -
zoj3540
2011-11-02 21:33 894/* 其实就是把总共的 放置次数减去不能放置的那些就行 ... -
poj1741(树的分治,基于边的 分治)
2011-11-02 20:25 3326/* 树基于边的分治算法,计算树中距离小于等于k的点 ... -
hdu2939
2011-10-29 18:36 823source: http://acm.hdu.edu.cn/s ... -
hdu3982
2011-10-29 10:20 909//半平面交+多边形和圆的交的面积 #include ...
相关推荐
ZOJ解题报告ZOJ解题报告ZOJ解题报告ZOJ解题报告
zoj题目简单归类zoj题目简单归类zoj题目简单归类
acm中zoj1002的可运行C++程序
包含了zoj700多道题目的源代码,在做题时可以参考
Problem Arrangement zoj 3777
ZOJ题目答案源码
学习ACM程序设计的朋友一定要看,这是训练必备的POJ ZOJ题目分类及解题思路
一个非常非常非常非常实用的zoj结题代码
浙大ZOJ题目分类,可以让你更方便快速锁定那你想要联系的题目,是自己快速提高·
zoj 1003 c语言的,要写这么多描述吗。。
本代码是zoj上AC的1951的代码,把双重循环简化为O(n),不过素数判断的改进还不够
ZOJ1805代码
ZOJ题解集合-截至2835。共1244个文件,C/C++,有重复
zoj1027解题指南和代码,还不错,是学校培训给的。
zoj吐血制作,希望大家喜欢
zoj4041正确题解源代码,以及运行程序
zoj 题库 详细解答 解题代码 acm
大学ACM竞赛,ZOJ 1733 运用递归(优化)的方法。ac的代码。
能AC 通过的c++代码,包括zoj1002,1091,1789
ZOJ完全解题报告,喜欢ACM的同学,欢迎下载