编程珠玑及续笔记

编程珠玑
15.2 最长重复子串

#include<iostream>
#include<algorithm>
using namespace std;    
const int MAXN=5000000;
char c[MAXN],*a[MAXN];

int comlen(char *p,char *q)
{
    int i=0;
    while(*p&&(*p++==*q++))++i;
    return i;
}

int pstrcmp(const void *a,const void *b)
{
    return strcmp(*(char **)a,*(char**)b);
}                                                    

int main()
{
    char ch;
    int n=0,maxlen=0,maxi=0;
    while(cin>>ch)
    {
        a[n]=&c[n];
        c[n++]=ch;
    }
    c[n]='\0';
    qsort(a,n,sizeof(char *),pstrcmp);
    for(int i=0;i<n-1;++i)
    {
        if(comlen(a[i],a[i+1])>maxlen)
        {
            maxlen=comlen(a[i],a[i+1]);
            maxi=i;
        }
    }
    cout<<maxlen<<' '<<a[maxi]<<endl;
    return 0;
}

qsort函数的用法：http://www.cnblogs.com/syxchina/archive/2010/07/29/2197382.html

12 取样问题
蓄水池抽样
先把1到k放入数组，从i=k+1开始，以k/i的概率选择i，若i被选择，则随机替换数组中的一个数。最终每个数被选中的概率均为k/n。
证：
i被选中的概率 = 选择i的概率 * (i后面元素不被选择的概率 + i后面元素被选择的概率 * 不替换i的概率)

$$\begin{array}{l} p = \frac{k}{i} \times \left[ {\left( {\frac{{i + 1 - k}}{{i + 1}} + \frac{k}{{i + 1}} \times \frac{{k - 1}}{k}} \right) \times \cdots \times \left( {\frac{{n - k}}{n} + \frac{k}{n} \times \frac{{k - 1}}{k}} \right)} \right]\\ = \frac{k}{i} \times \left( {\frac{i}{{i + 1}} \times \cdots \times \frac{{n - 1}}{n}} \right) = \frac{k}{n} \end{array}$$

vector<int> f(vector<int> vec,int k)
{
    srand((unsigned)time(0));
    vector<int> result;
    for(int i=0;i<k;++i)result.push_back(vec[i]);
    for(int i=k;i<vec.size();++i)
    {
        int j=rand()%(i+1);
        if(j<k)result[j]=vec[i];
    }
    return result;
}

概率随机问题【1】【2】【3】：http://blog.csdn.net/hxz_qlh/article/details/12843131

编程珠玑（续）
第1章（性能监视工具）

//程序P5
bool prime(int n)
{
    if(n%2==0)return n==2;
    else if(n%3==0)return n==3;
    else if(n%5==0)return n==5;
    for(int i=7;i*i<=n;i+=2)
    {
        if(n%i==0)return false;
    }
    return true;
}

//程序P6
#include<iostream>
#include<vector>
using namespace std;

bool prime(int n,vector<int> &x)
{
    for(int i=0;x[i]*x[i]<=n;++i)
    {
        if(n%x[i]==0)return false;
    }
    x.push_back(n);
    return true;
}

int main()
{
    vector<int> x(1,2);
    cout<<2<<endl;
    for(int i=3;i<=1000;++i)
    {
        if(prime(i,x))cout<<i<<endl;
    }
    return 0;
}

第13章（绝妙的取样）

//算法S
unordered_set<int> fun1(int M,int N)
{
    unordered_set<int> S;
    for(int size=0;size<M;)
    {
        int T=rand()%N+1;
        if(S.find(T)==S.end())
        {
            S.insert(T);
            ++size;
        }
    }
    return S;
}

//算法F1
unordered_set<int> Sample(int M,int N)
{
    if(M==0)
    {
        unordered_set<int> S;
        return S;
    }
    unordered_set<int> S=Sample(M-1,N-1);
    int T=rand()%N+1;
    if(S.find(T)==S.end())S.insert(T);
    else S.insert(N);
    return S;
}

//算法F2
unordered_set<int> fun2(int M,int N)
{
    unordered_set<int> S;
    for(int J=N-M+1;J<=N;++J)
    {
        int T=rand()%J+1;
        if(S.find(T)==S.end())S.insert(T);
        else S.insert(J);
    }
    return S;
}

//算法P
vector<int> f(int N,int M)
{
    vector<int> S;
    for(int J=N-M+1;J<=N;++J)
    {
        int T=rand()%J+1;
        vector<int>::iterator iter=find(S.begin(),S.end(),T);
        if(iter==S.end())S.insert(S.begin(),T);
        else S.insert(iter+1,J);
    }
    return S;
}