20160729noip模拟赛zld

首先显然有多少个奇数，就有多少个回文串是最优的（没有奇数时构造一个回文串

然后有了k个“核心”,把剩下的字符顺序安排到这些的两侧，最后最短的回文串长度就是答案

#include<map>
#include<stack>
#include<queue>
#include<cstdio>
#include<string>
#include<vector>
#include<cstring>
#include<complex>
#include<iostream>
#include<assert.h>
#include<algorithm>
using namespace std;
#define inf 1001001001
#define infll 1001001001001001001LL
#define FOR(i,a,n) for(int (i)=a;(i)<=(n);++(i))
#define FOR1(i,n) for(int (i)=1;(i)<=(n);++(i))
#define ll long long
#define dbg(vari) cerr<<#vari<<" = "<<(vari)<<endl
#define gmax(a,b) (a)=max((a),(b))
#define gmin(a,b) (a)=min((a),(b))
#define ios0 ios_base::sync_with_stdio(0)
#define Ri register int
#define gc getchar()
#define il inline
il int read(){
    bool f=true;
    Ri x=;char ch;
    while(!isdigit(ch=gc))if(ch=='-')f=false;
    while(isdigit(ch)){x=(x<<)+(x<<)+ch-'';ch=gc;}
    return f?x:-x;
}
#define gi read()
#define FO(x) freopen(#x".in","r",stdin),freopen(#x".out","w",stdout);
int n,a[];
int main(){
    FO(palindromic);
    int T=gi;
    while(T--){
        int cnt=,sum=;
        n=gi;
        FOR1(i,n){
            a[i]=gi;sum+=a[i];
            if(a[i]&)cnt++;
        }
//        printf("%d %d\n",cnt,sum);
        if(cnt)
        cout<<(sum-cnt)//cnt*+<<endl;
        else cout<<sum<<endl;
    }
    return ;
}

上次做这题：20160430ysy出题的时候

首先任意两点的最短路只有两种情况：曼哈顿距离，曼哈顿距离+2

那么我们考虑怎样的点对曼哈顿距离会被影响

显然是下面这样的：黑色是障碍，红色和蓝点距离加2

那么我们分行列统计这样的对数，就可以（

//fyb
//旅行
#include<cmath>
#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
typedef double ll;
ll tot;
ll gg;
int n,m;
char s[];
int row[],col[];
int fyb(int i,int j,int k){return *(i-)*(k-j);}
ll dorow(){
    ll ans=;
    for(int i=;i<n;i++)
        for(int j=i+;j<=n;j++)
            ans=ans+(m-(row[i]>))*(m-(row[j]>))*(j-i)/gg;
    for (int i=;i<=n;i++)
        if (row[i]){
            ans=ans+fyb(row[i],row[i],m)/gg;
            for(int j=i-;j>&&row[j]&&row[j]<row[j+];j--)    ans=ans + fyb(row[j],row[i],m)/gg;
            for(int j=i+;j<=n&&row[j]&&row[j]<row[j-];j++)ans=ans + fyb(row[j],row[i],m)/gg;
        }
    return ans;
}
ll docol(){
    ll ans=;
    for (int i=;i<m;i++)
        for (int j=i+;j<=m;j++)
            ans=ans+(n-(col[i]>))*(n-(col[j]>))*(j-i)/gg;
    for (int i=;i<=m;i++)
        if (col[i]){
            ans=ans+fyb(col[i],col[i],n)/gg;
            for(int j=i-;j>&&col[j]&&col[j]<col[j+];j--) ans= ans + fyb(col[j],col[i],n)/gg;
            for(int j=i+;j<=m&&col[j]&&col[j]<col[j-];j++) ans=ans + fyb(col[j],col[i],n)/gg;
        }
    return ans;
}
int main(){
    freopen("chess.in","r",stdin);
    freopen("chess.out","w",stdout);
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        tot=;memset(row,,sizeof(row));memset(col,,sizeof(col));
        for (int i=;i<=n;i++){
            scanf("%s",s);
            for (int j=;j<m;j++)
                if (s[j]=='G'){
                    tot++;
                    row[i]=j+;
                    col[j+]=i;
                    break;
                }
        }
        gg=n*m-tot;
        gg = gg * gg;
        double gg2 = dorow() + docol();
        gg2 = gg2 * 2.0;
        printf("%.4lf\n", gg2);
    }
    return ;
}

（出现了玄学的换行

设定ans[i][j]

表示当前时刻，从i-j最早什么时候到达

首先我们发现，加上一条边只会影响这条边的起点出发的边和到达终点的边

那么我们反向加边，每次更新一些ans[i][j]

查询就变得简洁很多了。

#include<map>
#include<stack>
#include<queue>
#include<cstdio>
#include<string>
#include<vector>
#include<cstring>
#include<complex>
#include<iostream>
#include<assert.h>
#include<algorithm>
using namespace std;
#define inf 1001001001
#define infll 1001001001001001001LL
#define FOR0(i,n) for(int (i)=0;(i)<(n);++(i))
#define FOR1(i,n) for(int (i)=1;(i)<=(n);++(i))
#define ll long long
#define dbg(vari) cerr<<#vari<<" = "<<(vari)<<endl
#define gmax(a,b) (a)=max((a),(b))
#define gmin(a,b) (a)=min((a),(b))
#define ios0 ios_base::sync_with_stdio(0)
#define Ri register int
#define gc getchar()
#define il inline
il int read(){
    bool f=true;
    Ri x=;char ch;
    while(!isdigit(ch=gc))if(ch=='-')f=false;
    while(isdigit(ch)){x=(x<<)+(x<<)+ch-'';ch=gc;}
    return f?x:-x;
}
#define gi read()
#define FO(x) freopen(#x".in","r",stdin),freopen(#x".out","w",stdout);
int n,m,q,u[],v[],ans[];
struct data{
    int a,b,c,d,id;
    void input(){
        a=gi;b=gi;c=gi;d=gi;
    }
}a[];
bool operator<(data a,data b){
    return a.a<b.a;
}
int dis[][];
int main(){
    FO(plan);
    n=gi;m=gi;q=gi;
    memset(dis,/,sizeof(dis));
    FOR1(i,n)    dis[i][i]=;
    FOR1(i,m)    u[i]=gi,v[i]=gi;
    FOR1(i,q)    a[i].input(),a[i].id=i;
    sort(a+,a+q+);
    int t=q;
    for(int i=m;i>=;i--){
        dis[u[i]][v[i]]=dis[v[i]][u[i]]=i;
        for(int j=;j<=n;j++)    dis[u[i]][j]=min(dis[u[i]][j],max(dis[v[i]][j],i));
        for(int j=;j<=n;j++)    dis[v[i]][j]=min(dis[v[i]][j],max(dis[u[i]][j],i));
        while(t&&a[t].a==i){
            if(dis[a[t].c][a[t].d]<=a[t].b)ans[a[t].id]=true;
            --t;
        }
    }
    for(int i=;i<=q;i++) if(ans[i]) puts("Yes"); else puts("No");
}