8月10号小练

网站：csust训练3

A 三分的题目，给出灯高，人高，以及和墙的距离，求最长的影子长度。   Light Bulb     ZOJ 3203    

本来我是三分角度，角度[0,atan(H/D)]，测试数据都是对的，但是WA了，我到现在还是不知道为什么........

后来是三分投在墙上的影子长度，[0,h]，投在墙上的影子最多是h。

代码：     0MS

#include <stdio.h>
#include <iostream>
#include <math.h>
using namespace std;
int main()
{
    double H,h,D,l,r,mid,midd,l1,l2;
    int T;
    scanf("%d",&T);
    while (T--)
    {
        scanf("%lf%lf%lf",&H,&h,&D);
        r=h;
        l=0;
        while(r-l>1e-9)
        {
            mid=(l+r)/2;
            midd=(mid+r)/2;
            l1=mid+D*(h-mid)/(H-mid);
            l2=midd+D*(h-midd)/(H-midd);
            if(l1<l2)
                l=mid;
            else
                r=midd;
        }
        printf("%.3lf
",l2);
    }
    return 0;
}

这道题还有个数学方法，直接求峰值，详情见：http://blog.csdn.net/freezhanacmore/article/details/9887115 （表示膜拜，用数学方法的都膜拜）

B   求最长回文子串的长度，如果暴力0(n^2）果断超时，要用Manacher算法。

优化原因见：  http://www.2cto.com/kf/201210/164253.html

代码： 360 ms/15000ms     10936kb/65536kb     这道题的时间限制真长......从没做过15s的题Orz

#include<iostream>
#include<cstring>
#include<cstdio> 
using namespace std;
#define N 1000008
char str[N], str1[N<<1];
int p[N<<1];
int maxx;
void Manacher()
{
    int i, j, mx, id;
    memset(str1, '#', sizeof(str1));
    memset(p, 0, sizeof(p));    
    for(i=0; str[i]!=''; i++)
        str1[(i+1)<<1]=str[i];    //把字符插入到#####中
    str1[(i+1)<<1]='';   
    mx=0;
    maxx=0;
    for(i=1; str1[i]!=''; i++)
    {
       if( mx>i ) 
           p[i]=min(mx-i, p[id*2-i]);    //见：http://www.cnblogs.com/lv-2012/archive/2012/11/15/2772268.html    PS:没有这一步会超时
        else 
            p[i]=1; 
        while( str1[i-p[i]]==str1[i+p[i]] )    //拓展半径
            p[i]++;        
        if( i+p[i]>mx )//mx代表当前回文串扩展的最大距离 
        {
            mx=i+p[i];  
            id=i;
        }         
        if( p[i]-1>maxx )    
            maxx=p[i]-1;    //更新最长长度
    }    
}
int main()
{
    int Case=0;
    while( scanf("%s", str) && strcmp(str, "END")!=0 )   //str！=END
    {
        maxx=0;
        Manacher();
        printf("Case %d: %d
",++Case,maxx); 
    }
} 

C    Minimum Inversion Number   HDU 1394   大意是给n个数，分别是0~n-1，前面的数比后面的数大的数对有多少个，算好之后把第一个数移到最后，又算.......最后求最小的数对是多少个。此题用线段树做，但是我还没看线段树，过一阵再补上，先讲下不用线段树的做法。

代码  ：   265ms

#include <stdio.h>
#include <iostream>
using namespcae std;
int main()
{
    int n,i,j,sum,minn,sum1
    while(~scanf("%d",&n))
    {
        sum=0;
        for(i=0;i<n;i++)
            scanf("%d",&a[i]);
        for(i=0;i<n-1;i++)
            for(j=i+1;j<n;j++)
                if(a[i]>a[j])
                    sum++;   //求出初数列的个数
        minn=sum;   //初始化最小值
        for(i=1;i<n;i++)
        {
            sum1=sum+n-1-2*a[i];     //sum1=sum(上一个数列的值）+（n-1-a[i])（移到后面的那个数和最大数的差值）-a[i](比a[i]小的数的个数）
            if(minn>sum1)
                minn=sum1;    //更新最小值
            sum=sum1;   //更新上一个数列的值
        }
        printf("%d
",minn);
    }
    return 0;
}

D     Line belt   HDU 3400   三分叠三分

代码：     0MS ，注意不能用while，用do....while，因为至少要执行一次，算AB,CD 间的距离

#include<stdio.h>
#include<iostream>
#include<math.h>
#include<iostream>
#include<algorithm>
#include<string.h>
using namespace std;
struct point
{
    double x,y;
};
double dis(point a,point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));    //求两点间距离
}
double p,q,r;
double find2(point a,point c,point d)//三分CD
{
    point left,right;
    point mid,midmid;
    double t1,t2;
    left=c;
    right=d;
    do                                 
    {
        mid.x=(left.x+right.x)/2;
        mid.y=(left.y+right.y)/2;
        midmid.x=(mid.x+right.x)/2;
        midmid.y=(mid.y+right.y)/2;
        t1=dis(a,mid)/r+dis(mid,d)/q;
        t2=dis(a,midmid)/r+dis(midmid,d)/q;    //ids(a,mid)是求AB,CD间的距离
        if(t1>t2)
           left=mid;
        else
            right=midmid;
    }while(dis(left,right)>1e-5);   
    return t2;
}
double find(point a,point b,point c,point d)   //三分AB
{
    point left,right;
    point mid,midmid;
    double t1,t2;
    left=a;
    right=b;
    do
    {
        mid.x=(left.x+right.x)/2;
        mid.y=(left.y+right.y)/2;
        midmid.x=(mid.x+right.x)/2;
        midmid.y=(mid.y+right.y)/2;
        t1=dis(a,mid)/p+find2(mid,c,d);
        t2=dis(a,midmid)/p+find2(midmid,c,d);
        if(t1>t2)
        left=mid;
        else
            right=midmid;
    }while(dis(left,right)>=1e-5);
    return t2;
}
int main()
{
    int T;
    point a,b,c,d;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%lf%lf%lf%lf",&a.x,&a.y,&b.x,&b.y);
        scanf("%lf%lf%lf%lf",&c.x,&c.y,&d.x,&d.y);
        scanf("%lf%lf%lf",&p,&q,&r);
        printf("%.2lf
",find(a,b,c,d));
    }
    return 0;
}

E   超级简单，母函数解决 Fruit   HDU 2152

代码：    0ms

#include <iostream>
#include <stdio.h>
#include <string.h>
using namespace std;
int a[105],b[105];
int c1[10005],c2[10005];
int main()
{
    int n,m,i,j,k;
    while(~scanf("%d%d",&n,&m))
    {
        for(i=1;i<=n;i++)
            scanf("%d%d",&a[i],&b[i]);
        memset(c1,0,sizeof(c1));
        memset(c2,0,sizeof(c2));
        c1[0]=1;
        for(i=1;i<=n;i++)
        {
            for(j=0;j<=m;j++)
                for(k=a[i];j+k<=m&&k<=b[i];k++)
                c2[j+k]+=c1[j];
            for(j=0;j<=m;j++)
                c1[j]=c2[j];
            memset(c2,0,sizeof(c2));
        }
        printf("%d
",c1[m]);
    }
    return 0;
}