前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >义乌中学暑假集训 2021.07.12 C

义乌中学暑假集训 2021.07.12 C

作者头像
yzxoi
发布2022-09-19 13:31:43
3550
发布2022-09-19 13:31:43
举报
文章被收录于专栏:OI

义乌中学暑假集训 2021.07.12 C

Description

有一个 2\cdot 10^9\times 2\cdot 10^9 的网格图,现要从 (x_1,y_1) 走到 (x_2,y_2),每次只能走上下左右四个方向且不能走到网格图外面。

假设当前位置为 (x,y)

  • 向南走一步(x 加一)的代价为 2xy^2+2y^2+x^2
  • 向北走一步(x 减一)的代价为 -2xy^2+2y^2+x^2
  • 向东走一步(y 加一)的代价为 2x^2y+2x^2+y^2
  • 向西走一步(y 减一)的代价为 -2x^2y+2x^2+y^2

问从 (x_1,y_1) 走到 (x_2,y_2) 的最小代价对 998244353 取模的结果是多少?

1\leq n\leq 5\times 10^4,1\leq x_1,y_1,x_2,y_2\leq 10^9

Solution

​ 的花费:

T=x_k^2y_k^2-x_1^2y_1^2+\sum_{i=1}^{k-1}x_i^2+y_i^2

我们把起点和终点有关的项 先拿掉,相当于我们要最小化:

\sum_{i=1}^kx_i^2+y_i^2

然后就可以无脑分类讨论了:

Code

代码语言:javascript
复制
#include<bits/stdc++.h>
#define I inline
#define W while
#define RI register int 
#define Cn const
#define CI Cn int& 
#define gc getchar
#define pc putchar
#define int long long
using namespace std;
I void read(int& x){RI f=1;char c=gc();x=0;W(!('0'<=c&&c<='9')) f=c^'-'?f:-1,c=gc();W('0'<=c&&c<='9') x=x*10+(c-'0'),c=gc();x*=f;}
I void write(int x){x<0&&(pc('-'),x=-x),x>=10&&(write(x/10),0),pc(x%10+'0');}
Cn int N=510,P=998244353,Inv6=(P+1)/6;
int Tt,sx,sy,tx,ty,Ans;
signed main(){
    read(Tt);W(Tt--){
        read(sx),read(sy),read(tx),read(ty);
        #define LS(x) ((x)*((x)+1)%P*(2*(x)+1)%P*Inv6%P)
        #define Sqr(x) ((x)*(x)%P)
        #define CX(sx,sy,tx,ty) (sy<ty?(LS(ty)-LS(sy)+Sqr(sy)-Sqr(ty)+Sqr(sx)*(ty-sy)%P):(LS(sy)-LS(ty)+Sqr(sx)*(sy-ty)))
        #define CY(sx,sy,tx,ty) (sx<tx?(LS(tx)-LS(sx)+Sqr(sx)-Sqr(tx)+Sqr(sy)*(tx-sx)%P):(LS(sx)-LS(tx)+Sqr(sy)*(sx-tx)))
        #define CZ(sx,sy,tx,ty) (sx<tx?(2*(LS(tx)-LS(sx))-Sqr(tx)+Sqr(sx)+2*(LS(tx)-LS(sx)-Sqr(tx)+Sqr(sx))%P):(2*(LS(sx)-LS(tx))-Sqr(sx)+Sqr(tx)+2*(LS(sx)-LS(tx))%P))
        Ans=(Sqr(tx)*Sqr(ty)%P-Sqr(sx)*Sqr(sy))%P;
        if(sx==tx){Ans+=CX(sx,sy,tx,ty);}else
        if(sy==ty){Ans+=CY(sx,sy,tx,ty);}else
        if(sx==sy&&tx==ty){Ans+=CZ(sx,sy,tx,ty);}else
        if(sx<tx){
            if(sy>=ty) (Ans+=(CX(sx,sy,sx,ty)+CY(sx,ty,tx,ty))%P)%=P;else
            if(sx>=sy&&tx>=ty) (Ans+=(sx<=ty?(CX(sx,sy,sx,sx)+CZ(sx,sx,ty,ty)+CY(ty,ty,tx,ty)):(CX(sx,sy,sx,ty)+CY(sx,ty,tx,ty)))%P)%=P;else
            if(sx< sy&&tx< ty) (Ans+=(sy<=tx?(CY(sx,sy,sy,sy)+CZ(sy,sy,tx,tx)+CX(tx,tx,tx,ty)):(CY(sx,sy,tx,sy)+CX(tx,sy,tx,ty)))%P)%=P;else
            if(sx>=sy&&tx< ty) (Ans+=(CX(sx,sy,sx,sx)+CZ(sx,sx,tx,tx)+CX(tx,tx,tx,ty))%P)%=P;else
            if(sx< sy&&tx>=ty) (Ans+=(CY(sx,sy,sy,sy)+CZ(sy,sy,ty,ty)+CY(ty,ty,tx,ty))%P)%=P;
        }else{
            if(sy<=ty) (Ans+=(CY(sx,sy,tx,sy)+CX(tx,sy,tx,ty))%P)%=P;else
            if(sx>=sy&&tx>=ty) (Ans+=(sy<=tx?(CY(sx,sy,tx,sy)+CX(tx,sy,tx,ty)):(CY(sx,sy,sy,sy)+CZ(sy,sy,tx,tx)+CX(tx,tx,tx,ty)))%P)%=P;else
            if(sx< sy&&tx< ty) (Ans+=(sx<=ty?(CX(sx,sy,sx,ty)+CY(sx,ty,tx,ty)):(CX(sx,sy,sx,sx)+CZ(sx,sx,ty,ty)+CY(ty,ty,tx,ty)))%P)%=P;else
            if(sx>=sy&&tx< ty) (Ans+=(CY(sx,sy,sy,sy)+CZ(sy,sy,ty,ty)+CY(ty,ty,tx,ty))%P)%=P;else
            if(sx< sy&&tx>=ty) (Ans+=(CX(sx,sy,sx,sx)+CZ(sx,sx,tx,tx)+CX(tx,tx,tx,ty))%P)%=P;
        }write((Ans+P)%P),pc('\n');
    }
}
本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2021-07-22 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 义乌中学暑假集训 2021.07.12 C
    • Description
      • Solution
        • Code
        领券
        问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档