zoj 3511 Cake Robbery(乱搞)

28 篇文章 0 订阅

题意:给出一个n凸多边形,现在切m刀,保证m刀在多边形内部没有交点,问最后剩下边数最多的多边形的边数。

思路:写完看了看网上的题解,有暴力模拟的,还有用线段树的Orz……自己写的基本是乱搞,把这个凸多边形看成一个图,那么我首先要找一个度数为2的点开始删除其所在的多边形,直到最后都删除掉。具体方法是这样的,对于一个度数为2的点,分别向两边扩展,直到找到度数不为2的两个点,那么如果这两个点之间被切了,也就是有边,那么这一块多边形可以直接拿掉,如果没有边怎么办呢?可以把找到的那一部分删掉,然后在这两个点之间连一条长度为找到的路径长度的边,这样就将规模缩减了,然后继续去找度数为2的点……实现的时候想到了之前学过DLX的删除操作,其实就是个双向循环链表,L[u],R[u]分别表示u左边和右边的点,如果要删掉v1,v2之间的路径,只需要L[v1]=v2;R[v2]=v1;就行了。边是用map标记的,虽然麻烦了点,但感觉解法还是挺有趣的~


代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
const int maxn=10000+10;
int L[maxn],R[maxn],degree[maxn],n,m;
map<pair<int,int> ,int>mp;
int solve()
{
    queue<int>q;
    for(int i=1;i<=n;++i)
        if(degree[i]==2) q.push(i);
    int u,maxv=0;
    pair<int,int>pii;
    while(!q.empty())
    {
        u=q.front();q.pop();
        if(degree[u]!=2) continue;
        degree[u]=0;
        int cnt1=0,cnt2=0,v1,v2;
        v1=R[u];
        cnt1=mp[make_pair(u,v1)];
        while(v1!=u)
        {
            if(degree[v1]!=2) break;
            degree[v1]=0;
            cnt1+=mp[make_pair(v1,R[v1])];
            v1=R[v1];
        }
        if(v1==u)
        {
            maxv=max(maxv,cnt1+mp[make_pair(v1,u)]);
            break;
        }
        v2=L[u];
        cnt2=mp[make_pair(u,v2)];
        while(v2!=u)
        {
            if(degree[v2]!=2) break;
            degree[v2]=0;
            cnt2+=mp[make_pair(v2,L[v2])];
            v2=L[v2];
        }
        int val=mp[make_pair(v1,v2)];
        if(val==0)
        {
            mp[make_pair(v1,v2)]=mp[make_pair(v2,v1)]=cnt1+cnt2;
            L[v1]=v2;R[v2]=v1;
        }
        else
        {
            maxv=max(maxv,cnt1+cnt2+1);
            mp[make_pair(v1,v2)]=mp[make_pair(v2,v1)]=1;
            L[v1]=v2;R[v2]=v1;
            degree[v1]--;degree[v2]--;
            if(degree[v1]==2) q.push(v1);
            if(degree[v2]==2) q.push(v2);
        }

    }
    return maxv;
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    while(~scanf("%d%d",&n,&m))
    {
        memset(degree,0,sizeof(degree));
        mp.clear();
        int x,y;
        for(int i=0;i<m;++i)
        {
            scanf("%d%d",&x,&y);
            mp[make_pair(x,y)]=1;
            mp[make_pair(y,x)]=1;
            degree[x]++;degree[y]++;
        }
        for(int i=1;i<=n;++i)
        {
            if(i<n) mp[make_pair(i,i+1)]=mp[make_pair(i+1,i)]=1;
            degree[i]+=2;
            L[i]=i-1;
            R[i]=i+1;
        }
        mp[make_pair(n,1)]=mp[make_pair(1,n)]=1;
        L[1]=n;R[n]=1;
        int ans=solve();
        printf("%d\n",ans);
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值