有些题目可能没做，如计算几何、恶心模拟。

高级打字机

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

早苗入手了最新的高级打字机。最新款自然有着与以往不同的功能，那就是它具备撤销功能，厉害吧。

请为这种高级打字机设计一个程序，支持如下3种操作：

1.T x：在文章末尾打下一个小写字母x。(type操作)

2.U x：撤销最后的x次修改操作。（Undo操作）

（注意Query操作并不算修改操作）

3.Q x：询问当前文章中第x个字母并输出。（Query操作）

文章一开始可以视为空串。

输入

第1行：一个整数n，表示操作数量。

以下n行，每行一个命令。保证输入的命令合法。

输出

每行输出一个字母，表示Query操作的答案。

输入示例

7

T a

T b

T c

Q 2

U 2

T c

Q 2

输出示例

b

c

其他说明

对于20%的数据 n<=200;

对于50%的数据 n<=100000；保证Undo操作不会撤销Undo操作。

<高级挑战>

对于100%的数据 n<=100000；Undo操作可以撤销Undo操作。

<IOI挑战>

必须使用在线算法完成该题。

用可持久化线段树来维护数组就可以了

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=1000010;const int maxnode=2000010;char S[maxnode];int ToT,root[maxn],len[maxn],ls[maxnode],rs[maxnode];void update(int& y,int x,int l,int r,int pos,char c) { y=++ToT; if(l==r) S[y]=c; else { int mid=l+r>>1;ls[y]=ls[x];rs[y]=rs[x]; if(pos<=mid) update(ls[y],ls[x],l,mid,pos,c); else update(rs[y],rs[x],mid+1,r,pos,c); }}void query(int y,int l,int r,int pos) { if(l==r) printf("%c\n",S[y]); else { int mid=l+r>>1; if(pos<=mid) query(ls[y],l,mid,pos); else query(rs[y],mid+1,r,pos); }}int main() { int n=read(),vs=0; while(n--) { char cmd[2],s[2]; scanf("%s",cmd); if(cmd[0]=='T') { scanf("%s",s);vs++;len[vs]=len[vs-1]+1; update(root[vs],root[vs-1],1,100000,len[vs],s[0]); } else if(cmd[0]=='U') { int x=read();vs++; root[vs]=root[vs-x-1]; len[vs]=len[vs-x-1]; } else query(root[vs],1,100000,read()); } return 0;}
           
          
         
        
       
      

 

 

不等数列

难度级别：C； 运行时间限制：1000ms； 运行空间限制：256000KB； 代码长度限制：2000000B

试题描述

将1到n任意排列，然后在排列的每两个数之间根据他们的大小关系插入“>”和“<”。问在所有排列中，有多少个排列恰好有k个“<”。答案对2012取模。

输入

第一行2个整数n,k。

输出

一个整数表示答案。

输入示例

5 2

输出示例

66

其他说明

对于30%的数据：n <= 10

对于100%的数据：k < n <= 1000， 

 设f[n][k]表示1到n的排列，有k个<号的方案。

考虑将n+1插入n+1个位置中，如果插入到<之间不会增加<的数量，如果插入到>之间增加1个<。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int mod=2012;const int maxn=1010;long long f[maxn][maxn];int main() { int n=read(),k=read(); f[1][0]=1; rep(i,1,n-1) rep(j,0,i-1) { f[i][j]%=mod; f[i+1][j]+=(j+1)*f[i][j]; f[i+1][j+1]+=(i-j)*f[i][j]; } printf("%lld\n",f[n][k]%mod); return 0;}
           
          
         
        
       
      

 

经营与开发

难度级别：C； 运行时间限制：1000ms； 运行空间限制：256000KB； 代码长度限制：2000000B

试题描述

4X概念体系，是指在PC战略游戏中一种相当普及和成熟的系统概念，得名自4个同样以“EX”为开头的英语单词。

eXplore（探索）

eXpand（拓张与发展）

eXploit（经营与开发）

eXterminate（征服）

——维基百科

 

今次我们着重考虑exploit部分，并将其模型简化：

你驾驶着一台带有钻头（初始能力值w）的飞船，按既定路线依次飞过n个星球。

 

星球笼统的分为2类：资源型和维修型。（p为钻头当前能力值）

1.资源型：含矿物质量a[i]，若选择开采，则得到a[i]*p的金钱，之后钻头损耗k%，即p=p*(1-0.01k)

2.维修型：维护费用b[i]，若选择维修，则支付b[i]*p的金钱，之后钻头修复c%，即p=p*(1+0.01c)

    注：维修后钻头的能力值可以超过初始值（你可以认为是翻修+升级）

 

请作为舰长的你仔细抉择以最大化收入。

输入

第一行4个整数n,k,c,w。

以下n行，每行2个整数type,x。

type为1则代表其为资源型星球，x为其矿物质含量a[i]；

type为2则代表其为维修型星球，x为其维护费用b[i]；

输出

一个实数（保留2位小数），表示最大的收入。

输入示例

5 50 50 10

1 10

1 20

2 10

2 20

1 30

输出示例

375.00

其他说明

对于30%的数据 n<=100

另有20%的数据 n<=1000；k=100

对于100%的数据 n<=100000; 0<=k,c,w,a[i],b[i]<=100；保证答案不超过10^9

设f[i]表示后i个星球最多获益值且必须采集第i个星球，DP维护一个最大值即可。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=100010;int n,k,c,w,t[maxn],A[maxn];double f[maxn],C,K,mx;int main() { n=read();k=read();c=read();w=read(); rep(i,1,n) t[i]=read(),A[i]=read(); C=(100+c)/100.0;K=(100-k)/100.0; dwn(i,n,1) { if(t[i]==1) f[i]=mx*K+w*A[i]; else f[i]=mx*C-w*A[i]; mx=max(mx,f[i]); } printf("%.2lf\n",mx); return 0;}
           
          
         
        
       
      

 

机器人

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

早苗入手了最新的Gundam模型。最新款自然有着与以往不同的功能，那就是它能够自动行走，厉害吧。

早苗的新模型可以按照输入的命令进行移动，命令包括‘E’、‘S’、‘W’、‘N’四种，分别对应东南西北。执行某个命令时，它会向对应方向移动一个单位。作为新型机器人，它可以执行命令串。对于输入的命令串，每一秒它会按命令行动一次。执行完命令串的最后一个命令后，会自动从头开始循环。在0时刻时机器人位于（0,0）。求T秒后机器人所在位置坐标。

输入

第1行：一个字符串，表示早苗输入的命令串，保证至少有1个命令

第2行：一个正整数T

输出

2个整数，表示T秒时，机器人的坐标。

输入示例

NSWWNSNEEWN

12

输出示例

-1 3

其他说明

【数据范围】

对于60%的数据 T<=500,000 且命令串长度<=5,000

对于100%的数据 T<=2,000,000,000 且命令串长度<=5,000

【注意】

向东移动，坐标改变改变为(X+1,Y); 

向南移动，坐标改变改变为(X,Y-1); 

向西移动，坐标改变改变为(X-1,Y); 

向北移动，坐标改变改变为(X,Y+1); 

这个就不说了233

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=5010;char s[maxn];int main() { scanf("%s",s);int n=strlen(s); int T=read(); int x0=0,y0=0; rep(i,0,n-1) { if(s[i]=='W') x0--; if(s[i]=='E') x0++; if(s[i]=='N') y0++; if(s[i]=='S') y0--; } int x=x0*(T/n),y=y0*(T/n); rep(i,0,(T%n)-1) { if(s[i]=='W') x--; if(s[i]=='E') x++; if(s[i]=='N') y++; if(s[i]=='S') y--; } printf("%d %d\n",x,y); return 0;}
           
          
         
        
       
      

 

数列

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

a[1]=a[2]=a[3]=1

a[x]=a[x-3]+a[x-1]  (x>3)

求a数列的第n项对1000000007（10^9+7）取余的值。

 

输入

第一行一个整数T，表示询问个数。

以下T行，每行一个正整数n。

输出

每行输出一个非负整数表示答案。

输入示例

3

6

8

10

输出示例

4

9

19

其他说明

对于30%的数据 n<=100；

对于60%的数据 n<=2*10^7；

对于100%的数据 T<=100，n<=2*10^9；

直接矩阵快速幂即可

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int MOD=1000000007;typedef long long ll;struct Matrix { ll A[3][3]; Matrix() {memset(A,0,sizeof(A));} Matrix operator * (const Matrix& b) const { Matrix c; rep(i,0,2) rep(j,0,2) { c.A[i][j]=0; rep(k,0,2) c.A[i][j]+=A[i][k]*b.A[k][j]; c.A[i][j]%=MOD; } return c; } void print() { rep(i,0,2) { rep(j,0,2) printf("%lld ",A[i][j]); puts(""); } }};void pow(Matrix& ans,int n) { Matrix t=ans;n--; while(n) { if(n&1) ans=ans*t; t=t*t;n>>=1; }}ll solve(int n) { if(n<=3) return 1; Matrix t; t.A[0][0]=1;t.A[0][2]=1; t.A[1][0]=1;t.A[2][1]=1; pow(t,n-3); return (t.A[0][0]+t.A[1][0]+t.A[2][0])%MOD;}int main() { int T=read(); while(T--) printf("%lld\n",solve(read())); return 0;}
           
          
         
        
       
      

 

虫洞

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

N个虫洞，M条单向跃迁路径。从一个虫洞沿跃迁路径到另一个虫洞需要消耗一定量的燃料和1单位时间。虫洞有白洞和黑洞之分。设一条跃迁路径两端的虫洞质量差为delta。

1.从白洞跃迁到黑洞，消耗的燃料值减少delta，若该条路径消耗的燃料值变为负数的话，取为0。

2.从黑洞跃迁到白洞，消耗的燃料值增加delta。

3.路径两端均为黑洞或白洞，消耗的燃料值不变化。

作为压轴题，自然不会是如此简单的最短路问题，所以每过1单位时间黑洞变为白洞，白洞变为黑洞。在飞行过程中，可以选择在一个虫洞停留1个单位时间，如果当前为白洞，则不消耗燃料，否则消耗s[i]的燃料。现在请你求出从虫洞1到N最少的燃料消耗，保证一定存在1到N的路线。

输入

第1行：2个正整数N,M

第2行：N个整数，第i个为0表示虫洞i开始时为白洞，1表示黑洞。

第3行：N个整数，第i个数表示虫洞i的质量w[i]。

第4行：N个整数，第i个数表示在虫洞i停留消耗的燃料s[i]。

第5..M+4行：每行3个整数，u,v,k，表示在没有影响的情况下，从虫洞u到虫洞v需要消耗燃料k。

输出

一个整数，表示最少的燃料消耗。

输入示例

4 5

1 0 1 0

10 10 100 10

5 20 15 10

1 2 30

2 3 40

1 3 20

1 4 200

3 4 200

输出示例

130

其他说明

【数据范围】

对于30%的数据： 1<=N<=100,1<=M<=500

对于60%的数据： 1<=N<=1000,1<=M<=5000

对于100%的数据： 1<=N<=5000，1<=M<=30000

其中20%的数据为1<=N<=3000的链

1<=u,v<=N, 1<=k,w[i],s[i]<=200

【样例说明】

按照1->3->4的路线。

将一个点拆成两个点，分别是奇数秒和偶数秒第i个点，然后各种连边最短路即可。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=10010;const int maxm=1000010;const int INF=1000000000;struct Dijkstra { int n,m,first[maxn],next[maxm]; struct Edge { int to,dist;}edges[maxm]; int done[maxn],d[maxn]; struct HeapNode { int u,d; bool operator < (const HeapNode& ths) const { return d>ths.d;} }; void init(int n) { this->n=n;m=0; memset(first,0,sizeof(first)); } void AddEdge(int u,int v,int w) { //printf("%d --> %d %d\n",u,v,w); edges[++m]=(Edge){v,w};next[m]=first[u];first[u]=m; } void solve(int s) { rep(i,1,n) d[i]=INF,done[i]=0; priority_queue
            
              Q; Q.push((HeapNode){s,0});d[s]=0; while(!Q.empty()) { int x=Q.top().u;Q.pop(); if(done[x]) continue;done[x]=1; ren { Edge& e=edges[i]; if(d[e.to]>d[x]+e.dist) { d[e.to]=d[x]+e.dist; Q.push((HeapNode){e.to,d[e.to]}); } } } }}sol;//1---n 0 n+1----2*n 1int n,m,type[maxn],w[maxn],s[maxn],u[maxn*3],v[maxn*3],k[maxn*3];int first[maxn],next[maxn*3];void Add(int x,int v) { next[v]=first[x];first[x]=v;}int main() { n=read();m=read(); rep(i,1,n) type[i]=read(); rep(i,1,n) w[i]=read(); rep(i,1,n) s[i]=read(); rep(i,1,m) u[i]=read(),v[i]=read(),k[i]=read(),Add(u[i],i); sol.init(n*2); rep(x,1,n*2) { if(x<=n) sol.AddEdge(x,x+n,type[x]*s[x]); else sol.AddEdge(x,x-n,(type[x-n]^1)*s[x-n]); if(x<=n) ren { if(type[x]==type[v[i]]) sol.AddEdge(x,v[i]+n,k[i]); else if(type[x]) sol.AddEdge(x,v[i]+n,k[i]+abs(w[x]-w[v[i]])); else sol.AddEdge(x,v[i]+n,max(0,k[i]-abs(w[x]-w[v[i]]))); } else for(int i=first[x-n];i;i=next[i]) { if(type[x-n]==type[v[i]]) sol.AddEdge(x,v[i],k[i]); else if(type[x-n]) sol.AddEdge(x,v[i],max(0,k[i]-abs(w[x-n]-w[v[i]]))); else sol.AddEdge(x,v[i],k[i]+abs(w[x-n]-w[v[i]])); } } sol.solve(1); printf("%d\n",min(sol.d[n],sol.d[n*2])); return 0;}
            
           
          
         
        
       
      

 

 

双色球

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

机房来了新一届的学弟学妹，邪恶的chenzeyu97发现一位学弟与他同名，于是他当起了善良的学长233

“来来来，学弟，我考你道水题检验一下你的水平……”

一个栈内初始有n个红色和蓝色的小球，请你按照以下规则进行操作

1.      只要栈顶的小球是红色的，将其取出，直到栈顶的球是蓝色

2.      然后将栈顶的蓝球变成红色

3.      最后放入若干个蓝球直到栈中的球数为n

以上3步骤为一次操作

如栈中都是红色球，则操作停止，请问几次操作后停止

chenzeyu97出完题发现他自己不能AC所以想请你帮忙

输入

第一行为一个整数n，表示栈的容量为n

第二行为一个字符串，第i个字符表示自顶向下的第i个球的颜色，R代表红色，B代表蓝色

输出

一个整数表示操作数

输入示例

10

BRRRRRRRRR

输出示例

1

其他说明

50%的数据，1<=n<=20

100%的数据，1<=n<=50

发现将RRRRRR--B变成RRRRRR--R需要2^n次操作。

             n个R

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}char s[55];int main() { int n=read(); long long ans=0; scanf("%s",s); rep(i,0,n-1) ans|=(1ll<
           
          
         
        
       
      

 

魔方

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

ccy（ndsf）觉得手动复原魔方太慢了，所以他要借助计算机。

ccy（ndsf）家的魔方都是3*3*3的三阶魔方，大家应该都见过。

（3的“顺时针”改为“逆时针”，即3 4以图为准。）

ccy（ndfs）从网上搜了一篇攻略，并找人翻译成了他自己会做的方法。现在告诉你他的魔方情况，以及他从网上搜到的攻略，请你求出最后魔方变成什么样子。

输入数据观察角度说明图：

 

输入

　第一行，一串数字，表示从网上搜到的攻略。

　　　下面6*3行，每行3个数字，每三行表示魔方一个面的情况，六个面的顺序是前、后、左、右、上、下。

输出

　6*3行，表示处理后的魔方，形式同输入。

输入示例

23

121

221

111

123

321

111

123

321

132

132

231

132

121

112

233

332

111

333

输出示例

123

222

113

212

121

211

323

321

132

121

333

121

131

213

112

331

111

331

其他说明

40%的数据，攻略的长度小于5且仅有4种操作的其中一种

100%的数据，攻略的长度小于100

恩，问LZJ吧。

#include
      
       #include
       
        #define front 0#define back 1#define left 2#define right 3#define up 4#define down 5using namespace std;char order[110];struct Cube{    int block[6][3][3];    void Rotate(int n)    {        Cube b=*this;        if(n==1)        {            b.block[up][2][2]=block[front][2][2];            b.block[up][2][1]=block[front][2][1];            b.block[up][2][0]=block[front][2][0];            b.block[front][2][2]=block[down][2][2];            b.block[front][2][1]=block[down][2][1];            b.block[front][2][0]=block[down][2][0];            b.block[back][0][0]=block[up][2][2];            b.block[back][0][1]=block[up][2][1];            b.block[back][0][2]=block[up][2][0];            b.block[down][2][0]=block[back][0][2];            b.block[down][2][1]=block[back][0][1];            b.block[down][2][2]=block[back][0][0];            b.block[right][2][0]=block[right][0][0];            b.block[right][0][0]=block[right][0][2];            b.block[right][0][2]=block[right][2][2];            b.block[right][2][2]=block[right][2][0];            b.block[right][1][0]=block[right][0][1];            b.block[right][2][1]=block[right][1][0];            b.block[right][1][2]=block[right][2][1];            b.block[right][0][1]=block[right][1][2];        }        else if(n==2)        {            b.Rotate(1);            b.Rotate(1);            b.Rotate(1);        }        else if(n==3)        {            b.block[left][2][0]=block[back][2][0];            b.block[left][1][0]=block[back][1][0];            b.block[left][0][0]=block[back][0][0];            b.block[back][2][0]=block[right][2][0];            b.block[back][1][0]=block[right][1][0];            b.block[back][0][0]=block[right][0][0];            b.block[right][2][0]=block[front][2][0];            b.block[right][1][0]=block[front][1][0];            b.block[right][0][0]=block[front][0][0];            b.block[front][2][0]=block[left][2][0];            b.block[front][1][0]=block[left][1][0];            b.block[front][0][0]=block[left][0][0];            b.block[up][0][0]=block[up][2][0];            b.block[up][0][2]=block[up][0][0];            b.block[up][2][2]=block[up][0][2];            b.block[up][2][0]=block[up][2][2];            b.block[up][0][1]=block[up][1][0];            b.block[up][1][0]=block[up][2][1];            b.block[up][2][1]=block[up][1][2];            b.block[up][1][2]=block[up][0][1];        }        else if(n==4)        {            b.Rotate(3);            b.Rotate(3);            b.Rotate(3);        }        *this=b;    }}C;int main(){    int c;    scanf("%s",order);    for(int i=0;i<6;i++)        for(int x=0;x<3;x++)            for(int y=0;y<3;y++)            {                c=getchar();                while(!isdigit(c))c=getchar();                C.block[i][y][x]=c-'0';            }    for(int i=0,len=strlen(order);i
       
      

 

czy的后宫

难度级别：C； 运行时间限制：3000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

czy要妥善安排他的后宫，他想在机房摆一群妹子，一共有n个位置排成一排，每个位置可以摆妹子也可以不摆妹子。有些类型妹子如果摆在相邻的位置（隔着一个空的位置不算相邻），就不好看了。假定每种妹子数量无限，求摆妹子的方案数。

输入

输入有m+1行，第一行有两个用空格隔开的正整数n、m，m表示妹子的种类数。接下来的m行，每行有m个字符1或0,若第i行第j列为1，则表示第i种妹子第j种妹子不能排在相邻的位置，输入保证对称。（提示：同一种妹子可能不能排在相邻位置）。

输出

输出只有一个整数，为方案数（这个数字可能很大，请输出方案数除以1000000007的余数。

输入示例

2 2

01

10

输出示例

7

其他说明

【样例说明】

七种方案为(空，空)、（空，1）、（1、空）、（2、空）、（空、2）、（1,1）、（2,2）。

【数据范围】

20%的数据，1＜n≤5,0＜m≤10。

60%的数据，1＜n≤200,0＜m≤100。

100%的数据，1＜n≤1000000000，0＜m≤100。

注：此题时限1.5s是因为本评测机跑太慢，大家正常做

但写的太丑可能T一俩个点

将空格视为一种新的MZ，设f[i][k]表示已放了i个MZ，最后一个MZ种类为k的方案数，然后用矩阵快速幂优化。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=110;const int MOD=1000000007;typedef long long ll;int n,m;struct Matrix { ll A[maxn][maxn]; Matrix() {memset(A,0,sizeof(A));} Matrix operator * (const Matrix& b) const { Matrix c; rep(i,1,m) rep(j,1,m) { c.A[i][j]=0; rep(k,1,m) (c.A[i][j]+=A[i][k]*b.A[k][j])%=MOD; } return c; } void print() { rep(i,1,m) { rep(j,1,m) printf("%lld ",A[i][j]); puts(""); } }};void pow(Matrix& ans,int n) { Matrix t=ans;n--; while(n) { if(n&1) ans=ans*t; t=t*t;n>>=1; }}char s[maxn];int main() { Matrix ans; n=read();m=read(); rep(i,1,m) { scanf("%s",s+1); rep(j,1,m) ans.A[i][j]=(s[j]=='0'); } m++;rep(i,1,m) ans.A[i][m]=ans.A[m][i]=1; pow(ans,n); ll res=0; rep(i,1,m) (res+=ans.A[m][i])%=MOD; printf("%lld\n",res); return 0;}
           
          
         
        
       
      

 

护花

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

约翰留下他的N(N<=100000)只奶牛上山采木．他离开的时候，她们像往常一样悠闲地在草场里吃草．可是，当他回来的时候，他看到了一幕惨剧：牛们正躲在他的花园里，啃食着他心爱的美丽花朵！为了使接下来花朵的损失最小，约翰赶紧采取行动，把牛们送回牛棚． 牛们从1到N编号．第i只牛所在的位置距离牛棚Ti(1≤Ti≤2000000)分钟的路程，而在约翰开始送她回牛棚之前，她每分钟会啃食Di(1≤Di≤100)朵鲜花．无论多么努力，约翰一次只能送一只牛回棚．而运送第第i只牛事实上需要2Ti分钟，因为来回都需要时间．    写一个程序来决定约翰运送奶牛的顺序，使最终被吞食的花朵数量最小．

输入

第1行输入N，之后N行每行输入两个整数Ti和Di

输出

一个整数，表示最小数量的花朵被吞食

输入示例

6

3 1

2 5

2 3

3 2

4 1

1 6

输出示例

86

其他说明

【样例解释】

 约翰用6，2，3，4，1，5的顺序来运送他的奶牛 

与国王游戏类似，用相邻交换法自己证吧。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}typedef long long ll;const int maxn=100010;struct Cow { int t,d; bool operator < (const Cow& ths) const { return d*ths.t>ths.d*t;}}A[maxn];int n;ll ans,all;int main() { n=read(); rep(i,1,n) A[i].t=read(),A[i].d=read(),all+=A[i].d; sort(A+1,A+n+1); rep(i,1,n) { all-=A[i].d; ans+=all*A[i].t; } printf("%lld\n",ans*2); return 0;}
           
          
         
        
       
      

 

修剪草坪

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

在一年前赢得了小镇的最佳草坪比赛后，FJ变得很懒，再也没有修剪过草坪。现在，

新一轮的最佳草坪比赛又开始了，FJ希望能够再次夺冠。

然而，FJ的草坪非常脏乱，因此，FJ只能够让他的奶牛来完成这项工作。FJ有N

(1 <= N <= 100,000)只排成一排的奶牛，编号为1...N。每只奶牛的效率是不同的，

奶牛i的效率为E_i(0 <= E_i <= 1,000,000,000)。

靠近的奶牛们很熟悉，因此，如果FJ安排超过K（1<=K<=N）只连续的奶牛，那么，这些奶牛就会罢工

去开派对:)。因此，现在FJ需要你的帮助，计算FJ可以得到的最大效率，并且该方案中

没有连续的超过K只奶牛。

输入

* 第一行：空格隔开的两个整数N和K

* 第二到N+1行：第i+1行有一个整数E_i

输出

* 第一行：一个值，表示FJ可以得到的最大的效率值。

输入示例

5 2

1

2

3

4

5

输出示例

12

其他说明

输入解释：

FJ有5只奶牛，他们的效率为1，2，3，4，5。他们希望选取效率总和最大的奶牛，但是

他不能选取超过2只连续的奶牛

FJ可以选择出了第三只以外的其他奶牛，总的效率为1+2+4+5=12。

设f[i]表示前i头奶牛最少浪费的效率和，转移时用个单调队列即可。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #define rep(s,t) for(int i=s;i<=t;i++)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() {    int x=0,f=1;char c=getchar();    for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;    for(;isdigit(c);c=getchar()) x=x*10+c-'0';    return x*f;}const int maxn=100010;int n,k,A[maxn];int q[maxn],l,r;long long f[maxn],ans,mn=1ll<<50;int main() {    n=read();k=read();    rep(1,n) A[i]=read(),ans+=A[i];    int l=1,r=0;    rep(1,n) {        while(l<=r&&f[q[r]]>=f[i-1]) r--;q[++r]=i-1;        while(l<=r&&q[l]
          
         
        
       
      

 

虫洞

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

John在他的农场中闲逛时发现了许多虫洞。虫洞可以看作一条十分奇特的有向边，并可以使你返回到过去的一个时刻（相对你进入虫洞之前）。John的每个农场有M条小路（无向边）连接着N （从1..N标号）块地，并有W个虫洞（有向边）。其中1<=N<=500,1<=M<=2500,1<=W<=200。 现在John想借助这些虫洞来回到过去（出发时刻之前），请你告诉他能办到吗。 John将向你提供F(1<=F<=5)个农场的地图。没有小路会耗费你超过10000秒的时间，当然也没有虫洞回帮你回到超过10000秒以前。

输入

* Line 1: 一个整数 F, 表示农场个数。

* Line 1 of each farm: 三个整数 N, M, W。

* Lines 2..M+1 of each farm: 三个数（S, E, T）。表示在标号为S的地与标号为E的地中间有一条用时T秒的小路。

* Lines M+2..M+W+1 of each farm: 三个数（S, E, T）。表示在标号为S的地与标号为E的地中间有一条可以使John到达T秒前的虫洞。

输出

* Lines 1..F: 如果John能在这个农场实现他的目标，输出"YES",否则输出"NO"。

输入示例

2

3 3 1

1 2 2

1 3 4

2 3 1

3 1 3

3 2 1

1 2 3

2 3 4

3 1 8

输出示例

NO

YES

判负环用DFS版的SPFA比较优美。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=510;const int maxm=7010;struct SPFA { int n,m,ok,first[maxn],next[maxm]; void init(int n) { this->n=n;m=0;ok=0; fill(first+1,first+n+1,0); } struct Edge { int to,dist;}edges[maxm]; void AddEdge(int u,int v,int w) { edges[++m]=(Edge){v,w};next[m]=first[u];first[u]=m; } int vis[maxn],d[maxn]; int spfa(int x) { vis[x]=1; ren { Edge& e=edges[i]; if(d[e.to]>d[x]+e.dist) { if(vis[e.to]) return 1; else { d[e.to]=d[x]+e.dist; if(spfa(e.to)) return 1; } } } vis[x]=0; return 0; } int solve() { rep(i,1,n) d[i]=vis[i]=0; rep(i,1,n) if(spfa(i)) return 1; return 0; }}sol;int main() { int T=read(); while(T--) { int n=read(),m=read(),w=read(),a,b,c; sol.init(n); rep(i,1,m) { a=read();b=read();c=read(); sol.AddEdge(a,b,c);sol.AddEdge(b,a,c); } rep(i,1,w) { a=read();b=read();c=read(); sol.AddEdge(a,b,-c); } puts(sol.solve()?"YES":"NO"); } return 0;}
           
          
         
        
       
      

 

水灾

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

已经连续下了几天雨，却还是没有停的样子。土豪CCY刚从外地赚完1e元回来，知道不久除了自己别墅，其他的地方都将会被洪水淹没。

CCY所在的城市可以用一个N*M(N,M<=50)的地图表示，地图上有五种符号：“. * X D S”。其中“X”表示石头，水和人都不能从上面经过。“.”表示平原，CCY和洪水都可以经过。“*”表示洪水开始地方（可能有多个地方开始发生洪水）。“D”表示CCY的别墅。“S”表示CCY现在的位置。

CCY每分钟可以向相邻位置移动，而洪水将会在CCY移动之后把相邻的没有的土地淹没（从已淹没的土地）。

求CCY回到别墅的最少时间。如果聪哥回不了家，就很可能会被淹死，那么他就要膜拜黄金大神涨RP来呼叫直升飞机，所以输出“ORZ hzwer!!!”。

输入

第一行：N和M

接下来N行，每行M个字符，表示地图。

输出

若能回家，输出一个整数，表示最少回家步数；否则输出“ORZ hzwer!!!”。

输入示例

3 3

D.*

...

.S.

输出示例

3

只要CCY比洪水到达的时间早就可以走，两遍BFS即可。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=55;const int mmx[]={ 1,-1,0,0};const int mmy[]={ 0,0,-1,1};char s[maxn][maxn];int n,m,X0,Y0;queue
            
              Q;int vis[maxn][maxn],Time[maxn][maxn];void bfs1() { memset(vis,0,sizeof(vis)); rep(i,0,n-1) rep(j,0,m-1) { Time[i][j]=1<<30; if(s[i][j]=='*') { vis[i][j]=1;Q.push(i*m+j);Time[i][j]=0; } else if(s[i][j]=='S') { X0=i;Y0=j; } } while(!Q.empty()) { int x=Q.front()/m,y=Q.front()%m;Q.pop(); rep(dir,0,3) { int mx=x+mmx[dir],my=y+mmy[dir]; if(mx>=0&&mx
             
              =0&&my
              
               =0&&mx
               
                =0&&my
                
               
              
             
            
           
          
         
        
       
      

 

某种数列问题

难度级别：C； 运行时间限制：1000ms； 运行空间限制：262144KB； 代码长度限制：2000000B

试题描述

众所周知，chenzeyu97有无数的妹子(阿掉！>_<)，而且他还有很多恶趣味的问题，继上次纠结于一排妹子的排法以后，今天他有非(chi)常(bao)认(cheng)真(zhe)去研究一个奇怪的问题。有一堆他的妹子站成一排，然后对于每个妹子有一个美丽度，当然美丽度越大越好，chenzeyu97妹子很多，但是质量上不容乐观，经常出现很多美丽度为负数的妹子(喜闻乐见)，chenzeyu97希望从一排妹子里找出3队连续的妹子，使她们的美丽度和最大。注意，一个妹子不能被编入多个队伍而且一定要拿出三队，不然czy会闲着没事做~。

简单滴说就是：

给定一个数列，从中找到3个无交集的连续子数列使其和最大。

输入

第一行一个数n，表示数列长度。

接下来有n行，每行一个数，第i行为第i个数。

输出

仅有一个数，表示最大和。

输入示例

10

-1

2

3

-4

0

1

-6

-1

1

-2

输出示例

7

其他说明

【样例说明】

第一队妹子取2，3。

第二队妹子取0，1。

第三队妹子取1。

【数据范围】

请大家放心，虽然chenzeyu97妹子无数，但是这次他叫来的个数n是有限的。=v=

对于30%的数据，妹子数不大于200。

对于60%的数据，妹子数不大于2000。

对于100%的数据，妹子数不大于1000000。

而且，由于chenzeyu97没有CCR那样的影响力，所以他的妹子选完的最大美丽度和不超过maxlongint。(注:CCR随便选就爆long long，因为他是把妹狂魔=V=)。

枚举最后一段的起点，那么只要算前缀选出2段即可，具体看代码。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}typedef long long ll;const int maxn=1000010;const ll INF=1ll<<60;int n;ll A[maxn],S[maxn];ll f[maxn],g[maxn];void getf() { ll mn=0;f[0]=-INF; rep(i,1,n) f[i]=max(f[i-1],S[i]-mn),mn=min(mn,S[i]);}void getg() { ll mx=0;g[0]=-INF; rep(i,1,n) { if(i!=1) g[i]=max(g[i-1],S[i]+mx); else g[i]=-INF; mx=max(mx,(f[i]-S[i])); //printf("%lld ",g[i]); } //puts("");}void getans() { ll mx=S[n],ans=-INF; dwn(i,n,1) { ans=max(ans,g[i-1]+mx-S[i-1]); mx=max(mx,S[i-1]); } printf("%lld\n",ans);}int main() { n=read(); rep(i,1,n) A[i]=read(),S[i]=S[i-1]+A[i]; getf();getg();getans(); return 0;}/*4-1 -2 1 -3*/
           
          
         
        
       
      

 

密码锁

难度级别：C； 运行时间限制：1000ms； 运行空间限制：262144KB； 代码长度限制：2000000B

试题描述

hzwer有一把密码锁，由N个开关组成。一开始的时候，所有开关都是关上的。当且仅当开关x1,x2,x3,...xk为开，其他开关为关时，密码锁才会打开。

他可以进行M种的操作，每种操作有一个size[i]，表示，假如他选择了第i种的操作的话，他可以任意选择连续的size[i]个格子，把它们全部取反。（注意，由于黄金大神非常的神，所以操作次数可以无限>_<）

本来这是一个无关紧要的问题，但是，黄金大神不小心他的钱丢进去了，没有的钱他哪里能逃过被chenzeyu97 NTR的命运？>_<  于是，他为了虐爆czy，也为了去泡更多的妹子，决定打开这把锁。但是他那么神的人根本不屑这种”水题”。于是，他找到了你。

你的任务很简单，求出最少需要多少步才能打开密码锁，或者如果无解的话，请输出-1。

 

输入

第1行，三个正整数N，K，M，如题目所述。

第2行，K个正整数，表示开关x1,x2,x3..xk必须为开，保证x两两不同。

第三行，M个正整数，表示size[i]，size[]可能有重复元素。

输出

输出答案，无解输出-1。

输入示例

10 8 2

1 2 3 5 6 7 8 9

3 5

输出示例

2

其他说明

对于50%的数据，1≤N≤20，1≤k≤5，1≤m≤3;

对于另外20%的数据，1≤N≤10000，1≤k≤5,1≤m≤30;

对于100%的数据，1≤N≤10000，1≤k≤10，1≤m≤100。

这应该是这几道题中质量最高的题了，本蒟蒻并没有自己想出来。

注意到题目中的是区间修改，把沿途的位置取反，这个可以看做是在模2意义下，给区间的加一操作。在我们通常的思路中，对于区间的操作，原本是要修改区间长度个的位置的情况，我们都可以通过考虑它的差分序列，使得要修改的位置个数变成2个，我们要求最少的修改，使得原序列变成全0。

所以对原序列进行差分，那么每次修改就是要你对i号位置和i+size[]模2意义下的加1。

差分后的序列中，数值为1的个数是不会超过2k个，即不会超过20个。

考虑每次对i和i+x改动的过程，如果原序列中，i号位置和i+x号位置都是0的话，我们这么改，没有任何必要。所以任意时刻，数值为1的位置个数是不会增加的，那么我们可以把每一个的1看成一个的石子，那么每次我们可以把石子往某个方向移动size[]步，如果移动之后的位置存在石子的话，就对对碰，消掉了。

因为是对对碰，石子之间的关系肯定是一个匹配的关系，我们不妨求出Dist[i][j]表示，石子i要走到石子j的位置，至少需要移动多少步，那么我们可以枚举每一个石子，然后进行一遍的bfs即可，这一部分的复杂度是O(2kmn)。

现在问题转化为有一个大小不超过20的完全图，我们想要求它的最小权最大匹配。

可以用状态压缩DP，设f[S]表示集合S的最小划分代价，每次选出最小的元素作为匹配元素之一，枚举另一个元素即可。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int INF=1<<29;int n,k,m,cnt,A[10010],num[25],e[25][25],s[110];queue
            
              Q;int vis[10010],d[10010];int f[1<<25];void bfs(int x,int id) { memset(vis,0,sizeof(vis)); Q.push(x);vis[x]=1;d[x]=0; while(!Q.empty()) { x=Q.front();Q.pop(); rep(i,1,m) { if((x+s[i])<=n&&!vis[x+s[i]]) { vis[x+s[i]]=1;d[x+s[i]]=d[x]+1;Q.push(x+s[i]); } if((x-s[i])>=1&&!vis[x-s[i]]) { vis[x-s[i]]=1;d[x-s[i]]=d[x]+1;Q.push(x-s[i]); } } } rep(i,1,cnt) { if(vis[num[i]]) e[id][i]=d[num[i]]; else e[id][i]=INF; } }int main() { n=read();k=read();m=read(); rep(i,1,k) A[read()]=1; rep(i,1,m) s[i]=read(); dwn(i,++n,1) A[i]^=A[i-1]; rep(i,1,n) if(A[i]) A[i]=++cnt,num[cnt]=i; rep(i,1,n) if(A[i]) bfs(i,A[i]); rep(S,1,(1<
            
           
          
         
        
       
      

 

抓牛

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

   农夫约翰被通知，他的一只奶牛逃逸了！所以他决定，马上出发，尽快把那只奶牛抓回来．

他们都站在数轴上．约翰在N(O≤N≤100000)处，奶牛在K(O≤K≤100000)处．约翰有两种办法移动，步行和瞬移：步行每秒种可以让约翰从x处走到x+1或x-1处；而瞬移则可让他在1秒内从x处消失，在2x处出现．然而那只逃逸的奶牛，悲剧地没有发现自己的处境多么糟糕，正站在那儿一动不动．

       那么，约翰需要多少时间抓住那只牛呢？

输入

仅有两个整数N和K

输出

最短时间

输入示例

5 17

输出示例

4

直接BFS即可。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=100010;queue
            
              Q;int d[maxn],vis[maxn];void bfs(int s,int t) { swap(s,t); vis[s]=1;Q.push(s); while(!Q.empty()) { int x=Q.front();Q.pop(); if(x==t) {printf("%d\n",d[x]);return;} if(x>=1&&!vis[x-1]) { d[x-1]=d[x]+1; vis[x-1]=1; Q.push(x-1); } if(x<=99999&&!vis[x+1]) { d[x+1]=d[x]+1; vis[x+1]=1; Q.push(x+1); } if(x*2<=100000&&!vis[x*2]) { d[x<<1]=d[x]+1; vis[x<<1]=1; Q.push(x<<1); } }}int main() { bfs(read(),read()); return 0;}
            
           
          
         
        
       
      

 

路面修整

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

FJ打算好好修一下农场中某条凹凸不平的土路。按奶牛们的要求，修好后的路面高度应当单调上升或单调下降，也就是说，高度上升与高度下降的路段不能同时出现在修好的路中。 整条路被分成了N段，N个整数A_1, ... , A_N (1 <= N <= 2,000)依次描述了每一段路的高度(0 <= A_i <= 1,000,000,000)。FJ希望找到一个恰好含N个元素的不上升或不下降序列B_1, ... , B_N，作为修过的路中每个路段的高度。由于将每一段路垫高或挖低一个单位的花费相同，修路的总支出可以表示为： |A_1 - B_1| + |A_2 - B_2| + ... + |A_N - B_N| 请你计算一下，FJ在这项工程上的最小支出是多少。FJ向你保证，这个支出不会超过2^31-1。

输入

第1行: 输入1个整数：N

第2..N+1行: 第i+1行为1个整数：A_i

输出

第1行: 输出1个正整数，表示FJ把路修成高度不上升或高度不下降的最小花费

输入示例

7

1

3

2

4

5

3

9

输出示例

3

其他说明

【样例解释】

FJ将第一个高度为3的路段的高度减少为2，将第二个高度为3的路段的高度增加到5，总花费为|2-3|+|5-3| = 3，并且各路段的高度为一个不下降序列 1,2,2,4,5,5,9。

先离散，然后DP两遍。

以算单调递增为例，设f[i][k]表示前i段路已符合要求，并且第i段路修到了k高度的最小代价。f[i][k]=min(f[i][k-1],f[i-1][j]+|Ai-Bk|)

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}typedef long long ll;const int maxn=2010;int n,A[maxn],B[maxn];ll f[maxn][maxn];int main() { n=read(); rep(i,1,n) B[i]=A[i]=read(),f[i][0]=1ll<<60; sort(B+1,B+n+1); rep(i,1,n) rep(j,1,n) f[i][j]=min(f[i][j-1],f[i-1][j]+abs(B[j]-A[i])); ll ans=f[n][n]; rep(i,1,n) rep(j,1,n) f[i][j]=min(f[i][j-1],f[i-1][j]+abs(B[n-j+1]-A[i])); printf("%lld\n",min(ans,f[n][n])); return 0;}
           
          
         
        
       
      

 

教主的魔法

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

教主最近学会了一种神奇的魔法，能够使人长高。于是他准备演示给XMYZ信息组每个英雄看。于是N个英雄们又一次聚集在了一起，这次他们排成了一列，被编号为1、2、……、N。

每个人的身高一开始都是不超过1000的正整数。教主的魔法每次可以把闭区间[L, R]（1≤L≤R≤N）内的英雄的身高全部加上一个整数W。（虽然L=R时并不符合区间的书写规范，但我们可以认为是单独增加第L（R）个英雄的身高）

CYZ、光哥和ZJQ等人不信教主的邪，于是他们有时候会问WD闭区间 [L, R] 内有多少英雄身高大于等于C，以验证教主的魔法是否真的有效。

WD巨懒，于是他把这个回答的任务交给了你。

输入

 第1行为两个整数N、Q。Q为问题数与教主的施法数总和。

 第2行有N个正整数，第i个数代表第i个英雄的身高。

 第3到第Q+2行每行有一个操作：

（1）若第一个字母为“M”，则紧接着有三个数字L、R、W。表示对闭区间 [L, R] 内所有英雄的身高加上W。

（2）若第一个字母为“A”，则紧接着有三个数字L、R、C。询问闭区间 [L, R] 内有多少英雄的身高大于等于C。

输出

 对每个“A”询问输出一行，仅含一个整数，表示闭区间 [L, R] 内身高大于等于C的英雄数。

输入示例

5 3

1 2 3 4 5

A 1 5 4

M 3 5 1

A 1 5 4

输出示例

2

3

其他说明

【输入输出样例说明】

原先5个英雄身高为1、2、3、4、5，此时[1, 5]间有2个英雄的身高大于等于4。教主施法后变为1、2、4、5、6，此时[1, 5]间有3个英雄的身高大于等于4。

【数据范围】

对30%的数据，N≤1000，Q≤1000。

对100%的数据，N≤1000000，Q≤3000，1≤W≤1000，1≤C≤1,000,000,000

分块打标记即可，经典问题。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           using namespace std; inline void read(int& x){    char ch=getchar();x=0;    while(!isdigit(ch)) ch=getchar();    while(isdigit(ch)) x=x*10+ch-'0',ch=getchar();}const int maxn=1000010;int n,m,SIZE,bl[maxn],st[maxn],en[maxn],A[maxn],B[maxn],add[maxn];void reset(int x){     for(int i=st[x];i<=en[x];i++) B[i]=A[i];     sort(B+st[x],B+en[x]+1);     /*for(int i=st[x];i<=en[x];i++) printf("%d ",B[i]);     printf("  %d\n",x);*/}int query(int x,int v){    if(B[st[x]]>v) return 0;    int L=st[x],R=en[x]+1;    while(R-L>1)    {        int M=L+R>>1;        if(B[M]<=v) L=M;        else R=M;    }    return L-st[x]+1;}int main(){    read(n);read(m);SIZE=(int)sqrt(n);    for(int i=1;i<=n;i++)     {        read(A[i]);A[i]=-A[i];        bl[i]=(i-1)/SIZE+1;        if(!st[bl[i]]) st[bl[i]]=i;        en[bl[i]]=i;    }    for(int i=1;i<=bl[n];i++) reset(i);    int L,R,v,ans;    while(m--)    {       char tp[3];scanf("%s",tp);       read(L),read(R),read(v);v=-v;       if(tp[0]=='A')       {            ans=0;           for(int i=bl[L]+1;i
          
         
        
       
      

 

吃豆豆

难度级别：C； 运行时间限制：5000ms； 运行空间限制：262144KB； 代码长度限制：2000000B

试题描述

两个PACMAN吃豆豆。一开始的时候，PACMAN都在坐标原点的左下方，豆豆都在右上方。PACMAN走到豆豆处就会吃掉它。PACMAN行走的路线很奇怪，只能向右走或者向上走，他们行走的路线不可以相交。

         请你帮这两个PACMAN计算一下，他们两加起来最多能吃掉多少豆豆。

输入

     第一行为一个整数N，表示豆豆的数目。接下来N行，每行一对正整数Xi，Yi，表示第i个豆豆的坐标。任意两个豆豆的坐标都不会重合。

输出

     仅有一行包含一个整数，即两个PACMAN加起来最多能吃掉的豆豆数量。

输入示例

8

8  1

1  5

5  7

2  2

7  8

4  6

3  3

6  4

输出示例

7

其他说明

对于30%的数据，1<=N<=25；

对于70%的数据，1<=N<=500；

对于100%的数据，1<=N<=2000,1<=Xi ，Yi <=200000  ；

似乎有DP的做法？？？

这里是费用流的解法：可以发现两个PACMAN的路线总不会交叉的，那么问题转化为点容量为1的费用流。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=4010;const int maxm=5000010;const int INF=1000000000;struct ZKW { int n,m,s,t,ans,cost; int first[maxn],next[maxm]; struct Edge { int from,to,flow,cost;}edges[maxm]; void init(int n) { this->n=n;m=0; memset(first,-1,sizeof(first)); } void AddEdge(int u,int v,int cap,int cost) { edges[m]=(Edge){u,v,cap,cost};next[m]=first[u];first[u]=m++; edges[m]=(Edge){v,u,0,-cost};next[m]=first[v];first[v]=m++; } int d[maxn],vis[maxn],inq[maxn]; int dfs(int x,int a) { if(x==t||!a) {ans+=a*cost;return a;} vis[x]=1;int flow=0,f; ren { Edge& e=edges[i]; if(!e.cost&&e.flow&&!vis[e.to]&&(f=dfs(e.to,min(a,e.flow)))) { e.flow-=f;edges[i^1].flow+=f; flow+=f;a-=f;if(!a) break; } } return flow; } int bfs() { queue
            
              Q;rep(i,1,n) d[i]=INF; d[t]=0;inq[t]=1;Q.push(t); while(!Q.empty()) { int x=Q.front();Q.pop();inq[x]=0; //printf("%d %d\n",x,d[x]); ren { Edge& e=edges[i^1]; if(e.flow&&d[e.from]>d[x]+e.cost) { d[e.from]=d[x]+e.cost; if(!inq[e.from]) inq[e.from]=1,Q.push(e.from); } } } rep(i,0,m-1) edges[i].cost+=d[edges[i].to]-d[edges[i].from]; cost+=d[s];return d[s]!=INF; } int solve(int s,int t) { this->s=s;this->t=t;ans=cost=0; while(bfs()) do memset(vis,0,sizeof(vis));while(dfs(s,INF)); return ans; }}sol;int x[maxn],y[maxn];int main() { int n=read(),s=n*2+1,t=n*2+2,S=s+2;sol.init(2*n+3); sol.AddEdge(S,s,2,0); rep(i,1,n) x[i]=read(),y[i]=read(),sol.AddEdge(i,i+n,1,-1),sol.AddEdge(s,i,1,0),sol.AddEdge(i+n,t,1,0); rep(i,1,n) rep(j,1,n) if(i!=j&&x[i]<=x[j]&&y[i]<=y[j]) sol.AddEdge(i+n,j,1,0); printf("%d\n",-sol.solve(S,t)); return 0;}
            
           
          
         
        
       
      

雷神领域

难度级别：C； 运行时间限制：1000ms； 运行空间限制：262144KB； 代码长度限制：2000000B

试题描述

Lwins_7k+ (GYZ)，Synophia大陆首屈一指的天才魔法师，创造了一个新魔法：雷神领域。

这个魔法会首先在地面上形成正方网格魔法阵列，然后在某些位置召唤雷电轴标。注意：一个位置只能有一个雷电轴标存在。雷电轴标总出现在正方网格魔法阵列的顶点处，所以我们可以用一个有序整数对(x_i,y_i)标记它的位置。

然后，如果存在三个雷电轴标A,B,C满足：x_A=x_B, y_A=y_C，则该魔法会立即召唤一个位置在(x_C,y_B)的雷电轴标，如此往复直至不存在满足条件的雷电轴标组为止。

最后，Lwins_7k+将选择一条正方网格魔法阵列上的路径并使用自然力激活它们，这时候雷神领域的魔法强度就被定义为路径上的雷电轴标所占用的行数/列数中的较小值。但是由于Lwins_7k+才刚创造这个魔法，如果选择的路径中出现U形或往返子路径，那么会给施展魔法的魔法师带来一定危险性。所以选择的路径不应该包含U形子路径或往返子路径。（类似 |__ __ __| 这样的广义U形路径也不行）

由于自然力必须从魔法师的体内调和，所以路径的起点应该是魔法师所站的位置，即正方网格魔法阵列的左下角(0,0)点。

现在他想要询问你——以扫地为生来隐藏自己真实身份的大陆有史以来最伟大的式神制造师Lwins_***，他所能释放的雷神领域的最高魔法强度。当然，路径由你来为他选择。虽然如此，作为一个魔法天才，其实他只关心最高魔法强度而已。

 

输入

输入文件field.in第一行包含一个正整数N，表示接下来的数据个数。

接下来N行，每行包含两个正整数x_i,y_i，表示(x_i,y_i)处存在一个雷电轴标。

输出

输出文件field.out仅包含一行，为可能的最高魔法强度。

输入示例

7

1 1

1 1

1 1

1 2

1 3

2 2

3 1

输出示例

3

其他说明

对于30%的数据：N<=10, 0<x_i,y_i<=5。

对于80%的数据：0<x_i,y_i<=2000。

对于100%的数据：N<=15000, 0<x_i,y_i<=5000。

首先明确雷电轴标是有传递性的，我们可以建一个并查集维护横纵坐标的连通性，而如果横坐标x和纵坐标y连通，那么在(x,y)就有一个雷电轴标。

然后5000*5000的DP即可。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=5000;int pa[maxn*2+50],f[maxn+50][maxn+50];inline int findset(int x) { return pa[x]==x?x:pa[x]=findset(pa[x]);}int main() { int n=read(); rep(i,1,10000) pa[i]=i; rep(i,1,n) { int x=read(),y=read(); pa[findset(x)]=findset(y+maxn); } rep(i,1,maxn) rep(j,1,maxn) if(findset(i)==findset(j+maxn)) f[i][j]=f[i-1][j-1]+1; else f[i][j]=max(f[i][j-1],f[i-1][j]); printf("%d\n",f[maxn][maxn]); return 0;}
           
          
         
        
       
      

魔术球问题弱化版

难度级别：C； 运行时间限制：3000ms； 运行空间限制：262144KB； 代码长度限制：2000000B

试题描述

假设有 n 根柱子，现要按下述规则在这 n 根柱子中依次放入编号为 1，2，3，…的球。

（1）每次只能在某根柱子的最上面放球。

（2）在同一根柱子中，任何 2 个相邻球的编号之和为完全平方数。

试设计一个算法，计算出在 n 根柱子上最多能放多少个球。例如，在 4 根柱子上最多可放 11 个球。

对于给定的 n，计算在 n 根柱子上最多能放多少个球。

输入

第 1 行有 1 个正整数 n，表示柱子数。

输出

一行表示可以放的最大球数。

输入示例

4

输出示例

11

其他说明

N<=60，时限为3s；比起原题还有弱化在不用打出方案，方案太坑了

先二分答案，将i号球看成一个点，如果i<j且i+j是个完全平方数，那么i->j连一条有向边，每一条路径上的球可以放在一根柱子上。问题转化成DAG最小路径覆盖，判断答案是否<=n即可。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=5010;const int maxm=100010;struct Dinic { int n,m,s,t,clo,d[maxn],vis[maxn],cur[maxn]; int first[maxn],next[maxm]; struct Edge { int from,to,flow;}edges[maxm]; void init(int n) { this->n=n;m=0; memset(first,-1,sizeof(first)); } void AddEdge(int u,int v,int cap) { edges[m]=(Edge){u,v,cap};next[m]=first[u];first[u]=m++; edges[m]=(Edge){v,u,0};next[m]=first[v];first[v]=m++; } int DFS(int x,int a) { if(x==t||!a) return a; int flow=0,f; for(int& i=cur[x];i!=-1;i=next[i]) { Edge& e=edges[i]; if(d[e.to]==d[x]+1&&(f=DFS(e.to,min(a,e.flow)))) { e.flow-=f;edges[i^1].flow+=f; flow+=f;a-=f;if(!a) break; } } return flow; } int BFS() { queue
            
              Q;Q.push(s);vis[s]=++clo; while(!Q.empty()) { int x=Q.front();Q.pop();cur[x]=first[x]; ren { Edge& e=edges[i]; if(e.flow&&vis[e.to]!=clo) { vis[e.to]=clo; d[e.to]=d[x]+1; Q.push(e.to); } } } return vis[t]==clo; } int solve(int s,int t) { this->s=s;this->t=t;int flow=0; while(BFS()) flow+=DFS(s,1e9); return flow; }}sol;int m;int is(int x) { int c=(int)sqrt(x)+0.5; return c*c==x;}int check(int n) { int s=n*2+1,t=n*2+2;sol.init(n*2+2); rep(i,1,n) rep(j,i+1,n) if(is(i+j)) sol.AddEdge(i,j+n,1); rep(i,1,n) sol.AddEdge(s,i,1),sol.AddEdge(i+n,t,1); return sol.solve(s,t)+m>=n;}int main() { m=read(); int l=1,r=2010; while(l+1
             
              >1; if(check(mid)) l=mid; else r=mid; } printf("%d\n",l); return 0;}
             
            
           
          
         
        
       
      

 

征兵

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

一个国王，他拥有一个国家。最近他因为国库里钱太多了，闲着蛋疼要征集一只部队要保卫国家。他选定了N个女兵和M个男兵，但事实上每征集一个兵他就要花10000RMB，即使国库里钱再多也伤不起啊。他发现，某男兵和某女兵之间有某种关系（往正常方面想，一共R种关系），这种关系可以使KING少花一些钱就可以征集到兵，不过国王也知道，在征兵的时候，每一个兵只能使用一种关系来少花钱。这时国王向你求助，问他最少要花多少的钱。

输入

第一行：T，一共T组数据。

接下来T组数据，

第一行包括N,M,R

接下来的R行 包括Xi,Yi,Vi 表示如果招了第Xi个女兵，再招第Yi个男兵能省Vi元（同样表示如果招了第Yi个男兵，再招第Xi个女兵能也省Vi元）

输出

共T行，表示每组数据的最终花费是多少（因为国库里的钱只有2^31-1，所以保证最终花费在maxlongint范围内）

输入示例

2

5 5 8

4 3 6831

1 3 4583

0 0 6592

0 1 3063

3 3 4975

1 3 2049

4 2 2104

2 2 781

5 5 10

2 4 9820

3 2 6236

3 1 8864

2 4 8326

2 0 5156

2 0 1463

4 1 2439

0 4 4373

3 4 8889

2 4 3133

输出示例

71071

54223

其他说明

数据保证T<=5 ,m,n<=10000,r<=50000,Xi<=m,Yi<=n,Vi<=10000,结果<=2^31-1

《挑战程序设计竞赛》有详细解答。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=20010;const int maxm=50010;struct Edge { int u,v,w; bool operator < (const Edge& ths) const { return w>ths.w; }}e[maxm];int pa[maxn];int findset(int x) { return x==pa[x]?x:pa[x]=findset(pa[x]);}int main() { int T=read(); while(T--) { int n=read(),m=read(),r=read(),ans=0; rep(i,1,n+m) pa[i]=i; rep(i,1,r) e[i].u=read()+1,e[i].v=read()+n+1,e[i].w=read(); sort(e+1,e+r+1); rep(i,1,r) { int x=findset(e[i].u),y=findset(e[i].v); if(x!=y) ans+=e[i].w,pa[x]=y; } printf("%d\n",(n+m)*10000-ans); } return 0;}
           
          
         
        
       
      

 

 

无线通讯网

难度级别：C； 运行时间限制：1000ms； 运行空间限制：51200KB； 代码长度限制：2000000B

试题描述

国防部计划用无线网络连接若干个边防哨所。2种不同的通讯技术用来搭建无线网络；每个边防哨所都要配备无线电收发器；有一些哨所还可以增配卫星电话。

任意两个配备了一条卫星电话线路的哨所（两边都拥有卫星电话）均可以通话，无论他们相距多远。而只通过无线电收发器通话的哨所之间的距离不能超过D，这是受收发器的功率限制。收发器的功率越高，通话距离D会更远，但同时价格也会更贵。

收发器需要统一购买和安装，所以全部哨所只能选择安装一种型号的收发器。换句话说，每一对哨所之间的通话距离都是同一个D。

你的任务是确定收发器必须的最小通话距离D，使得每一对哨所之间至少有一条通话路径（直接的或者间接的）。

输入

第1行：2个整数S（1<=S<=100）和P（S<P<=500），S表示可安装的卫星电话的哨所数，P表示边防哨所的数量。

     接下里P行，每行描述一个哨所的平面坐标（x,y），以km为单位，整数，0<=x,y<=10000。

输出

第1行：1个实数D，表示无线电收发器的最小传输距离。精确到小数点后两位。

输入示例

2 4

 0 100

 0 300

 0 600

 150 750

输出示例

212.13

其他说明

对于20%的数据  P=2，S=1

对于另外20%的数据  P=4，S=2

对于100%的数据  1<=S<=100，S<P<=500

二分答案，用并查集判判即可。

#include
      
       #include
       
        #include
        
         #include
         
          #include
          
           #include
           
            #define rep(i,s,t) for(int i=s;i<=t;i++)#define dwn(i,s,t) for(int i=s;i>=t;i--)#define ren for(int i=first[x];i!=-1;i=next[i])using namespace std;inline int read() { int x=0,f=1;char c=getchar(); for(;!isdigit(c);c=getchar()) if(c=='-') f=-1; for(;isdigit(c);c=getchar()) x=x*10+c-'0'; return x*f;}const int maxn=510;int n,m,cnt,x[maxn],y[maxn];struct Edge { int u,v,w;}e[maxn*maxn>>1];int pa[maxn];int findset(int x) { return x==pa[x]?x:pa[x]=findset(pa[x]);}int check(int x) { rep(i,1,n) pa[i]=i; rep(i,1,cnt) if(e[i].w<=x) pa[findset(e[i].u)]=findset(e[i].v); int res=0; rep(i,1,n) if(i==findset(i)) res++; return res==1||res<=m;}int main() { m=read();n=read(); rep(i,1,n) x[i]=read(),y[i]=read(); rep(i,1,n) rep(j,i+1,n) e[++cnt]=(Edge){i,j,(x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j])}; int l=0,r=300000000,mid; while(l
            
             >1)) r=mid; else l=mid+1; printf("%.2lf\n",sqrt(l)); return 0;}