UVa12633-Super Rooks on Chessboard-容斥+FFT

题目大意就是给你一个R*C的棋盘，上面有超级兵，这种超级兵会攻击 同一行、同一列、同一主对角线的所有元素，现在给你N个超级兵的坐标，需要你求出有多少方块是不能被攻击到的(R,C,N<50000)

遇到这种计数问题就要联想到容斥（组合数学太重要了），由容斥原理：

被攻击的方块数=行被攻击的方块数+列被攻击的方块数+主对角线被攻击的方块数-同时被行、列攻击的方块数-同时被行、对角线攻击的方块数-同时被列、对角线攻击的方块数+同时被行、列、对角线攻击的方块数

因为行列都不会很大，所以我们用三个数组分别记录行、列、对角线上含有超级兵的情况（从1开始）。对于对角线，我们不妨从右上角开始编起，从1-R+C-1,那么对于任意一块<x,y>，它所属的对角线为x-y+C。

开始计数：

行被攻击的方块数和列被攻击的方块数，这个很简单，就是扫一遍数组，然后如果有超级兵就加上

主对角线被攻击的方块数：我们需要知道对于任何一条主对角线，在R*C的棋盘里，它的主对角线有多长，根据推导，我们发现当对角线编号为1~C的时候，对角线都是从第一行开始的，它的长度按道理来讲应该是随着编号的增加递增的，但是当长度达到一定的时候就会收到行数的限制，因此长度为min(i,R)。可以理解为，i就是因为收到列数的限制的最大的长度，但是同时还需要收到行数的限制。当编号为C+1~R+C-1的时候，此时长度总体来讲会递减，主要收到行数的限制，但是同样也会收到列数的限制，因此长度为min(R+C-i,C)。得到长度以后我们就可以遍历一遍得到主对角线被攻击的方块数

同时被行、列攻击的方块数：这个很简单，就是总共被攻击的行数*被攻击的列数

同时被行、对角线攻击的方块数：我们只要知道对角线在行上的范围，然后再统计在这个范围内有多少行被攻击就可以了。为了快速得到一个范围内有多少行被攻击，我们可以用一个前缀数组，来计算某一个行区间内有多少被攻击。现在重点在于如何求对角线所占行的范围。同样的我们需要进行分类讨论：当序号为1~C的时候，每一个对角线都是从第1行开始的，我们只需要加上上面求得的长度就能够得到最后一行应该为min(i,C)。当序号为C+1~R+C-1的时候，每一行是从i-C+1开始的，到i-C+1+min(R+C-i,C)结束。

同时被列、对角线攻击的方块数：分析同上，我们这里只讨论每一条对角线列的范围额：当序号为1~C的时候，是C-i+1到C-i+1+min(i,C)，当序号为C+1~R+C-1的时候，是1到min(R+C-1,C)

同时被行、列、对角线攻击的方块数：这个可能就没有很简单了，我们需要求出所有既被行、列攻击，又被对角线攻击的块数。设X行Y列，则如果要同时被攻击，应该还存在一个对角线X-Y+C，观察这个式子，我们有什么方法迅速得到吗？好多个数相加、好多个数相减，我们还需要得到所有的结果，而且还不能遍历（会超时），我们就应该联想到FFT，加减运算转化为多项式乘法运算后，使用FFT就能nlogn解决问题。那么这两个多项式应该怎么构造呢？行多项式我们不妨就设为r[i]*xi,那么列多项式应该为c[i]*xC-i，这样运算以后的结果就是xi-j+C，得到的对角线的系数就是这条对角线同时被行列攻击的块数（因为i,j是不同的，所以是这个对角线上的不同的块）。如果i-j+C处恰好有一条对角线，接说明这些块需要加上。

可以看一哈我写的代码，应该还是比较清晰的。

AC代码：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<climits>
#include<algorithm>
#include<ctime>
#include<cstdlib>
#include<queue>
#include<set>
#include<map>
#include<cmath>using namespace std;
const int MAXN=1<<17;
const double PI=acos(-1.0);
typedef long long ll;int read()
{int x=0,sign=1; char c=getchar();while(c<'0' || c>'0') {if(c=='-') sign=-1; c=getchar();}while(c>='0' && c<='9'){x=x*10+c-'0'; c=getchar();}return x*sign;
}struct complex
{double r,i;complex(double _r=0,double _i=0):r(_r),i(_i){}complex operator +(const complex &b) {return complex(r+b.r,i+b.i);}complex operator -(const complex &b) {return complex(r-b.r,i-b.i);}complex operator *(const complex &b) {return complex(r*b.r-i*b.i,r*b.i+i*b.r);}
}A[MAXN],B[MAXN];void change(complex y[],int len)
{int i,j,k;for(i = 1, j = len/2;i < len-1;i++){if(i < j)swap(y[i],y[j]);k = len/2;while( j >= k){j -= k;k /= 2;}if(j < k)j += k;}
}void fft(complex y[],int len,int on)
{change(y,len);for(int h = 2;h <= len;h <<= 1){complex wn(cos(-on*2*PI/h),sin(-on*2*PI/h));for(int j = 0;j < len;j += h){complex w(1,0);for(int k = j;k < j+h/2;k++){complex u = y[k];complex t = w*y[k+h/2];y[k] = u+t;y[k+h/2] = u-t;w = w*wn;}}}if(on == -1)for(int i = 0;i < len;i++)y[i].r /= len;
}void FFT(int a[],int la,int b[],int lb)//la,lb分别是a,b数组的最高位+1 
{int len=1; while(len<la+lb) len<<=1;for(int i=0;i<la;++i) A[i]=complex(a[i],0);for(int i=la;i<len;++i) A[i]=complex(0,0);for(int i=0;i<lb;++i) B[i]=complex(b[i],0);for(int i=lb;i<len;++i) B[i]=complex(0,0);fft(A,len,1); fft(B,len,1);for(int i=0;i<len;++i) A[i]=A[i]*B[i];fft(A,len,-1);
}int R,C,D,N;
int r[MAXN],c[MAXN],d[MAXN];
int sumr[MAXN],sumc[MAXN];//行列前缀和 
ll ans;int main()
{//freopen("data.txt","w",stdout);int T;//T=read();scanf("%d",&T);for(int Case=1;Case<=T;++Case){//initans=0; sumr[0]=sumc[0]=0;for(int i=0;i<=R;i++) r[i]=0;for(int i=0;i<=C;i++) c[i]=0;for(int i=0;i<=D;i++) d[i]=0;//read //R=read(); C=read(); N=read();scanf("%d%d%d",&R,&C,&N);int u,v;while(N--){//u=read(); v=read();scanf("%d%d",&u,&v);r[u]=1; c[v]=1; d[u-v+C]=1;}//统计行列被攻击的块数,得到前缀和 for(int i=1;i<=R;++i){if(r[i]) ans+=C; sumr[i]=sumr[i-1]+r[i];}for(int i=1;i<=C;++i){if(c[i]) ans+=R; sumc[i]=sumc[i-1]+c[i];}//test//减去行列同时被攻击的块数ans-=sumr[R]*sumc[C]; //统计主对角线被攻击的块数//减去同时被行和对角线、列和对角线攻击的块数 //1~R-1+C  X-Y+CD=R-1+C;for(int i=1;i<=D;++i){if(d[i]){if(i<=C)//对角线块数递增，行数限制 {int dr1=0,dr2=min(i,R);	//行数始末 int dc1=C-i,dc2=dc1+dr2;//列数始末 ans+=dr2;ans-=sumr[dr2]-sumr[dr1];ans-=sumc[dc2]-sumc[dc1];}else//对角线块数递减，列数限制 {int dc1=0,dc2=min(R+C-i,C);//列数始末 int dr1=i-C,dr2=dr1+dc2;	//行数始末ans+=min(R+C-i,C);ans-=sumr[dr2]-sumr[dr1];ans-=sumc[dc2]-sumc[dc1];}}}//加上同时被行列对角线攻击的方块reverse(c,c+C+1);FFT(r,R+1,c,C+1);for(int i=1;i<=D;++i){if(d[i]){ans+=(ll)(A[i].r+0.5);	}}printf("Case %d: %lld\n",Case,(ll)R*C-ans);//if(Case!=T) printf("\n");}return 0;
}

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.mzph.cn/news/383760.shtml

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈email:809451989@qq.com，一经查实，立即删除！

相关文章

Linux网络编程---I/O复用模型之poll

Linux网络编程---I/O复用模型之poll

https://blog.csdn.net/men_wen/article/details/53456474Linux网络编程—I/O复用模型之poll 1.函数poll poll系统调用和select类似&#xff0c;也是在指定时间内轮询一定数量的文件描述符&#xff0c;以测试其中是否有就绪者。 #include <poll.h>int poll(struct pollfd…
阅读更多...
FFT模板

FFT模板

整理了一下&#xff0c;自己写了一下模板 const double PIacos(-1.0); struct complex {double r,i;complex(double _r0,double _i0):r(_r),i(_i){}complex operator (const complex &b) {return complex(rb.r,ib.i);}complex operator -(const complex &b) {return c…
阅读更多...
Linux网络编程---I/O复用模型之epoll

Linux网络编程---I/O复用模型之epoll

https://blog.csdn.net/men_wen/article/details/53456474 Linux网络编程—I/O复用模型之epoll 1. epoll模型简介 epoll是Linux多路服用IO接口select/poll的加强版&#xff0c;e对应的英文单词就是enhancement&#xff0c;中文翻译为增强&#xff0c;加强&#xff0c;提高&…
阅读更多...
POJ 1741tree-点分治入门

POJ 1741tree-点分治入门

学习了一下点分治&#xff0c;如果理解有误还请不吝赐教。 为了快速求得树上任意两点之间距离满足某种关系的点对数&#xff0c;我们需要用到这种算法。 点分治是树上的一种分治算法&#xff0c;依靠树和子树之间的关系进行分治从而降低复杂度。 和其他树上的算法有一些区别…
阅读更多...
基于单链表的生产者消费者问题

基于单链表的生产者消费者问题

『生产者与消费者问题分析』「原理」生产者生产产品&#xff0c;消费者消费产品。产品如果被消费者消费完了&#xff0c;同时生产者又没有生产出产品&#xff0c;消费者 就必须等待。同样的&#xff0c;如果生产者生产了产品&#xff0c;而消费者没有去消费&#x…
阅读更多...
C++智能指针（一）智能指针的简单介绍

C++智能指针（一）智能指针的简单介绍

https://blog.csdn.net/nou_camp/article/details/70176949C智能指针 在正式了解智能指针前先看一下下面的一段代码 #include<iostream> using namespace std; class A { public:A():_ptr(NULL), _a(0){}~A(){} public:int* _ptr;int _a; };void test() {A a;int *p1 ne…
阅读更多...
聪聪可可-点分治

聪聪可可-点分治

聪聪和可可是兄弟俩&#xff0c;他们俩经常为了一些琐事打起来&#xff0c;例如家中只剩下最后一根冰棍而两人都想吃、两个人都想玩儿电脑&#xff08;可是他们家只有一台电脑&#xff09;……遇到这种问题&#xff0c;一般情况下石头剪刀布就好了&#xff0c;可是他们已经玩儿…
阅读更多...
C++智能指针（二）模拟实现三种智能指针

C++智能指针（二）模拟实现三种智能指针

https://blog.csdn.net/nou_camp/article/details/70186721在上一篇博客中提到了Auto_ptr(C智能指针&#xff08;一&#xff09;)&#xff0c;下面进行模拟实现Auto_ptr 采用类模板实现 #include<iostream> using namespace std; template<class T> class Autoptr …
阅读更多...
Prime Distance On Tree-树分治+FFT

Prime Distance On Tree-树分治+FFT

题目描述 Problem description. You are given a tree. If we select 2 distinct nodes uniformly at random, what’s the probability that the distance between these 2 nodes is a prime number? Input The first line contains a number N: the number of nodes in this…
阅读更多...
C++智能指针（三）总结

C++智能指针（三）总结

https://blog.csdn.net/nou_camp/article/details/70195795 在上一篇博客中&#xff08;C智能指针&#xff08;二&#xff09;&#xff09;模拟实现了三种智能指针。 其中最好的就是shared_ptr,但是这并不代表它就是最完美的&#xff0c;它也有问题&#xff0c;这个问题就是循环…
阅读更多...
POJ2114-Boatherds-树分治

POJ2114-Boatherds-树分治

题目描述 Boatherds Inc. is a sailing company operating in the country of Trabantustan and offering boat trips on Trabantian rivers. All the rivers originate somewhere in the mountains and on their way down to the lowlands they gradually join and finally th…
阅读更多...
c++11 你需要知道这些就够了

c++11 你需要知道这些就够了

https://blog.csdn.net/tangliguantou/article/details/50549751c11新特性举着火把寻找电灯今天我就权当抛砖引玉&#xff0c;如有不解大家一起探讨。有部分内容是引用自互联网上的内容&#xff0c;如有问题请联系我。T&& 右值引用 std::move 右值引用出现之前我们只能…
阅读更多...
HDU5977-Garden of Eden-树分治+FWT

HDU5977-Garden of Eden-树分治+FWT

题目描述 When God made the first man, he put him on a beautiful garden, the Garden of Eden. Here Adam lived with all animals. God gave Adam eternal life. But Adam was lonely in the garden, so God made Eve. When Adam was asleep one night, God took a rib fro…
阅读更多...
C++11新特性学习

C++11新特性学习

https://blog.csdn.net/tennysonsky/article/details/778170481、什么是C11C11标准为C编程语言的第三个官方标准&#xff0c;正式名叫ISO/IEC 14882:2011 - Information technology -- Programming languages -- C。在正式标准发布前&#xff0c;原名C0x。它将取代C标准第二版I…
阅读更多...
C++ override 关键字用法

C++ override 关键字用法

override关键字作用&#xff1a; 如果派生类在虚函数声明时使用了override描述符&#xff0c;那么该函数必须重载其基类中的同名函数&#xff0c;否则代码将无法通过编译。举例子说明struct Base {virtual void Turing() 0;virtual void Dijkstra() 0;virtual void VNeumann…
阅读更多...
Gym - 101981I-MagicPotion-最大流

Gym - 101981I-MagicPotion-最大流

题目描述 There are n heroes and m monsters living in an island. The monsters became very vicious these days, so the heroes decided to diminish the monsters in the island. However, the i-th hero can only kill one monster belonging to the set Mi. Joe, the st…
阅读更多...
c++仿函数 functor

c++仿函数 functor

https://www.cnblogs.com/decade-dnbc66/p/5347088.html内容整理自国外C教材先考虑一个简单的例子&#xff1a;假设有一个vector<string>&#xff0c;你的任务是统计长度小于5的string的个数&#xff0c;如果使用count_if函数的话&#xff0c;你的代码可能长成这样&#…
阅读更多...
HDU4812-D Tree-树分治

HDU4812-D Tree-树分治

题目描述 There is a skyscraping tree standing on the playground of Nanjing University of Science and Technology. On each branch of the tree is an integer (The tree can be treated as a connected graph with N vertices, while each branch can be treated as a v…
阅读更多...
成为C++高手之实战项目

成为C++高手之实战项目

https://blog.csdn.net/niu_gao/article/details/51458721 在内存中模拟出一副牌&#xff0c;然后模拟洗牌&#xff0c;发牌等动作。 流程是这样的&#xff1a;构建一副牌保存到一个数组中—洗牌—创建玩家—向玩家发牌–输出每个玩家的牌。 #include <stdio.h> #include…
阅读更多...
C++中String类的实现

C++中String类的实现

https://www.cnblogs.com/zhizhan/p/4876093.html原文&#xff1a;http://noalgo.info/382.html String是C中的重要类型&#xff0c;程序员在C面试中经常会遇到关于String的细节问题&#xff0c;甚至要求当场实现这个类。只是由于时间关系&#xff0c;可能只要求实现构造函数、…
阅读更多...
最新文章

Copyright @ 2022~2023