ZOJ Monthly, June 2012

ZOJ Monthly, June 2012

http://acm.zju.edu.cn/onlinejudge/showContestProblems.do?contestId=338

B题   我用的是二分 + 树状数组 ，x&(-x) 写成了 x&（-x）弄了半天

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
const int maxn=100010;
#define lowbit(x) x&(-x)

ll tab[maxn],num[maxn],op[maxn];
int Max,tot;
char ord[maxn][10];
int cnt[maxn];

int bin1(ll key)
{
    int l=1,r=Max,mid;
    while(l<=r)
    {
        mid=(l+r)>>1;
        if(tab[mid]==key) return mid;
        else if(tab[mid] > key) r=mid-1;
        else l=mid+1;
    }
}
void insert(int pos,int d)
{
    for(int i=pos;i<=Max;i+=lowbit(i))
      num[i]+=d;
}
int query(int pos)
{
    int sum=0;
    for(int i=pos;i>0;i-=lowbit(i))
      sum+=num[i];
    return sum;
}
long long bin2(int s)
{
    int l=1,r=Max,mid,ans=Max;
    while(l<=r)
    {
        mid=(l+r)>>1;
        int tmp=query(mid);
        if(tmp>=s)
        {
            ans=mid,r=mid-1;
        }
        else  l=mid+1;
    }
    return tab[ans];
}
int main()
{
    int ca,n,k;
    scanf("%d",&ca);

    while(ca--)
    {
        scanf("%d",&n);
        for(int i=0;i<n;i++)
        {
            scanf("%s %lld",ord[i],&op[i]);
            tab[i+1]=op[i];
        }
        sort(tab+1,tab+n+1);
        Max=1;
        for(int i=2;i<=n;i++)
         if(tab[i]!=tab[i-1]) tab[++Max]=tab[i];
        memset(num,0,sizeof(num));
        memset(cnt,0,sizeof(cnt));
        tot=0;
        //cout<< Max <<endl;
        for(int i=0;i<n;i++)
        {
            int pos=bin1(op[i]);
            //cout<<"bug" << pos <<" "<< query(pos) <<endl;
            if(ord[i][0]=='a')
            {
                tot++;
                insert(pos,1);
                cnt[pos]++;
                if(tot%2)
                {
                    printf("%lld\n",bin2(tot/2+1));
                }
                else
                {
                    long long a=bin2(tot/2),b=bin2(tot/2+1);
                    //cout<< a <<" "<< b<< " " << query(3)<<endl;
                    if(((a+b)%2+2)%2==0) printf("%lld\n",(a+b)/2);
                    else printf("%.1lf\n",(a+b)*1.0/2);
                }
            }
            else
            {
                int now=query(pos),pre=query(pos-1);
                //cout<< pre <<" "<<now<<endl;
                if(cnt[pos]==0) puts("Wrong!");
                else
                {
                    tot--;
                    insert(pos,-1);
                    cnt[pos]--;
                    if(tot==0) puts("Empty!");
                    else
                    {
                        if(tot%2)
                        {
                            printf("%lld\n",bin2(tot/2+1));
                        }
                        else
                        {
                                long long a=bin2(tot/2),b=bin2(tot/2+1);

                                if(((a+b)%2+2)%2==0) printf("%lld\n",(a+b)/2);
                                else printf("%.1lf\n",(a+b)*1.0/2);
                        }
                    }
                }
            }
        }
    }
    return 0;
}
/*
8
add 1
add 1
add 1
add 3
add 3
add 4
add 2
remove 2
*/

J题  这个题很简单 ，很常见的状态压缩dp，d[state][node] 表示遍历了state 状态 到达node 节点所花费的最短时间

#include <iostream>
#include <cstring>
#include <cstdio>

using namespace std;
const int INF=100000000;
int n,m,T;
int value[10],start,end,dist[10][10];
int ans;
int dp[1024][11];

void floyed()
{
    /* for(int i=0;i<n;i++)
     for(int j=0;j<n;j++)
        dist[i][j]=map[i][j];*/
     for(int i=0;i<n;i++)
     for(int j=0;j<n;j++)
     for(int k=0;k<n;k++)
     {
         dist[i][j]=min(dist[i][j],dist[i][k]+dist[k][j]);
     }
}

int main()
{
    int a,b,c;
    while(scanf("%d%d%d",&n,&m,&T)==3)
    {
        scanf("%d%d",&start,&end);
        for(int i=0;i<n;i++) scanf("%d",&value[i]);
        for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
          dist[i][j]=INF;
        for(int i=0;i<m;i++)
        {
            scanf("%d %d %d",&a,&b,&c);
            if(dist[a][b]>c) dist[a][b]=dist[b][a]=c;
        }
        for(int i=0;i<n;i++) dist[i][i]=0;
        floyed();
        ans=0;
        for(int j=0;j<(1<<n);j++)
        for(int i=0;i<n;i++)
          dp[j][i]=INF;
        dp[1<<start][start]=0;
        int limit=1<<n;
        for(int i=0;i<limit;i++)
        for(int j=0;j<n;j++)
          if(dp[i][j]<=T)
          {
              for(int k=0;k<n;k++)
                dp[i|(1<<k)][k]=min(dp[i|(1<<k)][k],dp[i][j]+dist[j][k]);
          }
         //int ans
         for(int i=0;i<limit;i++)
          if(dp[i][end]<=T)
          {
              int sum=0;
              for(int j=0;j<n;j++)
                if(i&(1<<j)) sum+=value[j];
              if(ans<sum)  ans=sum;
          }
          printf("%d\n",ans);
    }
    return 0;
}

最蛋疼的就是K了，浪费了我这么多时间，还没有做出来，思路我想出来了，就是把题意理解的错误，0 ≤ s < 2^63, 2 ≤ k ≤ 62 ，这个s如果用十进制的组多2^63,我还以为long long 也会超，都没敢再做；

#include <iostream>
#include <vector>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
typedef unsigned long long ll;
const ll inf=1LL<<63;

int a[110];
vector<pair<int,int> > e[110];
ll sum[110];

int main()
{
    ll s;
    char str[110];
    int k;
   for(int i=2;i<=62;i++){
        if(a[i]==0){
            e[i].push_back(make_pair(i,1));
            for(int j=i*2;j<=62;j+=i){
                a[j]=1;
                int t=j,cnt=0;
                while(t%i==0){
                    t/=i;
                    cnt++;
                }
                e[j].push_back(make_pair(i,cnt));
            }
        }
    }
    //for(int i=0;i<e[62].size();i++)
    //      cout<< e[62][i].first<<" "<<e[62][i].second<<endl;
    while(cin>>str>>k)
    {
        int len=strlen(str);
        s=0;int p;
        for(int i=0;i<len;i++){
            if('0'<=str[i]&&str[i]<='9')
                p=str[i]-'0';
            else if('A'<=str[i]&&str[i]<='Z')
                p=str[i]-'A'+10;
            else if('a'<=str[i]&&str[i]<='z')
                p=str[i]-'a'+36;
            s=s*k+p;
        }
        //cout<< s <<endl;
        ll ans=inf;
        for(int i=0;i<e[k].size();i++)
        {
            int p=e[k][i].first;
            ll N=s,cnt=0;
            while(N)
            {
                cnt+=N/p;
                N/=p;
            }
            ans=min(ans,cnt/e[k][i].second);
        }
        cout<< ans <<endl;
    }
    return 0;
}

D题 还有没有ac，不知道怎么错了，先贴出了以后再看；

#include <iostream>
#include <cstdio>

using namespace std;
const int maxn=510;

int n,m;
double matrix[maxn][maxn],sum1[maxn][maxn],sum2[maxn][maxn];
double heig[maxn][maxn],heigh[maxn][maxn];
int que[maxn];

int main()
{
    int Q,a,b,ca=1;
    while(scanf("%d %d",&n,&m)==2)
    {
        for(int i=0;i<n;i++)
        for(int j=0;j<m;j++)
        {
           scanf("%lf",&matrix[i][j]);
           sum1[i][j]=sum2[i][j]=0;
        }
        for(int i=0;i<n;i++)
        {
            double tot=0;
            for(int j=0;j<m;j++)
            {
               tot+=matrix[i][j];
               sum1[i][j]= (i==0?0:sum1[i-1][j])+tot;
            }
        }
        for(int i=0;i<n;i++)
        {
            double tot=0;
            for(int j=0;j<m;j++)
            {
               tot+=matrix[i][j]*matrix[i][j];
               sum2[i][j]=(i==0?0:sum2[i-1][j])+tot;
            }
        }
       // cout<<sum1[1][1]<<" "<<sum2[2][0]<<endl;
        printf("Case %d:\n",ca++);
        scanf("%d",&Q);
        double ans=99999999;
        int sx,sy;
        while(Q--)
        {
            scanf("%d %d",&a,&b);
            ans=99999999;
            for(int i=0;i<n;i++)
            {
                int l=0,r=-1;
                for(int j=0;j<m;j++)
                {
                    while(l<=r&&matrix[i][que[r]] < matrix[i][j]) r--;
                    que[++r]=j;
                    if(j>=b-1) heig[i][j-b+1]=matrix[i][que[l]];
                    if(r-l+1>=b) l++;
                }
            }
            //cout<< heig[1][1] <<endl;
            for(int i=0;i<m;i++)
            {
                int l=0,r=-1;
                for(int j=0;j<n;j++)
                {
                    while(l<=r&&heig[que[r]][i] < heig[j][i]) r--;
                    que[++r]=j;
                    if(j>=a-1) heigh[j-a+1][i]=heig[que[l]][i];
                    if(r-l+1>=a) l++;
                }
            }
            //cout<< heigh[2][3] <<endl;
            for(int i=0;i<=n-a;i++)
            for(int j=0;j<=m-b;j++)
             {
                 double tmp1=(sum2[i+a-1][j+b-1]- (j==0?0:sum2[i+a-1][j-1]) - (i==0?0:sum2[i-1][j+b-1]) + ((i==0||j==0)?0:sum2[i-1][j-1]) -heigh[i][j]*heigh[i][j])/(a*b-1);
                 double tmp2=(sum1[i+a-1][j+b-1]- (j==0?0:sum1[i+a-1][j-1]) - (i==0?0:sum1[i-1][j+b-1]) + ((i==0||j==0)?0:sum1[i-1][j-1])-heigh[i][j])/(a*b-1);
                 tmp2*=tmp2;
                 if(tmp1-tmp2 < ans)
                 {
                     ans=tmp1-tmp2;
                     sx=i,sy=j;
                 }
                 else if(tmp1-tmp2 == ans&&sx>i)
                 {
                    // ans=tmp1-tmp2;
                     sx=i,sy=j;
                 }
                 else if(tmp1-tmp2 == ans&&sx==i&&sy>j)
                 {
                     //ans=tmp1-tmp2;
                     sx=i,sy=j;
                 }
             }
             printf("(%d, %d), %.2lf\n",sx+1,sy+1,ans);
        }
    }
    return 0;
}

D 题终于知道错在哪里了，一个傻逼错误wa的我半死不活了，bug半天也没看出来，跟ac程序对拍才发现

  if(r-l+1>=b) l++;

写错了，看来模版题还是要多敲呀多记呀