QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#133742#4942. Robust DefenseSolitaryDream#WA 511ms7072kbC++204.4kb2023-08-02 13:49:052023-08-02 13:49:08

Judging History

你现在查看的是最新测评结果

  • [2023-08-10 23:21:45]
  • System Update: QOJ starts to keep a history of the judgings of all the submissions.
  • [2023-08-02 13:49:08]
  • 评测
  • 测评结果:WA
  • 用时:511ms
  • 内存:7072kb
  • [2023-08-02 13:49:05]
  • 提交

answer

#include<bits/stdc++.h>
using namespace std;

#define int long long

using tii=tuple<int,int,int>;

const int N=511,M=998244353;

int qpow(int a,int b)
{
    int ret=1;
    while(b)
    {
        if(b&1)
            ret=1ll*ret*a%M;
        b>>=1;
        a=1ll*a*a%M;
    }
    return ret;
}

int n,m,S;

int f[N][5],val[N][4];

struct P{
    int x,y;
}a[N],b[N];

P operator +(const P &a,const P &b)
{
    return {a.x+b.x,a.y+b.y};
}

P operator -(const P &a,const P &b)
{
    return {a.x-b.x,a.y-b.y};
}

int operator *(const P &a,const P &b)
{
    return a.x*b.y-a.y*b.x;
}

int sgn(const P &a)
{
    if(a.x!=0)
        return a.x>0?-1:1;
    return a.y>0?-1:1;
}

bool operator <(const P &a,const P &b)
{
    if(sgn(a)!=sgn(b))
        return sgn(a)<sgn(b);
    return a*b>0;
}

vector<tii>E;

int ty(const P &a)
{
    if(a.x>0)
        return 1;
    if(a.x<0)
        return 3;
    if(a.y>0)
        return 2;
    return 4;
}

signed main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cin>>n>>m>>S;
    S=S*qpow(100,M-2)%M;
    for(int i=1;i<=n;i++)
        cin>>a[i].x>>a[i].y;
    for(int i=1;i<=m;i++)
        cin>>b[i].x>>b[i].y;
    sort(b+1,b+m+1,[&](const P &a,const P &b){
        return a.x<b.x||(a.x==b.x&&a.y<b.y);
    });
    for(int i=1;i<=m;i++)
        for(int j=1;j<=m;j++)
        {
            if(i==j)
                continue;
            P v=b[j]-b[i];
            bool ok=1;
            for(int k=1;k<=n;k++)
            {
                if(v*(a[k]-b[i])<0)
                    ok=0;
            }
            if(!ok)
                continue;
            int p=1;
            int lx=b[i].x,rx=b[j].x;
            if(lx<=rx)
                rx--;
            else
                rx++,swap(lx,rx);
            if(b[i].x!=b[j].x)
            {
                for(int k=1;k<=m;k++)
                {
                    if(k==i||k==j)
                        continue;
                    if(b[k].x<lx||b[k].x>rx)
                        continue;
                    if(v*(b[k]-b[i])>0)
                        continue;
                    p=p*(M+1-S)%M;
                }
            }
            else
            {
                int ly=b[i].x,ry=b[j].x;
                if(ry<=ry)
                    ry--;
                else
                    ry++,swap(ly,ry);
                for(int k=1;k<=m;k++)
                {
                    if(k==i||k==j)
                        continue;
                    if(b[k].y<ly||b[k].y>ry)
                        continue;
                    if(v*(b[k]-b[i])>0)
                        continue;
                    p=p*(M+1-S)%M;
                }

            }
            E.push_back({i,j,p});
        }
    sort(E.begin(),E.end(), [&](const tii &x,const tii &y) {
        auto [ai,aj,ap]=x;
        auto [bi,bj,bp]=y;
        auto da=b[aj]-b[ai],db=b[bj]-b[bi];
        if(da<db||db<da)
            return da<db;
        if(sgn(da)<0)
            return ai<bi;
        else
            return ai>bi;
    });
    for(int i=1;i<=m;i++)
    {
        for(int j=0;j<5;j++)
            val[i][j]=1;
        for(int j=1;j<=m;j++)
        {
            if(b[j].x<b[i].x&&b[j].y<=b[i].y)
                val[i][0]=(val[i][0]*(M+1-S))%M;
            if(b[j].x>=b[i].x&&b[j].y<b[i].y)
                val[i][1]=(val[i][1]*(M+1-S))%M;
            if(b[j].x>b[i].x&&b[j].y>=b[i].y)
                val[i][2]=(val[i][2]*(M+1-S))%M;
            if(b[j].x<=b[i].x&&b[j].y>b[i].y)
                val[i][3]=(val[i][3]*(M+1-S))%M;
        }
    }
    int ans=0;
    for(int i=1;i<=m;i++)
    {
        for(int t=0;t<5;t++)
            for(int j=1;j<=m;j++)
                f[j][t]=(t==0&&j==i);
        // cerr<<endl;
        for(auto [u,v,p]:E)
        {
            if(u<i||v<i)
                continue;
            int r=ty(b[v]-b[u]);
            // cerr<<r<<endl;
            for(int k=0;k<=r;k++)
            {
                int coe=1;
                for(int t=k;t<r;t++)
                    coe=coe*val[u][t]%M;
                f[v][r]=(f[v][r]+f[u][k]*coe%M*p%M*S)%M;
            }
            // f[v]=(f[v]+f[u]*p%M*S)%M;
            // f[v]=(f[v]+f[u]%M)%M;
        }
        // cerr<<b[i].x<<" "<<b[i].y<<" "<<f[i]<<endl;
        for(int j=3;j<5;j++)
            ans=(ans+f[i][j])%M;
    }
    cout<<ans<<endl;
}

詳細信息

Test #1:

score: 100
Accepted
time: 1ms
memory: 3460kb

input:

1 4 50
0 0
-1 0
3 0
0 1
2 -1

output:

686292993

result:

ok single line: '686292993'

Test #2:

score: 0
Accepted
time: 1ms
memory: 3464kb

input:

3 5 20
3 0
1 3
5 3
0 0
0 6
6 0
6 6
3 3

output:

771443236

result:

ok single line: '771443236'

Test #3:

score: 0
Accepted
time: 1ms
memory: 3496kb

input:

1 2 3
4 5
7 9
-2 -3

output:

184375732

result:

ok single line: '184375732'

Test #4:

score: -100
Wrong Answer
time: 511ms
memory: 7072kb

input:

500 500 47
7 19
16 17
20 13
1 10
17 9
5 23
12 2
15 12
16 8
11 8
8 12
3 2
11 13
23 0
3 23
13 10
9 12
11 5
8 18
6 0
6 20
3 9
1 21
13 18
5 11
9 15
8 17
6 18
1 8
4 24
7 14
11 11
2 9
8 9
23 3
17 15
21 10
19 7
13 16
0 10
0 7
6 17
11 9
9 4
1 15
21 12
1 24
20 7
21 7
20 0
10 3
3 24
2 12
18 11
20 5
14 20
10 4...

output:

630218927

result:

wrong answer 1st lines differ - expected: '963504722', found: '630218927'