久久永久免费人妻精品下载,久久精品视频网,国产精品99久久久久久人

pku3167 Cow Patterns KMP好题

yzhw — Sun, 05 Jun 2011 01:35:00 GMT

题意�Q?br />K长度的模式串�Q�N长度的主�Ԍ��求多��个位置能够匚w��?br />�q�里模式串以及主串里的值均为rank�Q�即15567�?6678认�ؓ是匹配的�?br />
解法�Q?br />对于�q�种rank�c�，判断是否相等��是判断在此��C��前大于这个数字的数字个数、等于这个数字的数字个数以及��于�q�个数字的数字个数是否均相等。接下来的过�E�就�c�M��KMP的匹配了�?br />另外�Q�KMP nxt数组本质是最长后�~�=最长前�~��Q?span style="color: red;">同时要求后缀不能为整个字�W�串。好久不写KMP�Q�竟然犯了个�p�涂的事情，nxt数组里算出来竟然�?,2,3,4.。。无语�?br />
代码�Q?br />

1 Source Code
2 Problem: 3167        User: yzhw
3 Memory: 24496K        Time: 813MS
4 Language: GCC        Result: Accepted
5
6     Source Code
7
8     # include <stdio.h>
9     # include <string.h>
10     # define N 100005
11     # define M 30
12     int n,k,s;
13     int cal1[N][M],cal2[N][M],nxt[N],d1[N],d2[N];
14     int equal(int *e1,int *s1,int *e2,int *s2,int val1,int val2)
15     {
16         if(e1[val1-1]-s1[val1-1]==e2[val2-1]-s2[val2-1]&&
17            (e1[val1]-e1[val1-1])-(s1[val1]-s1[val1-1])==(e2[val2]-e2[val2-1])-(s2[val2]-s2[val2-1])&&
18            (e1[s]-e1[val1])-(s1[s]-s1[val1])==(e2[s]-e2[val2])-(s2[s]-s2[val2]))
19         return 1;
20         else return 0;
21     }
22     int main()
23     {
24         int i,t,j,p,total=0;
25         //freopen("cpattern.9.in","r",stdin);
26         //freopen("ans.txt","w",stdout);
27
28         scanf("%d%d%d",&n,&k,&s);
29         memset(cal1,0,sizeof(cal1));
30         memset(cal2,0,sizeof(cal2));
31         for(i=1;i<=n;i++)
32         {
33             scanf("%d",&t);
34             d2[i]=t;
35             //cal2[i][t]++;
36             for(j=1;j<=s;j++)
37                 cal2[i][j]+=cal2[i-1][j];
38             for(j=t;j<=s;j++)
39                 cal2[i][j]++;
40
41         }
42         for(i=1;i<=k;i++)
43         {
44             scanf("%d",&t);
45             d1[i]=t;
46             //cal1[i][t]++;
47             for(j=1;j<=s;j++)
48                 cal1[i][j]+=cal1[i-1][j];
49             for(j=t;j<=s;j++)
50                 cal1[i][j]++;
51         }
52         nxt[1]=0;
53         for(i=2;i<=k;i++)
54         {
55             p=nxt[i-1];
56             while(p&&!equal(cal1[p+1],cal1[0],cal1[i],cal1[i-p-1],d1[p+1],d1[i]))
57                 p=nxt[p];
58             if(equal(cal1[p+1],cal1[0],cal1[i],cal1[i-p-1],d1[p+1],d1[i])) nxt[i]=p+1;
59             else nxt[i]=0;
60         }
61         p=0;
62         for(i=1;i<=n;i++)
63         {
64             while(p&&!equal(cal1[p+1],cal1[0],cal2[i],cal2[i-p-1],d1[p+1],d2[i])) p=nxt[p];
65             if(equal(cal1[p+1],cal1[0],cal2[i],cal2[i-p-1],d1[p+1],d2[i])) p++;
66             if(p==k)
67             {
68                 total++;
69                 p=nxt[p];
70             }
71         }
72         printf("%d\n",total);
73         p=0;
74         for(i=1;i<=n;i++)
75         {
76             while(p&&!equal(cal1[p+1],cal1[0],cal2[i],cal2[i-p-1],d1[p+1],d2[i])) p=nxt[p];
77             if(equal(cal1[p+1],cal1[0],cal2[i],cal2[i-p-1],d1[p+1],d2[i])) p++;
78             if(p==k)
79             {
80                 printf("%d\n",i-k+1);
81                 p=nxt[p];
82             }
83         }
84         return 0;
85     }
86
87

yzhw 2011-06-05 09:35 发表评论

pku1245-1251 Mid-Central USA 2002 比赛�ȝ��

yzhw — Thu, 20 Jan 2011 18:21:00 GMT

摘要: �W�二�ơ做比赛�Q�有了viaxl的加入，有点动力 SolvedIDTitle Ratio(AC/att) Yes Problem A Pr... 阅读全文

yzhw 2011-01-21 02:21 发表评论

pku 1816 wild words trie树的�zȝ��

yzhw — Thu, 13 Jan 2011 10:05:00 GMT

题意�Q?br>模板包括3�U�字�W�：
a-z�Q�标准字�W?br>通配�W?�Q�代表�Q意一个字�W?br>*�Q�代�?个或多个字符

先给出n个模板串�Q�m个字�W�串�Q�问每个字符串分别匹配那几个模板�?br>n<=1e6,m<=1e2�Q�strlen<20

解法�Q?br>用trie树构建模板串�Q�然后再dfs来匹�?br>应�ؓ模板串中有重复串�Q�一个简单的�Ҏ��是，开辟一个链表空��_��链接每个模板串的末节�?br>然后判断每个字符串时DFS一�ơ，��能到达的尾部节点做个标讎ͼ�然后�Ҏ��之前记录的链表扫一遍就知道到达�q�哪些模板串的根节点�Q�匹配了哪些模板�Ԍ��
DFS的时候要注意下几点：
1、状态设定：(节点�~�号、当前字�W�串位置)
2、如果当前节点中存在'?"的�{�U��\径，则�{�U�过�?br>3、如果当前节点存�?*'的�{�U��\径，枚�D*匚w��的长度，然后转移
��d��杂度1e7左右。。。不知道有什么好的方法，我的�E�序C++跑了400MS。。。不��快的。�?br>
另外�Q�这题徏议用动态内存分配，我开始开�?00W的buffer�Q�竟然MLE。。。�?br>代码�Q?br>

1# include <cstdio>
2# include <cstring>
3# include <cstdlib>
4
5using namespace std;
6struct node
7{
8    node *nxt[28];
9    bool in;
10    node()
11    {
12          memset(nxt,NULL,sizeof(nxt));
13          in=false;
14    }
15}*ans[100005],head;
16char map[255];
17int c=1;
18node* insert(char *str)
19{
20   node *p=&head;
21   for(int i=0;str[i]!='\0';i++)
22   {
23      if(p->nxt[map[str[i]]]==NULL)
24         p->nxt[map[str[i]]]=new node();
25      p=p->nxt[map[str[i]]];
26   }
27   return p;
28}
29void match(node *p,char *str)
30{
31     if(!p) return;
32     if(*str=='\0')
33     {
34        p->in=true;
35        match(p->nxt[27],str);
36     }
37     else
38     {
39         match(p->nxt[map[*str]],str+1);
40         match(p->nxt[26],str+1);
41         for(int i=0;i<=strlen(str);i++)
42           match(p->nxt[27],str+i);
43     }
44}
45int main()
46{
47    int n,m;
48    for(int i='a';i<='z';i++)
49       map[i]=i-'a';
50    map['?']=26;
51    map['*']=27;
52    scanf("%d%d",&n,&m);
53    for(int i=0;i<n;i++)
54    {
55       char str[10];
56       scanf("%s",str);
57       ans[i]=insert(str);
58    }
59    for(int i=0;i<m;i++)
60    {
61        for(int j=0;j<n;j++) ans[j]->in=false;
62        char str[128];
63        scanf("%s",str);
64        match(&head,str);
65        bool flag=false;
66        for(int j=0;j<n;j++)
67          if(ans[j]->in)
68           {
69              flag=true;
70              printf("%d ",j);
71           }
72        if(flag) printf("\n");
73        else printf("Not match\n");
74    }
75  // system("pause");
76    return 0;
77}
78

yzhw 2011-01-13 18:05 发表评论

yzhw — Thu, 13 Jan 2011 01:57:00 GMT

题意�Q?br>问不包含n个子串的长度为m的字�W�串的构造个数�?br>
解法�Q?br>构造trie图，然后DP求长度�ؓm的合法串个数
以前高精度都靠java�Q�这�ơ手写，各种错误。。。唉。�?br>
代码�Q?br>

  1 # include <iostream>
  2 # include <cstring>
  3 using namespace std;
  4 struct BigInteger
  5 {
  6     int bit[100];
  7     bool init;
  8     BigInteger()
  9     {
10         memset(bit,0,sizeof(bit));
11         init=true;
12     }
13     BigInteger operator+(const BigInteger &pos)
14     {
15         BigInteger res;
16         res.init=init;
17         for(int i=0;i<99;i++)
18         {
19             res.bit[i]+=bit[i]+pos.bit[i];
20             res.bit[i+1]+=res.bit[i]/10;
21             res.bit[i]%=10;
22         }
23         return res;
24     }
25     void print()
26     {
27         int i;
28         for(i=99;i>0&&!bit[i];i--);
29         for(int j=i;j>=0;j--)
30             cout<<bit[j];
31         cout<<endl;
32     }
33 };
34 struct node
35 {
36     node *nxt[51],*pre;
37     bool end;
38     void clear()
39     {
40         memset(nxt,NULL,sizeof(nxt));
41         end=false;
42         pre=NULL;
43     }
44     node()
45     {
46         clear();
47     }
48 }buffer[200];
49 int map[255];
50 int c=1,n,m,num;
51 void insert(char *str)
52 {
53     node *p=buffer;
54     for(int i=0;str[i]!='\0';i++)
55     {
56         if(!(p->nxt[map[str[i]]]))
57                 p->nxt[map[str[i]]]=&buffer[c++];
58         p=p->nxt[map[str[i]]];
59     }
60     p->end=true;
61 }
62 void make_per()
63 {
64     int s=-1,e=-1;
65     node *q[200];
66     node *p=buffer;
67     for(int i=0;i<n;i++)
68         if(p->nxt[i])
69         {
70             p->nxt[i]->pre=p;
71             q[++e]=p->nxt[i];
72         }
73         else
74             p->nxt[i]=p;
75     while(s!=e)
76     {
77         p=q[++s];
78         for(int i=0;i<n;i++)
79         {
80             node *pre=p->pre;
81             while(pre->nxt[i]==NULL) pre=pre->pre;
82             if(p->nxt[i])
83             {
84                 p->nxt[i]->pre=pre->nxt[i];
85                 p->nxt[i]->end=(p->nxt[i]->pre->end||p->nxt[i]->end);
86                 q[++e]=p->nxt[i];
87             }
88             else
89                 p->nxt[i]=pre->nxt[i];
90         }
91     }
92 }
93 BigInteger dp[200][55],zero,one;
94 BigInteger solve(int s,node *p)
95 {
96     if(p->end) return zero;
97     else if(s==m) return one;
98     else if(!(dp[p-buffer][s].init)) return dp[p-buffer][s];
99     else
100     {
101         dp[p-buffer][s].init=false;
102         for(int i=0;i<n;i++)
103         {
104             dp[p-buffer][s]=dp[p-buffer][s]+solve(s+1,p->nxt[i]);
105         }
106         return dp[p-buffer][s];
107     }
108
109 }
110 int main()
111 {
112     cin>>n>>m>>num;
113     char str[128];
114     cin>>str;
115     int tc=0;
116     for(int i=0;str[i]!='\0';i++)
117         map[str[i]]=tc++;
118     while(num--)
119     {
120         cin>>str;
121         insert(str);
122     }
123     make_per();
124     one.bit[0]=1;;
125     solve(0,buffer);
126     dp[0][0].print();
127     return 0;
128
129 }

yzhw 2011-01-13 09:57 发表评论

pku3691 DNA repair 自动�?DP

yzhw — Wed, 12 Jan 2011 05:21:00 GMT

摘要: 题意�Q�给��Z��个DNA序列�Q�和若干病毒DNA片段�Q�问最��修改DNA的次��C��得DNA序列中不包含病毒解法�Q�首先构造DFA�Q�再�ơ犯了和��州现场赛一��L��错误�Q�没�?�l�束标记状态的转移。。。。WA了N久）�Q�然后就是在状态机上的DP了��odp[i]][j]为当前在自动机第i个节点上时suffix(j)成�ؓ合法子串需要修改的最��次数。下面状态�{�U�d��单了。。（�q�种状态�{�U�通过方程和难描述�Q�大家看�E�序吧，我写的很�?.. 阅读全文

yzhw 2011-01-12 13:21 发表评论

yzhw — Thu, 06 Jan 2011 15:26:00 GMT

摘要: 题目很罗嗦，说到底干三�g事：1、统计句子的个数�Q�以".:;!?"为分隔符�Q�，�q�个只要�l�计分隔�W�的个数卛_��2、统计单词的个数�Q�这个需要注意一点，不能仅仅以空格来spilt字符�Ԍ��然后�l�计个数的方法。。只有除��L��点符号后非空的patten才能��一个word3、统计音节的个数�? -es, -ed and -e (except -le) endings are ignored, �Q�末��以es... 阅读全文

yzhw 2011-01-06 23:26 发表评论

pku 2945 Find the Clones 字符串hash

yzhw — Fri, 10 Dec 2010 08:55:00 GMT

题意�Q?br>�l�出n个字�W�串�Q�统计出�?�?�?.n�ơ的字符串各有多��种
解法�Q?br>字符串hash�Q�我用的强大无比的Java的HashMap�Q�非常给力～
代码�Q?br>

1 import java.io.*;
2 import java.util.*;
3 public class Main {
4
5     /**
6      * @param args
7      */
8     public static void main(String[] args) throws IOException{
9         BufferedReader in=new BufferedReader(new InputStreamReader(System.in));
10         HashMap<String,Integer> refer=new HashMap<String,Integer>();
11         int cal[]=new int[20001];
12         while(true)
13         {
14             String tmp[]=in.readLine().split(" ");
15             if(Integer.parseInt(tmp[0])==0&&Integer.parseInt(tmp[1])==0) break;
16             int n=Integer.parseInt(tmp[0]);
17             refer.clear();
18             for(int i=0;i<n;i++)
19             {
20                 String str=in.readLine();
21                 if(refer.containsKey(str)) refer.put(str,refer.get(str)+1);
22                 else refer.put(str, 1);
23             }
24             Arrays.fill(cal,0);
25             for(int i:refer.values())
26               cal[i]++;
27             for(int i=1;i<=n;i++)
28                 System.out.println(cal[i]);
29
30         }
31
32     }
33
34 }

yzhw 2010-12-10 16:55 发表评论

yzhw — Thu, 09 Dec 2010 13:09:00 GMT

题意�Q?br>�l�出一堆字�W�串�Q�设X��Z��个字�W�串�Q�其本��n或者其反�{传在每个字符串中均出现。问x的最镉K��度�?br>解法�Q?br>�׃��数据��，100个字�W�串�Q�最大长度�ؓ100�Q�我��q��了二�?暴力枚�D验证。�?br>二分字符串长度，然后枚�D其中一个字�W�串的�v始位�|�s�Q�验证以s��v始，长度为len的字�W�串�Q�或其反转串�Q�是否在所有字�W�串中都出现�Q�这步复杂度o(n3))
加上二分的复杂度�Q��d��杂度为o(n³logn)
正解应该是用后缀数组+二分。将所有原串和原串的反串连接�v来计��后�~�数组�Q�二分子串长度，然后用height数组分类来验证。。复杂度o(n²logn)

代码�Q?br>

1# include <cstdio>
2# include <cstring>
3using namespace std;
4char str[101][105];
5int n;
6void strreverse(char *str)
7{
8   int len=strlen(str);
9   for(int i=0;i<len/2;i++)
10   {
11      char tmp=str[i];
12      str[i]=str[len-1-i];
13      str[len-1-i]=tmp;
14   }
15}
16bool chk(int len,int pos)
17{
18   for(int s=0;s+len<=strlen(str[pos]);s++)
19   {
20      char tmp=str[pos][s+len];
21      bool flag=true;
22      str[pos][s+len]='\0';
23      for(int j=0;j<n;j++)
24      {
25         if(strstr(str[j],str[pos]+s)) continue;
26         else
27         {
28             strreverse(str[pos]+s);
29             char *p=strstr(str[j],str[pos]+s);
30             strreverse(str[pos]+s);
31             if(!p)
32             {
33               flag=false;
34               break;
35             }
36         }
37      }
38      str[pos][s+len]=tmp;
39      if(flag) return true;
40
41   }
42   return false;
43}
44int main()
45{
46   int test;
47   scanf("%d",&test);
48   while(test--)
49   {
50      int s=0,e=0xfffffff,pos;
51      scanf("%d",&n);
52      for(int i=0;i<n;i++)
53      {
54        scanf("%s",str[i]);
55        if(strlen(str[i])<e) e=strlen(str[i]),pos=i;
56      }
57      while(s<=e)
58      {
59         int mid=(s+e)>>1;
60         if(chk(mid,pos)) s=mid+1;
61         else e=mid-1;
62      }
63      printf("%d\n",s-1);
64   }
65   return 0;
66}
67
68

yzhw 2010-12-09 21:09 发表评论

yzhw — Mon, 06 Dec 2010 17:36:00 GMT

摘要: 题意�Q�给��Z��些病毒的特征�Ԍ��如果一个程序（或者将�E�序反�{�Q�中出现了某个病毒的特征�Ԍ��则该�E�序被这个病毒感染了。给��q�病毒串�Q�一个程序串�Q�问改程序被多少�U�病毒感染了�Q�解法：比赛的时候模板有bug�Q�WA到死�Q�竟然这个用了数月的模板之前�q�神奇的通过了N道自动机的题目，不可思议。。一道非常裸的自动机�Q�将�E�序串正反匹配一遍既得答案。自动机节点�l�构如下�Q?Code highlighting produce... 阅读全文

yzhw 2010-12-07 01:36 发表评论

yzhw — Sat, 27 Nov 2010 13:06:00 GMT

题意�Q?br>�l�出一个字�W�串�Q�求其所有��@环同构串中字典序最大的串以及最��的丌Ӏ��ƈ且计��这两个串在所有��@环同构串中出现的�ơ数

解法�Q?br>�W�一问，用经典的求串的最��表�C�的��法��可以了�?br>�W�二问，利用KMP��法前缀数组的性质�Q�即最大后�~��{�于最长前�~�的位�|�。从最��表�C�（最大表�C�）的位�|�开始计��next函数�Q�将start+len-1位置作合法标志，计算nxt数组的时候顺便计��标讎ͼ�与AC自动机方法相同）�Q�然后统计标��C��数即可～

代码�Q�（GCC�Q?br>

1 # include <string.h>
2 # include <stdio.h>
3 # include <stdbool.h>
4 char str[2050000];
5 int pre[2050000];
6 bool end[2050000];
7 int maxpos(int len)
8 {
9     int p1=0,p2=1,l=0,i;
10     while(p1<len&&p2<len&&l<len)
11     {
12         int cmp=str[p1+l]-str[p2+l];
13         if(!cmp)
14              l++;
15         else
16         {
17             if(cmp<0) p1+=l+1;
18             else p2+=l+1;
19             l=0;
20             p2+=(p2==p1);
21         }
22     }
23     return p1<p2?p1:p2;
24 }
25 int minpos(int len)
26 {
27     int p1=0,p2=1,l=0,i;
28     while(p1<len&&p2<len&&l<len)
29     {
30         int cmp=str[p1+l]-str[p2+l];
31         if(!cmp)
32              l++;
33         else
34         {
35             if(cmp>0) p1+=l+1;
36             else p2+=l+1;
37             l=0;
38             p2+=(p2==p1);
39         }
40     }
41     return p1<p2?p1:p2;
42 }
43 int gettimes(int start,int len)
44 {
45     int p,res=0;
46     memset(end,0,sizeof(end));
47     end[start+len-1]=1;
48     pre[start]=start-1;
49     for(p=start+1;p<2*len-1;p++)
50     {
51         pre[p]=pre[p-1];
52         while(pre[p]!=start-1&&str[pre[p]+1]!=str[p]) pre[p]=pre[pre[p]];
53         if(str[pre[p]+1]==str[p]) pre[p]++;
54         if(pre[p]!=start-1&&!end[p]) end[p]=end[pre[p]];
55         res+=end[p];
56     }
57     return res;
58 }
59 int main()
60 {
61     while(gets(str))
62     {
63         int len=strlen(str),i;
64         for(i=len;i<2*len-1;i++) str[i]=str[i-len];
65         str[2*len-1]='\0';
66         int p1=minpos(len);
67         int t1=gettimes(p1,len);
68         printf("%d %d ",p1+1,t1);
69         p1=maxpos(len);
70         t1=gettimes(p1,len);
71         printf("%d %d\n",p1+1,t1);
72     }
73     return 0;
74 }
75

yzhw 2010-11-27 21:06 发表评论

yzhw — Sat, 27 Nov 2010 11:59:00 GMT

把一个长为len的字�W�串围成一个圈�Q�然后以��L��一个字�W�作��v点，都会产生一个新的长为len的字�W�串�Q�字�W�串的最��表�C�就是所有新字符串中字典序最��的那个�?br>下面�q�个函数��是解决�q�个问题的，�q�回��gؓ字典序最��的串的在原串中的�v始位�|��?br>

1 int MinimumRepresentation(char *s,int l)//串s[0~l-1]的最��表�C�Z��|?/span>
2 {
3     int i = 0, j = 1, k = 0,t;
4     while (i < l && j < l && k < l)//找不到比它还��的或�?nbsp;完全匚w��
5     {
6         t = s[(i+k)%l] - s[(j+k)%l];
7         //if (s[(i+k) >= l ? i+k-l : i+k] == s[(j+k) >= l ? j+k-l : j+k])
8         if (t == 0)
9             k++;//相等的话,��长度加1
10         else
11         {
12             if (t > 0)//大于的话,s[i]为首的肯定不是最��表�C?最大表�C�就�?lt;
13                 i += k + 1;
14             else
15                 j += k + 1;
16             if (i == j)
17                 j++;
18             k = 0;
19         }
20     }
21     return min(i,j);
22 }

基本��x��是两个位置的字�W�比较，如果s[i+k] > s[j+k]那么i到i+k位置都不是最��表�C�的位置�Q�所以i直接跳k+1步，反之j直接跳k+1步�?br>本题代码�Q?br>

1 import java.io.*;
2 public class Main {
3     static int minpos(String str)
4     {
5         int p1=0,p2=1,len=0;
6         while(p1<str.length()&&p2<str.length()&&len<str.length())
7         {
8             int res=str.charAt((p1+len)%str.length())-str.charAt((p2+len)%str.length());
9             if(res==0)
10                 len++;
11             else
12             {
13                 if(res>0) p1+=len+1;//如果是最大表�C�，则p2+=len+1�Q�下面亦�?/span>
14                 else p2+=len+1;
15                 len=0;
16                 p2=p2+(p1==p2?1:0);
17             }
18         }
19         return Math.min(p1, p2)+1;
20     }
21     public static void main(String[] args) throws IOException{
22         BufferedReader in=new BufferedReader(new InputStreamReader(System.in));
23         int test=Integer.parseInt(in.readLine());
24         while((test--)!=0)
25             System.out.println(minpos(in.readLine()));
26     }
27
28 }
29

yzhw 2010-11-27 19:59 发表评论