2012 ACM/ICPC Asia Regional Online Warmup 部分代码和数据

总结见:http://blog.csdn.net/woshi250hua/article/details/7908090


1002代码:

#include <stdio.h>
#include <string.h>
#define MAX 1110000
#define int64 __int64//long long


struct point{

    int64 y,x;
}pot[MAX];
int n,m,arr[MAX];
int qu[MAX],head,tail;
int64 sum[MAX],dp[MAX];


int CheckIt(int x,int y,int z) {

    point p0 = pot[x],p1 = pot[y],p2 = pot[z];
    return (p0.x - p1.x) * (p0.y - p2.y) - (p0.y - p1.y) * (p0.x - p2.x) < 0;
}
int NotBest(int x,int y,int k) {

  
     point p0 = pot[x],p1 = pot[y]; 
     //printf("%lld,%lld   %lld,%lld\n",p0.x,p0.y,pot[1].x,p1.y);
     return p0.y - k * p0.x > p1.y - k * p1.x;
}
int GetInt() {

    char ch = ' ';
    while (ch < '0' || ch > '9')
        ch = getchar();
    int x = 0;
    while (ch >= '0' && ch <= '9')
        x = x * 10 + ch - '0',ch = getchar();
    return x;
}
int64 Solve_DP() {

    int i,j,k;


    head = 0,tail = 0;
    //pot[0].x = sum[1],pot[0].y = sum[1] * sum[1];
    for (i = 1; i <= n; ++i) {

        pot[i].x = sum[i];
        pot[i].y = dp[i-1] + sum[i] * sum[i];
        while (head <= tail - 1
                && CheckIt(qu[tail-1],qu[tail],i)) tail--;
        qu[++tail] = i;


        while (head + 1 <= tail
                && NotBest(qu[head],qu[head+1],2 * sum[i])) head++;
        k = qu[head];
        dp[i] = pot[k].y - 2 * sum[i] * pot[k].x + m + sum[i] * sum[i];
      
    }
    return dp[n];
}


int main()
{
    int i,j,k;


    while (scanf("%d%d",&n,&m),n+m) {

        for (i = 1; i <= n; ++i)
            sum[i] = GetInt();//scanf("%d",&sum[i]);


        int64 ans = Solve_DP();
        printf("%I64d\n",ans);
    }
}

1003代码:

#include <stdio.h>
#include <string.h>
#define MAX 10000
#define int64 __int64


int used[MAX],n,m;
int sum[MAX],col[MAX];


int64 Gcd(int64 x,int64 y) {

    int r = x % y;
    while (r) {

        x = y,y = r;
        r =  x % y;
    }
    return y;
}
int64 Lcm(int64 x,int64 y) {

    return x / Gcd(x,y) * y;
}
int Cal(int x) {

    int i,j,k,cnt;
    cnt = used[x] = 1;


    i = (x + m - 1) / m;
    j = x % m ? x % m : m;
    k = sum[j-1] + (col[j] - i + 1);
    //printf("%d\n",k);
    while (!used[k]) {

        used[k] = 1;
        cnt++;
        i = (k + m - 1) / m;
        j = k % m ? k % m : m;
        k = sum[j-1] + (col[j] - i + 1);
    }
    return cnt;
}

int main()
{
    int i,j,k;


    while (scanf("%d%d",&n,&m),n + m) {

        memset(used,0,sizeof(used));
        int64 ans = 1;
        for (i = 1; i <= m; ++i) {

            col[i] = n / m;
            if (i <= n % m) col[i]++;
            sum[i] = sum[i-1] + col[i];
        }


        for (i = 1; i <= n; ++i)
            if (!used[i])  {

                int64 tp = Cal(i);
                ans = Lcm(ans,tp);
            }


        printf("%I64d\n",ans);
    }
}


1005代码:

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <queue>
using namespace std;
#define MAX 2100
#define INF (1<<30)
#define min(a,b) ((a)<(b)?(a):(b))


int n,m,arr[MAX];
int dp[MAX][MAX];
int cost[MAX][MAX];
priority_queue<int> high;
priority_queue<int,vector<int>,greater<int> > low;


void Calculate(int beg) {

    int i,j,k,tp;
    int sumlow = 0,sumhigh = 0;


    cost[beg][beg] = 0;
    if (arr[beg] > arr[beg + 1]) {

        low.push(arr[beg]);
        high.push(arr[beg+1]);
    }
    else {

        low.push(arr[beg+1]);
        high.push(arr[beg]);
    }
    cost[beg][beg+1] = low.top() - high.top();
    sumlow += low.top(),sumhigh += high.top();


    for (i = beg + 2; i <= n; ++i) {

        if (arr[i] >=  high.top()) {

            low.push(arr[i]);
            sumlow += arr[i];
        }
        else {

            high.push(arr[i]);
            sumhigh += arr[i];
        }


        int hsize = high.size(),lsize = low.size();
        int htop = high.top(),ltop = low.top();
        if (hsize - lsize == 2) {

                sumlow += htop;
                sumhigh -= htop;
                lsize++,hsize--;
                low.push(htop),high.pop();
         }
        else if (low.size() - high.size() == 2) {

                sumlow -= ltop;
                sumhigh += ltop;
                lsize--,hsize++;
                high.push(ltop),low.pop();
         }


        if (lsize >= hsize) tp = low.top();
        else tp = high.top();
        cost[beg][i] = sumlow - tp * lsize;
        cost[beg][i] = cost[beg][i] + tp * hsize - sumhigh;
        //printf("%d,%d low%d high%d  pt %d cost%lld\n",beg,i,low.size(),high.size(),tp,cost[beg][i]);
    }
}
void Initial() {

    int i,j,k;


    for (i = 1; i < n; ++i) {

        while (!low.empty()) low.pop();
        while (!high.empty()) high.pop();
        Calculate(i);
    }
    for (i = 0; i <= n; ++i)
        for (j = 0; j <= m; ++j)
            dp[i][j] = INF;
}


int main()
{
    int i,j,k;


    while (scanf("%d%d",&n,&m),n) {

        for (i = 1; i <= n; ++i)
            scanf("%d",&arr[i]);
        Initial();


        dp[0][0] = 0;
        for (j = 0; j < m; ++j)
           for (k = 0; k < n; ++k)
               if (dp[k][j] != INF)
                    for (i = k + 1; i <= n ; ++i)
                        dp[i][j+1] = min(dp[i][j+1],dp[k][j]+cost[k+1][i]);


        printf("%d\n",dp[n][m]);
    }
}
/*
7 2
6 5 4 3 2 1 7
5 5
5 5 5 5 5
5 2
-5 -4 -3 -2 -1
7 3
6 5 4 3 2 1 7
*/

1007代码:

#include <stdio.h>
#include <string.h>
#define MIN 1100
#define MAX 1000000


struct node {

    int v,kind;
    node *next;
}tree[MAX*2],*head[MIN];
int ptr,n,m,k,flag;
int fa[MIN],rank[MIN];


void Initial() {

    ptr = flag = 0;
    memset(head,NULL,sizeof(head));
    for (int i = 1; i <= n; ++i)
        fa[i] = i,rank[i] = 1;
}
void AddEdge(int x,int y,int kind) {

    tree[ptr].v = y,tree[ptr].kind = kind;
    tree[ptr].next = head[x],head[x] = &tree[ptr++];
}
int Getfa(int x) {

    int r = x,q;
    while (r != fa[r]) r = fa[r];
    while (r != x) {

        q = fa[x];
        fa[x] = r;
        x = q;
    }
    return r;
}
int Solve() {

    int i,j;
    int mmax = 0,mmin = 0;
    node *p;


    for (i = 1; i <= n; ++i) {

        p = head[i];
        while (p != NULL) {

            if (p->kind) {

                int u = i,v = p->v;
                int ta = Getfa(u),tb = Getfa(v);
                if (ta != tb)
                    fa[ta] = tb,mmax++;
            }
            p = p->next;
        }
    }


    for (i = 1; i <= n; ++i)
        fa[i] = i,rank[i] = 1;
    for (i = 1; i <= n; ++i) {

        p = head[i];
        while (p != NULL) {

            if (p->kind == 0) {

                int u = i,v = p->v;
                int ta = Getfa(u),tb = Getfa(v);
                if (ta != tb)
                    fa[ta] = tb,mmin++;
            }
            p = p->next;
        }
    }

    
    mmin = n - 1 - mmin;
    //printf("%d %d\n",mmin,mmax);
    if (mmin <= k && mmax >= k) return 1;
    else return 0;
}


int main()
{
    int i,j;


    while (scanf("%d%d%d",&n,&m,&k),n+m+k) {

        Initial();
        for (i = 1; i <= m; ++i) {

            int b,c;
            char s[2];
            scanf("%s %d%d",s,&b,&c);
            AddEdge(b,c,s[0]=='B');
            AddEdge(c,b,s[0]=='B');
        }


        int flag = Solve();
        printf("%d\n",flag);
    }
}

未完,待续...

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值