Noip 2013day 1 被虐报告

第一题:
刚开始的时候一下子傻了眼,居然不是白痴题?!
然后就去写T3了。后来才写了个O(n)的求循环节算法,后来才发现可以快速幂水过
Ans=(x+10^k m)%n=(x+10^k%n)%n
然后10^k%n 可以快速幂过

第二题:
看都没看据说是Greedy+bit

首先,将两个序列各自排序,然后大对大,小对小,这个贪心可以保证(ai-bi)^2最小,证明:

sigma(ai-bi)^2=sigma ai^2 + sigma bi^2 – 2sigma aibi

这时,只要sigma aibi最大就可以了,那么证明在该贪心下这个多项式最大:

若存在a>b,c>d,且ac+bd<ad+bc,则a(c-d)<b(c-d),则a<b,与a>b矛盾,所以若a>b,c>d,则ac+bd>ad+bc

将此式子进行推广:

当a1<a2<a3<…<an ,b1<b2<…<bn的情况下 sigma(ai*bi)最大,即sigma(ai-bi)^2最小。

明显,我们同时移动2个序列中的数与移动一个序列中的数是等效的,所以,我们可以保持一个序列不变,只移动另外一个序列,这样,我们就可以在排序后重定义另外一个序列的数的优先级,然后求逆序对个数,就是最小交换个数

第三题:
Kruskal求最大生成树。然后求树上路径最值,接着就什么都可以搞啦,树上倍增,lct或者树链剖分随便上就可以a了。

代码(在smart oj的非官方数据上过了,应该没问题吧):

T1(快速幂):

#include <cstdio>

#include <cstring>



using namespace std;

 

int k;

long long ans;

int n,m,x;

 

long long Exp(int y){

    if(!y)return 1;

    if(y==1)return 10%n;

    long long X=Exp(y>>1);

    if(y&1){

        return((X*X%n)*10)%n;

    }else return(X*X)%n;

}

 

int main(){

    scanf("%d%d%d%d",&n,&m,&k,&x);

    ans=Exp(k);

    ans*=m;

    ans%=n;

    ans+=x;

    ans%=n;

    printf("%I64d",ans);

    fclose(stdin),fclose(stdout);

    return 0;

}


T2(整理两个数组成一个,然后求逆序对):

#include<cstdio>
#include<algorithm>
#include<cstring>
 
using namespace std;
 
#define MAXN 100010
#define lowbit(x)(((~(x))+1)&x)
#define MAX 99999997
 
struct saver{
    int v,t;
};
saver a[MAXN],b[MAXN];
 
bool cmp(saver x,saver y){
    return x.v<y.v;
}
 
int n,r[MAXN],ans=0;
int t[MAXN];
 
void Add(int x){for(int i=x;i<=n;i+=lowbit(i)) t[i]++;}
int Sum(int x){
    int rec=0;
    for(;x;x-=lowbit(x)) rec+=t[x];
    return rec;
}
 
int main(){
    scanf("%d",&n);
    for(int i=0;i++<n;) scanf("%d",&a[i].v),a[i].t=i;
    for(int i=0;i++<n;) scanf("%d",&b[i].v),b[i].t=i;
    sort(a+1,a+n+1,cmp),sort(b+1,b+n+1,cmp);
    for(int i=0;i++<n;) r[a[i].t]=b[i].t;
    for(int i=n;i;i--) ans+=Sum(r[i]),Add(r[i]),ans%=MAX;
    printf("%d\n",ans);
    return 0;
}

T3(Kruskal+树上倍增):

#include <cstdio>
#include <algorithm>
#include <cstring>
 
using namespace std;
 
#define MAXN 10010
#define MAXM 50010
#define MAXQ 30010
#define MAXD 20
#define clear(x) memset(x,0,sizeof(x))
#define AddEdge(s,t,d) Add(s,t,d),Add(t,s,d)
#define inf 0x7fffffff
 
struct saver {
    int s,t,d;
} e[MAXM];
 
bool cmp(saver x,saver y) {
    return x.d>y.d;
}
 
int father[MAXN],n,m,q,Q[MAXQ][2];
 
int Find(int x) {
    int i,j=x;
    for (i=x;father[i];i=father[i]) ;
    while (father[j]) {
        int k=father[j];
        father[j]=i;
        j=k;
    }
    return i;
}
 
int up[MAXN][MAXD],Min[MAXN][MAXD],h[MAXN];
bool f[MAXN];
 
struct edge {
    edge *next;
    int t,d;
    edge () {
        next=NULL;
    }
} *head[MAXN];
 
void Add(int s,int t,int d) {
    edge *p=new(edge);
    p->t=t,p->d=d,p->next=head[s];
    head[s]=p;
}
 
void BuildTree(int v) {
    f[v]=false;
    for (edge *p=head[v];p;p=p->next) if (f[p->t]) {
        up[p->t][0]=v,Min[p->t][0]=p->d,h[p->t]=h[v]+1;
        BuildTree(p->t);
    }
}
 
int Move(int &v,int H) {
    int rec=inf;
    for (int i=MAXD-1;i>=0;i--) {
        if (h[up[v][i]]>=H) {
            rec=min(rec,Min[v][i]);
            v=up[v][i];
        }
    }
    return rec;
}
 
int Query(int u,int v) {
    if (Find(u)!=Find(v)) return -1;
    int rec=inf;
    if (h[u]!=h[v]) rec=h[u]>h[v]?Move(u,h[v]):Move(v,h[u]);
//    printf("%d\n",rec);
    if (u==v) return rec;
    for (int i=MAXD-1;i>=0;i--) {
        if (up[u][i]!=up[v][i]) {
            rec=min(rec,min(Min[u][i],Min[v][i]));
            u=up[u][i],v=up[v][i];
        }
    }
    rec=min(rec,min(Min[u][0],Min[v][0]));
    return rec;
}
 
int main() {
    clear(father),clear(head),clear(up);
    scanf("%d%d",&n,&m);
    for (int i=0;i<m;i++) scanf("%d%d%d",&e[i].s,&e[i].t,&e[i].d);
    sort(e,e+m,cmp);
    for (int i=0;i<m;i++) if (Find(e[i].s)!=Find(e[i].t)) {
        father[Find(e[i].s)]=Find(e[i].t);
        AddEdge(e[i].s,e[i].t,e[i].d);
    }
    memset(f,true,sizeof(f));
    for (int i=0;i++<n;) if (f[i]) h[i]=0,BuildTree(i),Min[i][0]=inf,up[i][0]=i;
    for (int i=0;i++<MAXD-1;) {
        for (int j=0;j++<n;) {
            up[j][i]=up[up[j][i-1]][i-1];
            Min[j][i]=min(Min[j][i-1],Min[up[j][i-1]][i-1]);
        }
    }
    scanf("%d",&q);
    while (q--) {
        int u,v;
        scanf("%d%d",&u,&v);
        printf("%d\n",Query(u,v));
    }
    return 0;
}

    原文作者:AmadeusChan
    原文地址: https://www.jianshu.com/p/a2bc639a279b
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞