C++博客-250的oi之路-随笔分类-oihttp://www.cppblog.com/wwy250/category/9749.htmlzh-cnTue, 30 Jun 2009 11:26:26 GMTTue, 30 Jun 2009 11:26:26 GMT60这里的最后一篇http://www.cppblog.com/wwy250/archive/2009/06/29/88765.html250250Mon, 29 Jun 2009 05:41:00 GMThttp://www.cppblog.com/wwy250/archive/2009/06/29/88765.htmlhttp://www.cppblog.com/wwy250/comments/88765.htmlhttp://www.cppblog.com/wwy250/archive/2009/06/29/88765.html#Feedback1http://www.cppblog.com/wwy250/comments/commentRss/88765.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/88765.html
现在说回我吧
关注我的blog的朋友可能发现 我这里已经好久没更新了
这还要从CTSC&APIO说起
话说那两次比赛考得很伪 究其原因 恐怕还是学习方法的问题
在这方面JL的oier有着得天独厚的劣势
所以 我经常说 作为一个弱省的oier 不要怕走弯路 这一定是必经的过程
在那次回来之后 我就再也没有写blog 原因大概就是对自己今年Au 失去了像CTSC之前那么强的信心 甚至怕省选落选 每天就做做各省省选题 心想Ag就不错了
后来省选又结束了 我以第二名的成绩进队 其实也是情理之中的
之后就是做做NOI 好多题都不能AC 但是要是把力所能及的点都做上 还不Au线高不少 这让我又有了一点点幻想 但我深知基本不会有人在NOI吧力所能及的点都答上的
这两天去参加中考 其实就是考这玩 我上高中一定是要靠OI的成绩的
参加中考见到了许多同学 半年没见到这么多同学了 还是这些同学好
在那个圈子里 没人会认为你是神牛 没人会认为你可能进国家队甚至IOI Au 也没人把你当小朋友 即使相信也不会有不断的orz 不会吝啬批评和鄙视甚至是鸡蛋里挑骨头
虽然少了写前辈与众心肠的话 但是却显得特别的真实
由于我们班是最BT的班 所以我的同学们也就都选择了最BT的高中 合适JL省太不重视OI了 所以我一定要去的Au才能保证高中也和他们在一起
为了同学们 为了不辜负当年手把手教我oi的lpq、给我讲了无数题的zys、在我孤独是陪伴我的wh 以及好多好多人
就让我再修炼28天 在NOI取得Au吧

这将是我这个blog的最后一篇blog 这次NOI之前先不写blog了 等这次NOI之后 重新开个blog
另外看了几个我们同学的blog发现原来blog不止是写解题报告的 尤其是这个写的太NB了

不Au就要离开他们了 5555555555555555555
不过即使不在一起上学 哪个学校的oi老师都会欢迎我的 想见还是能见到 而且可以经常见






250 2009-06-29 13:41 发表评论
]]>CTSC&APIO 09http://www.cppblog.com/wwy250/archive/2009/05/03/81785.html250250Sun, 03 May 2009 12:14:00 GMThttp://www.cppblog.com/wwy250/archive/2009/05/03/81785.htmlhttp://www.cppblog.com/wwy250/comments/81785.htmlhttp://www.cppblog.com/wwy250/archive/2009/05/03/81785.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/81785.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/81785.html虽然我很重视这两个比赛(尤其是前者) 但我并没有特意为他们准备 这是按照至NOI09的计划进行着
不敢奢望能取得什么太好的成绩 但总不能让我自己掏比赛费可车票(我们学校规定 Ag或以上成绩可以报销)
下面将写下我的比赛成绩由于事情还没有发生 所以现在还不能写 期望回来的时候不要让我不好意思写


250 2009-05-03 20:14 发表评论
]]>NOI 09 还有3个月http://www.cppblog.com/wwy250/archive/2009/04/27/81274.html250250Mon, 27 Apr 2009 14:53:00 GMThttp://www.cppblog.com/wwy250/archive/2009/04/27/81274.htmlhttp://www.cppblog.com/wwy250/comments/81274.htmlhttp://www.cppblog.com/wwy250/archive/2009/04/27/81274.html#Feedback3http://www.cppblog.com/wwy250/comments/commentRss/81274.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/81274.html在于何林的聊天中不止一次表示3个月足以变得很强
他还鼓励我今年进集训队 叫我自信
现在足以改变一切的3个月终于到来了
希望借何林神牛的吉言 用3个月 原了我的Au梦


250 2009-04-27 22:53 发表评论
]]>后缀数组http://www.cppblog.com/wwy250/archive/2009/04/14/79874.html250250Tue, 14 Apr 2009 04:29:00 GMThttp://www.cppblog.com/wwy250/archive/2009/04/14/79874.htmlhttp://www.cppblog.com/wwy250/comments/79874.htmlhttp://www.cppblog.com/wwy250/archive/2009/04/14/79874.html#Feedback2http://www.cppblog.com/wwy250/comments/commentRss/79874.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/79874.html09年罗XX的论文比那篇容易理解的多
后缀数组果然是一个很强大的东西 有关字符串的问题 只要学好自动机和后缀数组基本都能解决了
与其说是后缀数组强大 不如说height数组强大 后缀数组的作用也就是方便求解height数组
使用height数组 大致有如下你个常用方法
1、分组:将height值大于等于k的分置一组 使得同组内最长公共前缀>=k
2、二分:根据具体情况 一般都是二分答案
3、连接：将所有涉及到的字符串连在一起处理
还有一个神奇的性质:当循环节长度确定时 保证覆盖是s[0],s[l],s[l*2]～～中的某连续两个


250 2009-04-14 12:29 发表评论
]]>病毒的DNAhttp://www.cppblog.com/wwy250/archive/2009/04/13/79835.html250250Mon, 13 Apr 2009 15:33:00 GMThttp://www.cppblog.com/wwy250/archive/2009/04/13/79835.htmlhttp://www.cppblog.com/wwy250/comments/79835.htmlhttp://www.cppblog.com/wwy250/archive/2009/04/13/79835.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/79835.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/79835.html贴一下代码：
#include<iostream>
#define Rank(x) rank[(x)<?(n+1)]

using namespace std;

char a[100000],c[5000],b[5000],d[100000];
int sa[5001],height[5001],rank[5001],n,sum[5001],h[5001],_sa[5001],_rank[5001],m,ans[20][5001],tr[5001],tl[5001],maxk,maxi,lft[100001],rght[100001],leftt[100001],rightt[100001];

struct node
{
    int th;
    char s;
}r[10000];

int cmps(const void *a,const void *b)
{
    return (*(node *)a).s-(*(node *)b).s;
}

int cmpth(const void *a,const void *b)
{
    return (*(node *)a).th-(*(node *)b).th;
}

void getsa()
{
    qsort(r+1,n,sizeof(node),cmps);
    int p=0;
    for(int i=1;i<=n;++i)
        if(r[i].s!=r[i-1].s)
        {
            rank[r[i].th]=++p;
            sa[i]=r[i].th;
        }
        else
        {
            rank[r[i].th]=p;
            sa[i]=r[i].th;
        }
     qsort(r+1,n,sizeof(node),cmpth);
     for(int l=1;p<n;l<<=1)
     {
         memset(sum,0,sizeof(sum));
         memset(h,0,sizeof(h));
         for(int i=1;i<=n;++i)
             ++sum[rank[i]+1];
        for(int i=1;i<=p;++i)
            sum[i]+=sum[i-1];    
         for(int i=n-l+1;i<=n;++i)
            _sa[sum[rank[i]]+(++h[rank[i]])]=i;
        for(int i=1;i<=n;++i)
             if(sa[i]>l)
                 _sa[sum[rank[sa[i]-l]]+(++h[rank[sa[i]-l]])]=sa[i]-l;
        memcpy(sa,_sa,sizeof(_sa));
        p=0;
        for(int i=1;i<=n;++i)
            _rank[sa[i]]=((rank[sa[i]]==rank[sa[i-1]])&&(Rank(sa[i]+l)==Rank(sa[i-1]+l)))?p:++p;
        memcpy(rank,_rank,sizeof(_rank));            
     }
}

void getans()
{
    for(int i=1;i<=n;++i)
        ans[0][i]=height[i];
    for(int i=1;1<<i<=n;++i)
        for(int j=1;j+i-1<=n;++j)
            ans[i][j]=ans[i-1][j]<?ans[i-1][j+(1<<(i-1))];
}

int askRMQ(int s,int t)
{
    for(int i=0;;++i)
        if(t-s<1<<i)
            return ans[i-1][s]<?ans[i-1][t-(1<<(i-1))+1];
}

void getheight(bool flag) 
{
    for(int k=0,i=1,j;i<=n;height[rank[i++]-1]=k)
        for(k?--k:0,j=sa[rank[i]-1];r[i+k].s==r[j+k].s;++k);
     getans();
     if(flag)
     {
         tr[1]=m;
         for(int i=2;i<=m;++i)
             tr[i]=askRMQ(rank[1]<?rank[i-1],rank[1]>?rank[i-1]);
     }
     else
     {
         tl[1]=m;
         for(int i=2;i<=m;++i)
             tl[i]=askRMQ(rank[1]<?rank[i-1],rank[1]>?rank[i-1]);            
     }
}

int main()
{
    scanf("%s",a);
    scanf("%s",b);
    m=strlen(a);
    n=strlen(b);
    for(int i=0;i<n;++i)
    {
        r[i+1].s=b[i];
        r[i+1].th=i+1;
    }
    getsa();
    getheight(1);
    maxk=-1;
    for(int i=0;i<m;++i)
    {
        if(maxk<i)
        {
            for(int j=0;;++j)
                if(j==n||i+j==m||a[i+j]!=b[j])
                {
                    maxk=i+j-1;
                    maxi=i;
                    rght[i]=j;
                    break;
                }
        }
        else
        {
            if(tr[maxi-i+1]>=maxk-i+1)
                for(int j=maxk-i+1;;++j)
                    if(j==n||i+j==m||a[i+j]!=b[j])
                    {
                        maxk=i+j-1;
                        maxi=i;
                        rght[i]=j;                        
                    }
            else
                rght[i]=tr[maxi-i+1];
        }
    }
    maxk=-1;
    for(int i=0;i<n;++i)
    {
        c[i]=r[i+1].s=b[n-i-1];
        r[i+1].th=i+1;
    }
    getsa();
    getheight(0);
    for(int i=0;i<m;++i)
        d[i]=a[m-i-1];
    for(int i=0;i<m;++i)
    {
        if(maxk<i)
        {
            for(int j=0;;++j)
                if(j==n||i+j==m||d[i+j]!=c[j])
                {
                    maxk=i+j-1;
                    maxi=i;
                    lft[i]=j;
                    break;
                }
        }
        else
        {
            if(tl[maxi-i+1]>=maxk-i+1)
                for(int j=maxk-i+1;;++j)
                    if(j==n||i+j==m||d[i+j]!=c[j])
                    {
                        maxk=i+j;
                        maxi=i;
                        lft[i]=j;                        
                    }
            else
                lft[i]=tl[maxi-i+1];
        }
    }
    for(int i=m-1;i>=0;--i)
        if(lft[i]==n&&i+n!=m)
            leftt[i]=lft[i]+leftt[i+n];
        else
            leftt[i]=lft[i];
    for(int i=m-1;i>=0;--i)
        if(rght[i]==n&&i+n!=m)
            rightt[i]=rght[i]+rightt[i+n];
        else
            rightt[i]=rght[i];
    int maxl=0;
    for(int i=0;i<m;++i)
        maxl>?=rightt[i]+leftt[m-i];
    
    if(maxl<n)
        puts("0");
    else
        printf("%f\n",double(maxl)/m);
    return 0;
}

这道题是03年饶向荣论文里的一道题 有一定难度
除了T数组的求解和论文中不同外 没有什么不同
论文中说的方法没看懂 期望有牛人能讲一下
我的方法很简单就是通过后缀数组完成的 但大大增加了代码长度这好象是我写oi题目写的最长的一道了(我太弱了) 一个不错的开始期望以后每天都能AC并且要多


250 2009-04-13 23:33 发表评论
]]>100 book 完成＆近期计划http://www.cppblog.com/wwy250/archive/2009/04/05/79033.html250250Sun, 05 Apr 2009 11:09:00 GMThttp://www.cppblog.com/wwy250/archive/2009/04/05/79033.htmlhttp://www.cppblog.com/wwy250/comments/79033.htmlhttp://www.cppblog.com/wwy250/archive/2009/04/05/79033.html#Feedback4http://www.cppblog.com/wwy250/comments/commentRss/79033.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/79033.html这里的题目难度都非常大 都可以作为NOI的试题 甚至更难
其贪心、构造、调整法颇多 还有一些没听说过的 例如差分约束系统、最小限度生成树、区间图判定、最小表示法等等 还涉及到许多数论、几何知识
除了前30题 后面的题 有许多都是NP问题 或者说是无法在要去时间内出解的问题 对提交答案式问题去颇帮助
其中一些搜索题目的技巧性很高 可以达到一想不到的效果 甚至比某些多项式算法的运行时间还快
多说无意毕竟这套资料是以前集训队留下来的 大家说的一定比我强

下面3周我将作一些专项训练 大概会使用到WC论文和集训队作业 另外这次的训练要多写代码 逢题必AC(虽然不一定是自己想的)

再下面的一周也就是CTSC前的一周我将计时完成以往的CTSC、WC、APIO 从而找到手感迎接即将到来的CTSC、APIO


250 2009-04-05 19:09 发表评论
]]>100 book若干题总结http://www.cppblog.com/wwy250/archive/2009/03/30/78423.html250250Mon, 30 Mar 2009 15:12:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/30/78423.htmlhttp://www.cppblog.com/wwy250/comments/78423.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/30/78423.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/78423.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/78423.html下面对前25道总结一下 在这25道中我学会了
差分约束系统
限制最小生成树
调整法
区间图判定
最小表示法

对构造、贪心有了更进一步的认识
还学了一点数论、几何知识

另外 我发想 有些题实在不好证明就不要证明 尤其是贪心、构造、调整法 有时候即使是错的也可能又不错的收效
要尝试 一个题写多个程序 小数据搜索 大数据写一个不太完善的算法 （当然要是能写出完善的最好）
还有我的代码能力实在有待提高



250 2009-03-30 23:12 发表评论
]]>找回状态http://www.cppblog.com/wwy250/archive/2009/03/30/78357.html250250Mon, 30 Mar 2009 04:56:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/30/78357.htmlhttp://www.cppblog.com/wwy250/comments/78357.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/30/78357.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/78357.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/78357.html我们学校把科技创新等一些杂项都交给了oier  今年科技创新老师又不管 从头到尾都是自己做的 写了一个论文却写成了WC风格的
接下来的几天oi的进展都不大 没有状态100book里的每一个题基本都要看解题报告 有的看解题报告都没看懂 除非出原题 要是让我自己做一定做不上
话又说回来这两天看的题都挺偏的 坐不上也有情可原 多见一些没准就会了
老师又想让我参加吉大的ACM队 其实我不想去 吉大的ACM很水 而且又不把我分到1队 去了只是提高代码能力(虽然我的代码能力有待提高) 回头再和老师商量一下

说了这么都希望能带走不好的状态 下一阶段要全力以赴oi 离CTSC APIO就剩5周了
在这5周里 我希望能完成 100book剩下的30+ 以及以往的CTSC WC APIO 要是可能就再看一点集训队作业


250 2009-03-30 12:56 发表评论
]]>100book 0017http://www.cppblog.com/wwy250/archive/2009/03/21/77369.html250250Fri, 20 Mar 2009 20:36:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/21/77369.htmlhttp://www.cppblog.com/wwy250/comments/77369.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/21/77369.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/77369.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/77369.html一道来自 IPSC的题目
IPSC这个比赛很有意思 所有题目都是提交答案式的 可惜我英文不好 还要组队
有意与我组队请与我联系 我的邮箱：wwy250@gmail.com
还有通过这个比赛的排名还是看出中国的教育存在着巨大的问题
学生组具有垄断地位 而成人组就不行了
言归正传
很容易想到的是二分图最大匹配的问题 可以在比赛的时间内跑出来
还有一种更优的方法 对于每种S 可以在O(1)时间内接触c
对于每一种S:{a1,a2,～～,a(n+1)/2} (a1<a2<～～<a(n+1)/2) c=a((a1+a2+～～+a(n+1)/2)%((n+1)/2))
这个可以用反证法证明：

若存在两个方案删数后的方案相同，设为 A ， B （集合），

另 a1<a2< …… <am,b1<b2< …… <bm ，

若 A 中删去元素 ai ， B 中删去 bj ，不妨设 i<j ，显然有 ai<bj

根据条件， j-i=(bj-ai) mod m,

所以 bj － ai ＝ j-i 或者 j-i+m

  若 bj － ai ＝ j － i ，因为 ai+1...aj=bi...bj-1, 所以 bj-ai-1>=j-i, 所以不可能

  若 bj － ai ＝ j － i+m,

   因为若 A ， B 存在，则必满足 (m-j)+(i-1)<=(n-bj)+(ai-1)

   因为 n ＝ m ＋ m-1, 所以 bj-ai<=j-i+m-1, 该情况也不能

 

综上所述，不可能存在两个方案删数后的方案相同。




250 2009-03-21 04:36 发表评论
]]>FJOI 06http://www.cppblog.com/wwy250/archive/2009/03/17/76822.html250250Mon, 16 Mar 2009 19:41:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/17/76822.htmlhttp://www.cppblog.com/wwy250/comments/76822.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/17/76822.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/76822.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76822.html题目
第一题是一道数学问题
学过数学竞赛的人一定知道对于这个问题
分得的所有数一定不是2就是3
而且2的个数<3
证明如下
显然分成含有1的显然不是最优的
若含有p(p>3)则可将其分成由2、3的和组成的序列保证其不差于此当前方案
第二题 一道维护决策单调性的DP
定义f[i][j]在前i个点建j个站前i个点以及建站的费用最小为多上
s[i][j] 为f[i][j]的方案
显然s[i-1][j]<=s[i][j]<=s[i][j+1]
剩下的就不说了




250 2009-03-17 03:41 发表评论
]]>ZJOI 08 day1http://www.cppblog.com/wwy250/archive/2009/03/14/76608.html250250Sat, 14 Mar 2009 15:21:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/14/76608.htmlhttp://www.cppblog.com/wwy250/comments/76608.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/14/76608.html#Feedback3http://www.cppblog.com/wwy250/comments/commentRss/76608.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76608.html题目
这套题做的很囧
第一题 题目没什么说的就是模拟 不过很麻烦
但是题目里的游戏很好玩 真的很好玩
第二题 一开始看像04年 ACM上海赛区的H题 田忌赛马 记得集训队资料里的解法是O(n^2)
但是发现 n<=100000 后来发现这题与田忌赛马是不一样的
田忌赛马 是求最大|小分差 那么怎么贪心呢
以求最高分为例 提出一个策略：
设我方的选手集为A  对方为B
若Amax>Bmax
则A-=Amax,B-=Bmax  ans+=2
否则 A-=Amin,B-=Bmax 如果Amin==Bmax  ans+=1
下面是证明
若Amax>Bmax
假设有一种方案 使得Bmax不与Amax交战&得分>当前方案
设于Bmax、Amax交战的分别为a,b
则将Bmax与Amax交战 a与b交战  其余与该方案相同 易证此方案不亚于 该方案&此方案得分=原方案
若Amax<=Bmax 同上述方法可证 这里就不多说了

值得一说的是第4题 我想了2天 实在没有思路将一些想法记在下面并将它添加到未解决问题中
首先想到的是将它想LCA->RMQ一样搞出一个欧拉序列 通过维护这个序列解题
那么借助什么数据结构好呢 线段树？
这好像不可能 倒不是得到答案的问题
关键是每次更改都要不止更改一个或常数个
看来 搞成一个序列是没戏 那么仍保持树状结构呢
这回 更改时好办了 但怎么得到答案哪
我又标程 不过我这个人最不擅长就是读程序 但可以看出标程用到平衡树
实在是不会 等oibh好了到哪顶上问问应该会有结果幻灯片 20

250 2009-03-14 23:21 发表评论
]]>JSOI07 完成情况http://www.cppblog.com/wwy250/archive/2009/03/13/76458.html250250Fri, 13 Mar 2009 07:08:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/13/76458.htmlhttp://www.cppblog.com/wwy250/comments/76458.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/13/76458.html#Feedback4http://www.cppblog.com/wwy250/comments/commentRss/76458.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76458.html题目
终于知道JLOI为什么是5题4h了 因为JSOI也是 而JS给JL出题 风格当然一样
而且也会有一些比较偏的题目 还有一点是数据弱他还不告诉你 比如最后一题朴素快排就能90分 如果是NOI的会 一定会说90%的数据n<=?的而且不会是90%的 最多是40% 没办法省选又没人赞助谁给你好好出题(好像今年NOI就没有所以WC的题目所有'<='都打成了‘=’)
言归正传
这套题目好题还是有的
比如第一题 虽然我至今没搞明白 但是我知道他要求的是：从A中选取最少的点 使得B中所有点都在A中选取点的凸包内  这个变化十分巧妙
第二题 枚举和牌和对子是必然的趋势 那么剩下的判断是否为和就只能在线性时间内解决了
也就是说题目只给了我们扫一次(或常数次)的机会 而且是能按n扫
这么近的时间不得让我们想到贪心 如果对于一张牌 可以组成顺子 也可以组成刻牌 这个时候一定要有一种固定的选择
假设选择顺子 很显然若是111234 本来可以和的牌 就不胡了
那如果是刻字呢 经反复试验没有找到反例 在时间紧张的比赛中 不一定一定要证明 于是我写了一下
AC
看来我的感觉还可以 但光靠感觉是不行的 证明如下
若经过上诉贪心方法的到的答案是和牌则 这副牌一定是和牌
所以只需证明经上述算法得到的答案为非和时 这副牌一定非和 下面的证明均在经上述算法得到的答案为非和前提下进行
假设有一种方案能使这副牌和
则一定有至少一处 原方案为刻字而新方案为顺子
将每个组合按最小、较小、最大3个关键字顺次按有小到大排序
找到第一次这样的地方
显然之前的牌组合的方式两种方案是一样的
所以当将原方案中的刻字转化为顺子后 如果该方案为和牌则另两张在原方案为刻字的牌也与其后面两张组成顺子 与组成3个刻字等效 所以假设不成立(这样和在我们吉林打法还大呢)
第3题 我认为是一道比较偏的题目 看了解题报告仍觉得比较偏
一个被逼无奈的贪心 结果竟是AC 在这里我不想多说了 有兴趣的同学看这里吧
第4题 比较常见的DP 好像在URAL上做过 就是搞一个f[i][j]表示前i个字符 后缀为前缀j(这里的j只在我们预先搞好的trie里的编号)不含有可识别单词的个数 重点维护f 总之很麻烦 但好想 我就不想说了(我的表达能力容易把自己说糊涂了)
第5题 赤裸裸的后缀数组 只要将原串加倍即可 我的倍增可以AC而解题报告说这么做会超时 是不是他用string了 不超时才怪呢


250 2009-03-13 15:08 发表评论
]]>JSOI01 完成情况http://www.cppblog.com/wwy250/archive/2009/03/12/76378.html250250Thu, 12 Mar 2009 14:54:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/12/76378.htmlhttp://www.cppblog.com/wwy250/comments/76378.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/12/76378.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/76378.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76378.html题目
由于是01年的题目
难度自然比较低
前3题都是搜索/模拟题 在这里就不多累述
第4题一开始被他数据小的特点蒙骗了
搜索|状态压缩的DP 好像都不行 一时间没了头绪
后来想到了二分图 其实早应该想到二分图
以横向为例 显然对于每一条线段 如果线段上没有"墙" 则线段上对多只能有1个车
纵向同理
所以先遍历一次这个矩形 求出所有上述线段 以及所有非墙格子所在的横纵线段
将所有有相交的线段之间连一条边 求二分图最大匹配即可
对于某些所求为XX最多 每个XX影响两个元素的题目 二分图往往能够起到作用


250 2009-03-12 22:54 发表评论
]]>JLOI 告一段落http://www.cppblog.com/wwy250/archive/2009/03/12/76292.html250250Wed, 11 Mar 2009 19:26:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/12/76292.htmlhttp://www.cppblog.com/wwy250/comments/76292.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/12/76292.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/76292.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76292.html明天做江苏不过江苏我只有两届的题
但是是江苏给我们出省选试题 不做不行
各位游客如有江苏试题希望能共享资源


250 2009-03-12 03:26 发表评论
]]>JL08 棋局定式 解题报告http://www.cppblog.com/wwy250/archive/2009/03/12/76291.html250250Wed, 11 Mar 2009 19:21:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/12/76291.htmlhttp://www.cppblog.com/wwy250/comments/76291.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/12/76291.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/76291.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76291.html题目
对于这到题目很容易看出是一个多字符串匹配问题 显然是要用到trie图
不过这个题目要求将所有定式 一次trie恐怕是做不了了
难道要用其他的方法么 虽然单纯的trie无法满足题目要求 但是没有比他再像的了
继续思考trie图可以完成的事将所有匹配成功的位置以及与谁匹配返回
但是每个匹配成功的位置只能返回一个 与其成功匹配的字符串 显然这会有遗漏
比如
dd
bbdd
去匹配bbdd则最多只会有一个被记录
显然如果一个两个串需要同时被记录当且仅当一个串是另一个串的后缀并且较长串被匹配
所以每次记录长度最长的字符串即可
在这次匹配之后建一颗trie将所有字符串的逆串插入
查找所有被记录的字符串将路径上所有的字符串记录即可


250 2009-03-12 03:21 发表评论
]]>JLOI总结http://www.cppblog.com/wwy250/archive/2009/03/11/76170.html250250Tue, 10 Mar 2009 19:43:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/11/76170.htmlhttp://www.cppblog.com/wwy250/comments/76170.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/11/76170.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/76170.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76170.html题目有兴趣就看一下吧
其实JLOI的题目不怎么好 题目不是简单就是偏 好题很少 但我毕竟是JL人 JLOI还是要做的
连续做了7年的JLOI 从中看出了JLOI的成长 这一点让我很欣慰 毕竟整体水平的提高会带动个人的提高 这叫做水涨船高
但是不得不说JLOI的导向不好 水题多&时间少
强烈建议加大难度 改成two day 5h3题 这样才像NOI
以前的省选只要拿到基础分就可以进省队了 这个导向很不好 应该让真正能做上难题的人进才对
不过做在不好的题目也会有收获 下面是我的收获：
08 05 ：通过这到题从新认识了一下AC自动机 原来他不是向treap一样只要写对了就没事了
他有一些变化 所以我决定 明天看一下《多串匹配算法及其启示》《Trie图的构建、活用与改进》
05 03 ：知道了二分图最小点覆盖这个模型
另外还有了一些骗分、那部分分的心得 只可意会不可言传
另外由于时间关系08 04、06 03没有好好研究一下明天同一研究





250 2009-03-11 03:43 发表评论
]]>我未解决的问题http://www.cppblog.com/wwy250/archive/2009/03/11/76168.html250250Tue, 10 Mar 2009 18:19:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/11/76168.htmlhttp://www.cppblog.com/wwy250/comments/76168.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/11/76168.html#Feedback1http://www.cppblog.com/wwy250/comments/commentRss/76168.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76168.html这片文章被长期置顶 内容是从此帖发布以来我一直不能解决&值得研究的问题 每当我遇到不会的问题或解决了其中的某个问题都会更新它 期望大家能帮助我解 等那天得空把以前积累下来的问题也发上来

运动会

虽然拿到了满分 但是那是因为数据若 我用的是先想解决2—SAT问题贪心一个解在随机化搜索的方法

这个题有没有多项式算法？

双调路径

没看懂题题目中的"好"/"佳"没搞明白什么意思

但感觉应该是一道值得一做的题目 

块状链表

这个大概思想知道 有没有一个向讲treap一样详细的资料把他的具体实现方式以及所有的功能诠释

最小费用流 的负环情况的的处理方法

二维凸包 


合金
建筑抢修

上面两个题我都有解题报告不幸的是没看懂 大家可以去看看 看看会不会有什么启示










250 2009-03-11 02:19 发表评论
]]>oibh.org/bbs 挂了？http://www.cppblog.com/wwy250/archive/2009/03/11/76166.html250250Tue, 10 Mar 2009 17:58:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/11/76166.htmlhttp://www.cppblog.com/wwy250/comments/76166.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/11/76166.html#Feedback3http://www.cppblog.com/wwy250/comments/commentRss/76166.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76166.html经OIBH管理人员证实确实挂了


250 2009-03-11 01:58 发表评论
]]>JLOI04完成情况http://www.cppblog.com/wwy250/archive/2009/03/11/76164.html250250Tue, 10 Mar 2009 17:10:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/11/76164.htmlhttp://www.cppblog.com/wwy250/comments/76164.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/11/76164.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/76164.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76164.html题目
越往前做题目就越简单了 而且只有4道
通过省选试题难度的变化 看来JLOI这几年确实进步了不少
这次的题目全都做上了
第一题是一道数学题 在这里就不多说了
第二题他让你求的是一个点(x,y)使得sigema(wi*((x-xi)^2+(y-yi)^2))(1<=i<=n 以下同此)最小
初看这个式子里又有x 又有 y 很复杂不好权衡
可是乘法是符合结合率的
所以原式=sigema(wi*(x-xi)^2)+sigema(wi*(y-yi)^2)
所以当sigema(wi*(x-xi)^2)、sigema(wi*(y-yi)^2)都取道最小时既为所求
则就好办了就是一个带权平均数问题(这个好像有学名不过忘了)
以x为例x=sigema(wi*xi)/sigema(xi) y同理
第3题：
如果枚举任意2个为c的方块然再bfs或并查集找最长曲线显然是要超时的(O((n*m)^3)，1<=n,m<=30)
发现由于只改动两个方块而每次都求一边最长曲线会造成许多重复运算
如果现将整个矩阵先用并查集做一遍然后枚举该边哪两个为c的方块再算最长曲线就可以在常数时间内完成最长曲线的求解
如果你总是觉得用什么不对劲的地方 没关系(其实我一开始也是这样的)
我还有一种人那你放心的方法：先枚举一个c然后做一次并查集然后再枚举另一个c再用并查集求解 这回不乱了吧(能合并的最多有2条曲线) 这个方法只是将常数增大了 而让算法清晰了 不失为一种不错的选择
第4题
想了好久差点放弃了 因为只想到了O((n+m)*(n*m)^2)的算法
但后来一看1<=n,m<=10这个时间复杂度是绰绰有余的
显然任意一个被切下来的矩形都与其他矩形再无干系 这就使得这个模型没有后效性
提到后效性不由得让我们想到DP
状态为f[x1][y1][x2][y2]表示把左上角为(x1,y1) 右下角为(x2,y2)的矩形分割成题目中要求的形态最少需要的切割长度 转移只需要枚举切割线就好了


250 2009-03-11 01:10 发表评论
]]>JLOI05完成情况http://www.cppblog.com/wwy250/archive/2009/03/10/76134.html250250Tue, 10 Mar 2009 09:41:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/10/76134.htmlhttp://www.cppblog.com/wwy250/comments/76134.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/10/76134.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/76134.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76134.html题目
这套题之前没有做过 不过还是没有NOI什么也做不上的感觉
1、2、4题居然全是DP
分别是以当前位置+路径长度、以当前节点为根结点的子树+使用机器数、当前节点+当前时间为状态
其中2题的状态转移又是一次DP 方法类似背包问题
第5题 我只想到了m*n*logn^2的算法：枚举m通过二分+树状数组查找下一步的位置 虽然不能AC但打表还是可以的(最大的点用时13s)
后来在OIBH上看到了一个n^m的解法：枚举m 将问题划归成长度为n-（被删除的城市个数）北京的位置随之改变 这样既可以忽略被删除的城市
特别注意m不应定<=n
第3题 这道题总数后点收获 二分图最小点覆盖问题 可惜我从来没听说过 现在听说也不晚如果你还没听说过就去看看这个吧：http://www.matrix67.com/blog/archives/116
现在看来很简单了 搞n+m各点表示每行/列 读入01矩阵 若该点是1则将所在行与所在列连一条边
剩下的工作就是求二分图最小点覆盖了 也就是它的最大匹配


250 2009-03-10 17:41 发表评论
]]>JLOI06完成情况http://www.cppblog.com/wwy250/archive/2009/03/09/76052.html250250Mon, 09 Mar 2009 15:48:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/09/76052.htmlhttp://www.cppblog.com/wwy250/comments/76052.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/09/76052.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/76052.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76052.html题目
我发现做虽然我是JL人但是做JLOI很没有意义原因在于简单&题目偏&我很多都做过
这样好了 我就不限时做了 把认为在比赛中应该做的写出来好了
这套题好象是170+就能进省队 不过那是当年 现在300怎么也可以了吧
显然2、4是可以拿200分的
那下一个100该怎么拿呢
由于第3题叙述不清+我的文学水平较差所以我没看懂题
第5题我同样没看懂
可是第1题用2—SAT的方法初始化+随机化搜索居然拿到了满分 真是天无绝人之路



250 2009-03-09 23:48 发表评论
]]>JLOI07完成情况http://www.cppblog.com/wwy250/archive/2009/03/09/76043.html250250Mon, 09 Mar 2009 15:00:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/09/76043.htmlhttp://www.cppblog.com/wwy250/comments/76043.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/09/76043.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/76043.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76043.html题目
由于这套题很简单并且做过就没有一一去做只是把不会的做了一下
1：模拟
2：加分二叉树
3：枚举
4：数学题 维护一种平衡
当年AC的程序如下：
#include<fstream>

using namespace std;

ifstream cin("smiley.in");
ofstream cout("smiley.out");

double y64(double k)
{
    if(k<64) return k;
    else return y64(k/64);
}

int main()
{
    double i,n,m;
    int ans;
    cin>>m>>n;
    ans=0;
    i=y64(n/m);
    if((1<=i)&&(i<2))ans=0;
    if((2<=i)&&(i<4))ans=2;
    if((4<=i)&&(i<8))ans=4;
    if((8<=i)&&(i<16))ans=8;
    if((16<=i)&&(i<32))ans=16;
    if(32<=i)ans=32;
    cout<<ans;
    return 0;
}
5.唯一的有一定难度的题
其实也不难
搞一个小根堆
队中的元素为一些线段 维护向这个线段注水从开始到溢出的时间 要求每个线段的左端挡板与右端挡板均为线段中最高的
每次取根结点将它与它将溢出方向上的第一个线段合并 直至该线段为最左或最又的线段


250 2009-03-09 23:00 发表评论
]]>JLOI08完成情况http://www.cppblog.com/wwy250/archive/2009/03/09/76033.html250250Mon, 09 Mar 2009 14:09:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/09/76033.htmlhttp://www.cppblog.com/wwy250/comments/76033.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/09/76033.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/76033.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/76033.html得分 364 用时 4h
由于做过分高一点
1.模拟没什么说的 得分100
2.动态规划 得分100
令f[c][l]=为匹配部分为字符串c的前缀l使得后面的字符串有两种分解方法的最小长度
这句话说的我自己都听不懂请大家结合下图看
o_JL08 codes.bmp
所以答案即为min(f[i][0]|1<=i<=n)
转移方程为：f[c][l]=min(min(getf(i,strlen(s[c])-l)+l|后缀l+1为i的子    串),min(getf(c,l+strlen(s[i]))|i为后缀l+1的子串))
至于输出方案就不用我说了吧
3.按N皇后搜的 得分30
4.贪心+随机化+卡时 得分54
每次对于能直接拓展的每个节点有一个权值既通过该节点可直接拓展的节点数×一个随机数（0<rand()<1）
从小到大一次拓展全职由小到大的节点+上一个卡时就是这个分数
5.KMP做的 得分80分
想到 AC自动机可是我只能判断字符串集里是否有棋局的子串 而如可将所有子串求出无从下手


250 2009-03-09 22:09 发表评论
]]>线段树&DPhttp://www.cppblog.com/wwy250/archive/2009/03/09/75952.html250250Sun, 08 Mar 2009 20:09:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/09/75952.htmlhttp://www.cppblog.com/wwy250/comments/75952.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/09/75952.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/75952.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/75952.html由于这是我很长时间以前做的事了
所以只放出我搜索到的资料（大部分看了 还有一些没看懂）
这是1月前我写给自己看的 所以很乱 好多都是我硬盘里的地址
大家就将就着看吧
  阅读全文

250 2009-03-09 04:09 发表评论
]]>基础代码http://www.cppblog.com/wwy250/archive/2009/03/09/75951.html250250Sun, 08 Mar 2009 19:57:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/09/75951.htmlhttp://www.cppblog.com/wwy250/comments/75951.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/09/75951.html#Feedback0http://www.cppblog.com/wwy250/comments/commentRss/75951.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/75951.htmlhttp://www.cppblog.com/Files/wwy250/%E5%9F%BA%E7%A1%80%E4%BB%A3%E7%A0%81.rar
由于以后还有机会用到就没有一一测试
但思路一定是正确的
如果发现错误期望提出
里面 trie图错了
这样就对了

#include<iostream>

using namespace std;

struct node
{
    bool match;
    node *faild,*chaild[26];
}*trie=new node(),*super=new node(),*q[100000];
int h,l;

void insert(node *t,char *c)
{
    if(!*c)
    {
        t->match=1;
        return ;
    }
    if(!t->chaild[*c-'a'])
        t->chaild[*c-'a']=new node();
    insert(t->chaild[*c-'a'],c+1);    
}

void build()
{
    for(int i=0;i<26;++i)
        super->chaild[i]=trie;
    trie->faild=super;
    q[l=1]=trie;
    for(;h++!=l;)
    {
        node *p=q[h];
        for(int i=0;i<26;++i)
                  if(p->chaild[i])
                  {
                      p->chaild[i]->faild=p->faild->chaild[i];
                      p->chaild[i]->match|=p->chaild[i]->faild->match;
                      q[++l]=p->chaild[i];
                  }
                  else
                      p->chaild[i]=p->faild->chaild[i];
    }
}

char c[1000];

int match(node *t,int th)
{
    if(t->match)
        cout<<th<<endl;
    if(!c[th])
        return 0;
       return match(t->chaild[c[th]-'a'],th+1);
}
树状数组里ask的变量名打错了
正确的应是
#include<iostream>
#define lowbit(x) ((x)&(-(x)))

using namespace std;

int c[1000001],n;

inline void add(int p,int v)
{
    for(int i=p;i<=n;i+=lowbit(i))
        c[i]+=v;
}

inline int ask(int p)
{
    int r=0;
    for(int i=p;i>0;i-=lowbit(i))
        r+=c[i];
    return r;
}



250 2009-03-09 03:57 发表评论
]]>至NOI 09要做的事http://www.cppblog.com/wwy250/archive/2009/03/09/75950.html250250Sun, 08 Mar 2009 19:48:00 GMThttp://www.cppblog.com/wwy250/archive/2009/03/09/75950.htmlhttp://www.cppblog.com/wwy250/comments/75950.htmlhttp://www.cppblog.com/wwy250/archive/2009/03/09/75950.html#Feedback1http://www.cppblog.com/wwy250/comments/commentRss/75950.htmlhttp://www.cppblog.com/wwy250/services/trackbacks/75950.html 后缀数组
网络流&二分图
线段树
DP
基础代码:treap、树状数组（1、2维）、后缀数组、有/无（上下界的）（费用）最大流、强联通分量、kmp、AC自动机、LCA——RMQ、heep+dis、 块状链表 、高斯消元、leftheap+topsort+欧拉路 +以前写过的
3h内不看答案  最后写出来  每天3个题（可用基础代码）
100book
论文
集训队作业
sgu
ctsc
wc
预计用时：14week
熟练基础代码
预计用时：1week
第一遍 ：做题——每题记时做 4h内 记录通过调试过样例时间、通过初步测试时间、通过对拍时间 并分别保存测试记录
第二遍 ：非AC 每题记录 通过调试过样例时间、通过初步测试时间 并分别保存测试记录（时间×1.5，得分×70%）
第三遍 ：看解题报告
每天2个题
制定NOI2009方案
noi（不可用基础代码）
预计用时：4week
强化线段树、dp&调整作息时间
预计用时：1week

我会在完成每项后放出所用到的资料名称以及我的原创资料 最前面已经是一个月之前的事了我放不了那么全了
注释：绿色为已完成


250 2009-03-09 03:48 发表评论
]]>