洛谷八月月赛Round1凄惨记

个人背景：

上午9:30放学，然后因为学校举办读书工程跟同学去书城选书，中午回来开始打比赛，下午又回老家，中间抽出一点时间调代码，回家已经8:50了

也许是7月月赛时“连蒙带骗”AK的太幸运然而因同学的id评测过判代码雷同扣100分后while(true) rp--;本次是一个凄惨.....

我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了

我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了

我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了

我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了

我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了

我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了我太弱了

－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－

题目：

http://scarlet.is-programmer.com/user_files/Scarlet/File/COOL's%20Life.pdf

－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－

过程：

第一题在去书城的公交车上就开始想，对于需要快速找到行x,y坐标都大的点，分别找x和y然后取交集怎么样？貌似还是n^2

唉，范围20*20的方格，可以每个宿舍保存在哪里的妹子啊，那样n＊400

然后，调试时出现迷之问题，用了vector二维数组，好像是溢出或是什么问题，造数据的循环卡在刚结束，怎么搞都不好

【2016-08-12在win下测试，把频繁的debug输出该少些，n=90000也出结果了，但很慢，洛谷仍re】

至少会了double取模 呵呵呵

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<cmath>
#include<cstring>
using namespace std;
typedef long double ld;
const int N=,rp=;

int n,cnt=;
double rxa,rxc,rya,ryc,rza,rzc;
double ans=;
vector<int> map[][];    //num
//int map2[21][21][N];
struct people{
    double x,y;
    int i,j;
    people(){x=;y=;i=;j=;}
    //int num;
}p[N];

double dmod(double x,int m){
    int t=x;
    int j=t/m;
    return x-j*m;
}
void init(){
    double x=,y=,z=;
    for(int i=;i<=n;i++){
        x=dmod((y*rxa+rxc),rp);
        y=dmod((z*rya+ryc),rp);
        z=dmod((x*rza+rzc),rp);
        p[i].x=dmod(x,)+dmod(y,)/+dmod(z,)/; p[i].i=p[i].x;

        x=dmod((y*rxa+rxc),rp);
        y=dmod((z*rya+ryc),rp);
        z=dmod((x*rza+rzc),rp);
        p[i].y=dmod(x,)+dmod(y,)/+dmod(z,)/; p[i].j=p[i].y;

        //if(i>89000||i%1000==0) printf("%d i %d j %d\n",i,p[i].i,p[i].j);
        map[p[i].i][p[i].j].push_back(i);
        //printf("%d %lf %lf\n",i,p[i].x,p[i].y);printf("hello");
    }
    return;
}

inline double cal(people &a,people &b){
    return fabs(a.x-b.x)+fabs(a.y-b.y);
}

void solve(){
    for(int i=;i<=n;i++){
        people &now=p[i];//if(i>89000||i%1000==0) printf(" i %d\n",i);
        for(int x=now.i+;x<;x++)
            for(int y=now.j+;y<;y++)
               for(int z=;z<map[x][y].size();z++){
                    ans+=cal(now,p[map[x][y][z]]); cnt++;
                }
    }
}

int main(){//system("pause");
    cin>>n;
    cin>>rxa>>rxc>>rya>>ryc>>rza>>rzc;
    init();

    //cout<<"p";
    solve();
    if(cnt==) printf("0.00000");
    printf("%.5f",ans/cnt);

    //cout<<"\n\n";
    //for(int i=1;i<=n;i++) printf("p  %d  %lf %lf %d %d\n",i,p[i].x,p[i].y,p[i].i,p[i].j);
}

用正解的思路写了一遍，依旧30，好像init()就一直不对

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<cmath>
#include<cstring>
using namespace std;
typedef long double ld;
const int N=,rp=;

int n,cnt=;
double rxa,rxc,rya,ryc,rza,rzc;
double ans=;

struct node{
    int num;
    double sum;
    node(int a=,double b=):num(a),sum(b){}
};
node map[][];    

struct people{
    double x,y;
    int i,j;
    people(){x=;y=;i=;j=;}
    //int num;
}p[N];

double dmod(double x,int m){
    int t=x;
    int j=t/m;
    return x-j*m;
}
void init(){
    double x=,y=,z=;
    for(int i=;i<=n;i++){
        x=dmod((y*rxa+rxc),rp);
        y=dmod((z*rya+ryc),rp);
        z=dmod((x*rza+rzc),rp);
        p[i].x=dmod(x,)+dmod(y,)/+dmod(z,)/; p[i].i=p[i].x;

        x=dmod((y*rxa+rxc),rp);
        y=dmod((z*rya+ryc),rp);
        z=dmod((x*rza+rzc),rp);
        p[i].y=dmod(x,)+dmod(y,)/+dmod(z,)/; p[i].j=p[i].y;

        map[p[i].i][p[i].j].num++; map[p[i].i][p[i].j].sum+=p[i].x-p[i].i+p[i].y-p[i].j;
       // printf("%d i %d j %d  sum %f\n",i,p[i].i,p[i].j,map[p[i].i][p[i].j].sum);
    }
    return;
}

void solve(){
    for(int i=;i<;i++)
        for(int j=;j<;j++){
            for(int ni=i+;ni<;ni++)
                for(int nj=j+;nj<;nj++){
                    node &node1=map[i][j],&node2=map[ni][nj];
                    //if(node1.num||node2.num)printf("map %d %f %d %f\n",node1.num,node1.sum,node2.num,node2.sum);
                    cnt+=node1.num*node2.num;
                    ans+=node1.num*node2.num*(ni-i+nj-j)-node2.num*node1.sum+node1.num*node2.sum;
                //    printf("ans %.5f\n",ans);
                }
        }
}
int main(){//system("pause");
    cin>>n;
    cin>>rxa>>rxc>>rya>>ryc>>rza>>rzc;
    init();
    solve();
    //printf("%f %d\n",ans,cnt);

    printf("%.5f",ans/cnt);

    //cout<<"p";

    //cout<<"\n\n";
    //for(int i=1;i<=n;i++) printf("p  %d  %lf %lf %d %d\n",i,p[i].x,p[i].y,p[i].i,p[i].j);
}

第二题，数学。

我只想到了对于每个数根号n枚举约数快速幂暴力，预处理前缀和，当l==r时改为直接算f(i)

莫名其妙给快速幂加了个记忆化，最终骗了40分

//
//  main.cpp
//  luogu8_2
//
//  Created by abc on 16/8/11.
//  Copyright © 2016年 abc. All rights reserved.
//

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
typedef long long ll;
const int MOD=1e9+,N=1e7+;
inline ll read(){
    char c=getchar(); ll x=,f=;
    while(c<''||c>''){if(c=='-')f=-;c=getchar();}
    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}
    return x;
}

ll n,d,q,l,r;
ll f[N],s[N];

ll pm[N];
inline ll powMod(ll a,ll b){
    if(pm[a]) return pm[a]; ll t=a;
    ll ans=;
    for(;b;b>>=,a=(a*a)%MOD)
        if(b&) ans=(ans*a)%MOD;
    return pm[t]=ans;
}

inline ll divisorPow(ll x){
    if(f[x]) return f[x];
    int m=sqrt(x+0.5);
    ll ans=; //
    for(int i=;i<=m;i++) if(x%i==){
        ans=(ans+powMod(i,d))%MOD;
        if(i*i!=x) ans=(ans+powMod(x/i,d))%MOD;
    }
    if(x!=) ans=(ans+powMod(x,d))%MOD;
    return f[x]=ans;
}

bool flag=;
void init(){
    if(flag) return;
    flag=;
    for(int i=;i<=n;i++){
        f[i]=divisorPow(i);s[i]=(s[i-]+f[i])%MOD;
    }
}

int main(int argc, const char * argv[]) {
    n=read();d=read();
    q=read();
    for(int i=;i<=q;i++){
        l=read();r=read();
        if(l==r) printf("%lld\n",divisorPow(l));
        else{init();printf("%lld\n",(s[r]-s[l-]+MOD)%MOD);}
    }
    //cout<<f[719]<<" "<<f[720]<<" "<<s[720];
    return ;
}

第三题，关于树，想了想，可以枚举神秘对象，然后对于他枚举所有链，利用lca统计经过他的链(具体方法见代码)，排下序取最大最小值即可。

复杂度O(n*m*logn) 要排序和用lca

然而只拿三十分，感觉n=5000,m=500也可以过啊？难道是rp