Thu, 22 Mar 2007 15:14:00 GMT

题目�Q?br />
�?到N(100000)中�Q意拿掉两个数�Q�把剩下�?9998个数��序打�ؕ�Q��ƈ且放入数�l�A中。要求只扫描一遍数�l�，把这两个数找出来。可以��用最��C��过5个局部变量，不能用数�l�变量，�q�且不能改变原数�l�的倹{�?br />
思�\�Q?br />
遍历一�ơ数�l�，求出�q�两个数的和a+b 与积a*b
a+b = 1+2+3+4+...+N- sum(A[]);    (1)
a*b = 1*2*3*4*...*N / multi(A[]);   (2)

主要解决sum与multi的溢出问�?br />
(1) 可化�?(N-A[0]) + (N-1-A[1]) + ...+ (3-A[N-3]) + 2 + 1

(2) 可以用对数来代替原数�q�行求积的等仯��?避免溢出的问题，但是�q�种�Ҏ(gu��)��?x��)��生一些精度上的问题，不知道大家有什么更好的�Ҏ(gu��)��Q?br />先求出log(a*b) :
                         = log(1*2*3*4*....*N)/log(A[0]*A[1]*A[2]*...*A[N-3])
                         = log(N)-log(A[0]) + log(N-1)-log(A[1]) + ... +log(3)-log(A[N-3]) + log(2) + log(1)

知道了两数的和与�U?由此��可以计��出a跟b的值来.

代码如下�Q?br />

#include <iostream>
#include <Ctime>
#include <Cmath>
using namespace std;

#define N 100000

//生成不同的随机数的数�l?/span>
void GetDiffRandomNum(int A[], int n)
{
    srand(unsigned(time(NULL)));
    int i=0;

    for(int index = n-1; index > 0; index--)
    {
        i = rand() % index;
        swap(A[i], A[index]);
    }

}

int main()
{

    int A[N]={0};
    for(int i=0; i<N; i++)
    {
        A[i] = i+1;
    }
    GetDiffRandomNum(A, N);
    //DISPLAY(A, N);

    unsigned int sum = 0;
    double logSum = 0;

    for(i=0; i<N-2; i++)
    {
        sum += N-i-A[i];
        logSum += log(N-i)-log(A[i]);
    }
    sum += 2 + 1;
    logSum += log(2)+log(1);

    double multi = exp(logSum);

    //两数的和与积
    cout<<int(sum)<<'\t'<<int(multi)<<endl;

    //求出两数
    for(i=1; i<=N; i++)
    {
        double temp = i*(sum-i);
        if(multi-0.5<=temp && temp <= multi+0.5)
            cout<<i<<'\t'<<int(sum-i)<<endl;
    }

    return 0;
}

PS�Q�谢谢枝~的帮助）(j��)请大家指�?br />
//................................
通过大家的帮助：(x��)
得到另一个写法，不会(x��)产生�_�ֺ�问题

(1+N)*N /2 - S = a + b
1/6 * n*(n + 1)*(2n + 1) - X = a*a + b*b

注：(x��)
1/6 * n*(n + 1)*(2n + 1)=1*1 + 2*2 + 3*3 +...+N*N
X = A[0]*A[0] + A[1]*A[1] +...A[N-3]*A[N-3]

==>
a + b = m
a*a + b*b = n

�׃��可解出a,b

    unsigned int sum = 0;
    unsigned int sqrSum = 0;

    for(i=0; i<N-2; i++)
    {
        sum += N-i-A[i];
        sqrSum += ((N-i)*(N-i)) - ((A[i])*A[i]);

    }
    sum += 2 + 1;
    sqrSum += 2*2 + 1*1;

猪头�?/a> 2007-03-22 23:14 发表评论

KMP��法资料

Thu, 09 Nov 2006 17:51:00 GMT

看了两三天的KMP��法,一直看的迷�q�L(f��ng)��p�的.现在把这些资料脓(chu��ng)在这�?..以备日后之需

1.串的模式匚w��的改�q�算�?�q�个�|�站�Ҏ(gu��)��的理解帮助很�?特别是右边的那块说明部分,以前自己脑筋老是转不�q�来) http://cist.dhu.edu.cn/kejian/%CA%FD%BE%DD%BD%E1%B9%B9%BE%AB%C6%B7%BF%CE%B3%CC/%D4%DA%CF%DF%D1%A7%CF%B0/text/chapter04/section3/c5.htm

2.KMP ��法的注�?http://www.cublog.cn/u/20/showart_136705.html

3.KMP��法中推导next[],nextval[]--手记 http://jiasimon040510.t8log.ccut.cn/blog-htm-do-showone-tid-6983.html

4.��法原理�Q?

在匹配过和中�Q�当��M��中第i个字�W�与模式串中�W�j个字�W�“失配”时�Q�s[i]!=t[j])�Q�将模式串尽量向右移动，让模式串中第k�Q�k

要让�q�个条�g成立�Q�那么在k之前的k个t字符[0 �?k-1]必须在i之前的k个s字符[i-k �?i-1]相匹配即�Q?/p>

t[0, 1, 2...k-1] == s[i-k, i-k+1, i-k+2...i-1] ---(1)

而由之前的部分匹配成功的�l�果可知�Q?br />
   t[0, 1, 2...j-1] == s[i-j, i-j+1, i-j+2...i-1]     ---(2)
==>
   t[j-k, j-k+1, j-k+2...j-1] == s[i-k, i-k+1, i-k+2...i-1]   --(3)

�?1)�?3)可得�Q?/p>

t[0, 1, 2...k-1] == t[j-k, j-k+1, j-k+2...j-1] ---(4)

求出k�?��是next[j]的��g��

��M��,相对我来�?��法不是很好�?但是大家看到我这么笨的�h到最后都能明白一�?大家��更没有理由看不懂了,��大家成功附上我的测试源�?

#include < iostream >

using namespace std;

void GetNext( char t[], int next[])
{
     int j = 0 ;
     int k = - 1 ;
    next[j] = k;
     int tlen = strlen(t);

     while (j < tlen)
     {
         if (k == - 1 || t[j] == t[k])
         {
            j ++ ;
            k ++ ;
             if (t[j] == t[k])
             {
                next[j] = next[k];
            }
             else
                next[j] = k;
        }
         else
         {
            k = next[k];
        }
    }
}

int KMP( char s[], char t[], int pos, int next[])
{
     int slen = strlen(s);
     int tlen = strlen(t);
     int i = 0 ;
     int j = 0 ;

     while (i < slen && j < tlen)
     {
         if (j == - 1 || s[i] == t[j])
         {
            i ++ ;
            j ++ ;
        }
         else
         {
            j = next[j];
        }
    }

     if (j == tlen)
     {
         return i - tlen;
    }
     else
         return - 1 ;
}

int main ( int argc, char ** argv)
{

     char s[] = " aaaabaabaaabaaabaaaaabaaabaaabaaabaaabaaabaaabaaabaaabaaabacb " ;

     char t[] = " aabaaa " ;

     int next[ 20 ] = { 0 } ;
    GetNext(t, next);
     for ( int i = 0 ; i < 20 ; i ++ )
        cout << " next[ " << i << " ]:   " << next[i] << endl;

    cout << KMP(s, t, 0 , next) << endl;
}

猪头�?/a> 2006-11-10 01:51 发表评论

国产综合av,亚洲国产成人av好男人在线观看,美女网站在线免费欧美精品

KMP���法资料

KMP��法资料